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Abstract 


The  radio  frequency  (RF)  spectrum  is  a  limited  resource.  Spectrum  allotment  disputes 
stem  from  this  scarcity  as  many  radio  devices  are  confined  to  a  fixed  frequency  or 
frequency  sequence.  One  alternative  is  to  incorporate  cognition  within  a  reconfigurable 
radio  platform,  therefore  enabling  the  radio  to  adapt  to  dynamic  RF  spectrum  environments. 
In  this  way,  the  radio  is  able  to  actively  observe  the  RF  spectrum,  orient  itself  to  the  current 
RF  environment,  decide  on  a  mode  of  operation,  and  act  accordingly,  thereby  sharing  the 
spectrum  and  operating  in  more  flexible  manner.  This  research  presents  a  novel  framework 
for  incorporating  several  techniques  for  the  purpose  of  adapting  radio  operation  to  the 
current  RF  spectrum  environment.  Specifically,  this  research  makes  six  contributions  to  the 
field  of  cognitive  radio:  (1)  the  framework  for  a  new  hybrid  hardware/software  middleware 
architecture,  (2)  a  framework  for  testing  and  evaluating  clustering  algorithms  in  the  context 
of  cognitive  radio  networks,  (3)  a  new  RF  spectrum  map  representation  technique,  (4) 
a  new  RF  spectrum  map  merging  technique,  (5)  a  new  method  for  generating  a  random 
key-based  adaptive  frequency-hopping  waveform,  and  (6)  initial  integration  testing  toward 
implementing  the  proposed  system  on  a  field-programmable  gate  array  (FPGA). 
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AN  ARCHITECTURE  FOR  COEXISTENCE  WITH  MULTIPLE  USERS  IN 


FREQUENCY  HOPPING  COGNITIVE  RADIO  NETWORKS 

1  Introduction 

Today’s  radio  frequency  (RF)  spectrum  is  increasingly  congested.  In  America, 
the  Federal  Communications  Commission  (FCC)  partitions  the  RF  spectrum  into 
frequency  ranges  based  on  the  needs  and  capabilities  of  radio  equipment.  Originally, 
the  FCC  performed  this  partitioning  according  to  input  from  public  hearings.  Over  time, 
spectrum  allocation  became  the  product  of  systems  such  as  lotteries  and  spectrum  auctions. 
Currently,  the  pervasiveness  of  equipment  that  uses  the  RF  spectrum  (radio  and  television 
stations,  cell  phones,  wireless  internet,  etc.)  demands  more  flexibility  of  the  spectrum 
allocation  process.  In  the  future,  the  ever-expanding  wireless  footprint  implies  that  devices 
will  need  to  adapt  their  transmission  and  reception  capabilities  to  the  RF  spectrum  in  which 
they  communicate.  This  practice  is  known  as  “dynamic  spectrum  access”  (DSA).  However, 
even  with  an  adaptive  communication  technique,  the  RF  spectrum’s  density  (especially  the 
overlapping  of  spread- spectrum  transmissions)  will  still  likely  detract  from  an  adaptive 
waveform’s  transmission. 

Additionally,  some  spectrum  assignments  are  applied  to  all  times  during  the  day,  yet 
many  are  only  used  consistently  for  brief  periods  of  time.  Frequencies  which  may  be  used 
at  any  given  time  by  a  specific  user,  but  which  are  used  only  periodically,  significantly 
contribute  to  wasted  spectrum  assignment.  This  issue  is  extensible  to  the  military’s  need  for 
radios  which  can  function  regardless  of  their  geographic  location,  the  time  of  day,  and  the 
local  spectrum  policy.  Were  military  (or  civilian)  radios  able  to  operate  in  unused  spectrum 
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Figure  1.1:  The  OODA  loop  as  a  “cognitive  cycle”  [3]. 


during  times  of  limited  actual  usage,  the  spectrum  could  be  used  much  more  efficiently. 
This  research  presents  an  architecture  through  which  this  need  can  be  met. 

It  is  believed  a  cognitive  radio  (CR)  can  be  used  to  implement  DSA  in  solving  the 
aforementioned  questions  as  related  to  RF  spectrum  congestion  and  spatial/temporal  policy 
conflicts.  CRs  are  a  subset  of  the  now-common  radio  implementation,  the  software-defined 
radio.  The  CR’s  power  is  twofold.  First,  by  using  a  field-programmable  gate  array  (FPGA), 
the  CR  becomes  a  high-speed,  circuit-based  yet  reconfigurable  device  [1].  (For  this  reason, 
some  refer  to  CRs  as  “firmware-defined  radios”).  Second,  CRs  implement  Colonel  John 
Boyd’s  famed  Observe-Orient-Decide-Act  (“OODA”)  loop  as  a  means  by  which  to  make 
informed,  adaptive  decisions  on  transmission  and  reception  policy.  In  summary,  the  CR  is 
DSA’s  enabling  technology  [2]. 

Colonel  Boyd’s  OODA  loop  models  the  initial  problem  set.  First,  a  radio  senses  its 
local  RF  spectrum  (observation).  Next,  the  software  contained  on  the  radio  constructs  a 
map  of  the  RF  spectrum  (a  radio  environment  map,  or  REM)  from  which  to  use  available 
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frequencies  (orientation).  Third,  the  radio  assembles  a  randomized  frequency  hopping 
pattern  (the  hopset),  for  distribution  to  neighboring  radios  (decision).  Finally,  the  radio  uses 
the  hopset  for  adaptive  communication  (action).  There  are  various  namings,  interpretations, 
and  implementations  of  this  model,  such  as  that  of  [3]  (e.g.,  the  “cognitive  cycle”),  but  each 
distills  down  to  the  OODA  loop.  Figure  1 . 1  shows  how  the  OODA  loop  maps  directly  to  the 
interpretation  from  [3].  The  mapping  of  colors  to  OODA  loop  functions  remains  constant 
in  the  related  diagrams  throughout  this  document. 

Figure  1.2  shows  how  the  proposed  system  is  follows  the  OODA  model  over  time  in 
the  presence  of  interference.  It  is  worthwhile  to  note  that  this  system  is  not  focused  on  the 
transmission  or  reception  of  specific  data — only  adapting  to  the  current  RF  environment. 
At  period  A,  the  system  constantly  observes  the  environment.  Period  B  shows  notional 
expected  performance  while  the  spectrum  is  changing,  the  existing  adaptations  are  under 
duress,  and  the  radio  is  reorienting  itself.  Period  C  is  when  the  system  must  decide  how 
to  react  to  the  new  RF  environment.  The  system  (re)acts  using  its  decision  at  period  D 
and,  finally,  resumes  observation  at  some  improved  level  of  performance  (three  feasible 
outcomes  are  shown). 

The  natural  extension  of  this  process  is  implementing  it  in  a  cognitive  radio  network 
(i.e.,  two  or  more).  In  this  case,  the  observation  step  involving  spectrum  sensing  expands 
to  include  spectrum  observation  dispersion  within  the  network.  If  those  radios  which  must 
communicate  are  aware  of  the  RF  spectra  at  those  other  nodes  within  the  network,  then 
observation  is  complete  within  the  network. 

1.1  Motivation 

1.1.1  Air  Force  Technology  Horizons  2010. 

The  vision  in  [4]  delivers  a  two-decade  (2010-2030)  outlook  of  the  Air  Force’s  current 
technological  capabilities,  an  assessment  of  future  challenges,  and  a  plan  for  the  molding 
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Figure  1.2:  Expected  system  operation. 


of  Air  Force  technology  to  meet  such  challenges.  Produced  by  the  office  of  the  Air  Force 
Chief  Scientist,  it  provides  a  high  level  review  to  the  service’s  senior  staff,  specifically  the 
Secretary  of  the  Air  Force  and  Chief  of  Staff  of  the  Air  Force.  “Technology  Horizons” 
presents  12  Air  Force  “Service  Core  Functions”  (SCFs),  analogous  to  “strategic”-level 
functions  needed  by  the  service  to  accomplish  its  missions.  Also  identified  are  30  “Potential 
Capability  Areas”  (PCAs).  PCAs  are  defined  as  “credibly  achievable  within  the  time 
horizon  addressed  by  this  study.”  Additionally,  each  PCA  is  mapped  to  each  SCF  such 
that  every  SCF  is  supported  by  multiple  PCAs.  Under  the  PCAs,  and  at  a  finer  level  of 
granularity,  the  document  also  lists  110  “Key  Technology  Areas”  (KTAs).  KTAs  compose 
a  large,  overlapping  set  of  those  technologies  that  make  each  PCA  possible.  In  paraphrase, 
these  are  the  focus  areas  the  Air  Force  believes  it  can  and  will  technologically  conquer  by 
the  year  2030  [4]. 

PCA  number  seven  is  “Frequency  Agile  Spectrum  Utilization”  (FASU).  Within  FASU 
is  Dynamic  Spectrum  Utilization  (DSA),  the  core  tactical  goal  of  this  thesis.  The  fact  that 
the  Air  Force  has  identified  FASU,  and  by  extension,  DSA,  as  a  major  tenet  of  the  service’s 
near-,  mid-,  and  long-term  technical  superiority  calls  for  additional  research,  development, 
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Figure  1.3:  Conflicting  LightSquared  and  GPS  spectrum  assignments. 


and  fielding  of  systems  in  this  key  area.  In  the  context  of  [4],  FASU  is  featured  in  five 
PCAs.  Of  larger  impact  is  the  high  number  of  SCFs  to  which  the  FASU  PCA  maps:  PCA 
number  seven  maps  to  11  out  of  12  SCFs;  of  all  30  PCAs,  only  five  map  to  an  equal  or 
greater  amount  of  SCFs.  The  bottom  line  is  this:  FASU  and  DSA  are  projected  to  have  a 
significant  impact  on  Air  Force  strategy  and  operations  over  the  next  two  decades. 

1.1.2  LightSquared  versus  the  FCC. 

Spectrum  access  and  sharing  issues  are  also  of  relevance  in  the  civilian  sector. 
In  recent  years,  the  company  LightSquared  has  sought  to  build  a  nation-wide  wireless 
broadband  network  using  previously  vacant  spectrum  (Figure  1.3,  darkened  within  the 
Mobile  Satellite  Band)  near  the  GPS  band  (Figure  1.3,  darkened  within  the  Global  Satellite 
Navigation  Band).  In  fact,  the  effort  was  originally  sanctioned  (and  even  propelled)  by 
the  FCC.  It  was  eventually  realized  through  testing,  though,  that  LightSquared’s  spectrum 
allotment  and  the  energy  generated  within  that  spectrum  would  bleed  over  into  the  global 
satellite  navigation  band.  Specifically,  the  GPS  sub-band  would  be  impacted,  if  not 
drowned  out,  in  certain  locations  [5]. 

As  a  result  of  this  conflict,  the  FCC  directed  LightSquared  to  form  a  working  group 
comprised  of  LightSquared,  government,  and  industry  GPS  engineers.  The  working 
group’s  decision  was  issued  in  a  report  in  June  2011  with  the  declaration  that  the  current 
course  of  action  would  indeed  degrade  GPS  receiver  performance.  More  damning  was  a 
comment  from  the  group’s  aviation  committee:  “For  the  originally  defined  LightSquared 
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spectrum  deployment  scenarios,  GPS-based  operations  are  expected  to  be  unavailable  over 
entire  regions  of  the  country  at  any  normal  operational  aircraft  altitude'’.  In  other  words, 
in  many  areas  aircraft  flying  below  2,000  feet  above  ground  level  would  have  little  or  no 
GPS  guidance.  This  is  particularly  dangerous  because  many  aircraft,  from  general  aviation 
to  airliners  to  military  aircraft,  use  a  high-precision  version  of  GPS  to  take  off  and  land. 
Degraded  GPS  in  these  environments  can  be  deadly,  especially  if  it  is  unexpected. 

The  obvious  source  of  contention  is  spectrum  allocation  policy  as  it  pertains  to 
governing  transmitters.  However,  [5]  poses  the  following  question:  Is  it  time  for  the 
FCC  to  focus  on  protecting  receivers ?  Unfortunately,  the  LightSquared/FCC  debate  is 
too  advanced  redesign  the  underlying  hardware,  but  in  future  and  with  receiver-oriented 
regulations,  a  situation  such  as  the  LightSquared  incident  could  be  resolved  on  the  drawing 
board.  For  example,  if  GPS  receivers  were  protected,  the  issues  would  be  oriented 
toward  the  end-user  from  the  very  start  of  system  design — it  is  reasonable  to  assume 
the  LightSquared  issue  could  have  been  avoided  early  in  the  process  by  designing  around 
receiver  protections.  Although  there  are  no  FCC  regulations  in  place  to  protect  receivers, 
this  research  is  a  start  toward  solving  the  problem  of  dynamic  receiver  behavior  in  the 
presence  of  interference. 

The  LightSquared  incident  is  a  clear  example  of  the  need  for  DSA  and  FASU  on  a 
national  infrastructure  scale.  Furthermore,  CR  is  an  ideal  method  for  implementing  both 
concepts.  Given  a  CR’s  natural  ability  to  avoid  spectrum  conflicts,  it  becomes  a  highly 
attractive  solution  to  conflicts  ranging  from  low-profile,  low-impact  personal  networks  to 
high-profile,  high-impact  situations  such  as  the  predicament  of  LightSquared. 

1.1.3  Tactical  DSA. 

The  Defense  Advanced  Research  Projects  Agency  (DARPA)  has  also  identified  the 
need  for  DSA  in  tactical  scenarios.  Currently,  ground-based  troops  have  no  means  for 
detecting  spectrum  usage  for  the  purpose  of  avoiding  such  interference.  In  theater,  troops 
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are  often  faced  with  spectrum  congestion.  In  a  tactical  environment  where  a  mission  is  at 
stake  and  lives  are  at  risk,  arguably  the  worst  problem  for  a  commander  is  a  breakdown 
or  total  loss  of  communication  with  his  troops.  DARPA  seeks  to  alleviate  this  issue 
by  providing  commanders  and  their  troops  with  a  means  for  identifying  occupied  and 
unoccupied  RF  spectrum.  This  capability  does  not  focus  on  the  content  of  transmissions, 
but  solely  on  classifying  the  energy  present  in  the  local  RF  spectrum  [6,  7]. 

Given  a  digital  analysis  of  the  spectrum,  this  research  presents  a  method  in  which 
a  map  of  the  radio  spectrum  intelligible  by  frequency-hopping  hardware  is  produced.  The 
representation  is  not  tailored  for  visual  inspection  or  human  interface  as  would  be  DARPA’ s 
proposed  map.  It  does,  however,  give  the  components  implemented  in  this  research  the 
capability  to  identify  unoccupied  portions  of  the  RF  spectrum  for  the  purpose  of  avoiding 
unnecessary  congestion  and  interference. 

Additionally,  DARPA  previously  researched  a  new  method  for  taking  advantage  of 
unused  RF  spectrum-conceivably  for  use  with  such  a  map  primarily  for  public  DSA  [8]. 
The  neXt  Generation  (XG)  communication  program  was  proven  in  2006  in  a  six-node 
network  to  be  able  to  opportunistically  use  otherwise  wasted  spectrum  [9].  The  proposed 
system  has  the  same  basic  function  regarding  the  ability  to  use  open  spectrum  space.  In 
addition  to  building  on  the  XG  radio’s  functionality  and  network  size,  this  research  involves 
a  fully-scalable  network  (to  be  tested  in  this  work  at  100  nodes)  and  a  frequency-hopping 
waveform  in  place  of  a  single-frequency  operation. 

1.1.4  Cognitive  Radio  Network  Reliability. 

If  a  radio  is  truly  cognitive,  there  is  a  well-defined  set  of  parameters  by  which  the  radio 
evaluates  its  situation  and  reacts  in  a  way  that  enables  it  to  continue  operation  using  a  new 
configuration.  We  believe  reliability  in  a  cognitive  radio  network  is  the  result  of  effective 
resource  usage  after  it  adapts  to  its  environment.  Our  cognitive  radio  uses  a  frequency 
hopping  waveform  for  communication.  Therefore,  the  system’s  cognitive  functionality 
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stems  from  its  ability  to  evaluate  occupied  and  unoccupied  frequencies  within  a  network, 
form  an  adaptive  hopset,  and  use  the  selected  waveform.  By  implementing  an  adaptive 
process,  the  CR  becomes  a  direct  solution  to  the  questions  of  FASU  and  DSA  in  that 
is  reconfigurable  to  multiple  frequencies  and  can  access  the  spectrum  in  a  self-propelled 
manner. 

Reliability  in  our  system,  then,  is  a  function  of  the  number  of  frequencies  from  which 
the  radio  selects  its  hopset.  A  large  set  of  open  frequencies  allows  the  radio  to  traverse 
more  of  the  RF  spectrum,  thereby  decreasing  the  likelihood  of  significant  data  loss  due  to 
interference.  It  is  expected  geographically  partitioning  a  CRN  yields  groups,  or  clusters, 
with  similar  observed  RF  spectra.  This  experiment  investigates  which  clustering  heuristic, 
if  any,  partitions  the  network  such  that  overall  network  complexity  is  reduced  for  the 
purpose  of  minimizing  extraneous  network  traffic  when  transmitting  data  (i.e.,  spectrum 
maps,  waveforms,  actual  data,  etc.).  If  network  complexity  is  reduced,  it  is  expected 
the  network  becomes  more  reliable  as  data  is  transmitted  more  efficiently  and  with  fewer 
collisions  and/or  delays. 

1.2  Problem  Statement 

The  following  operations  govern  the  functionality  of  the  whole  system.  Brief 
operation  descriptions  are  given  here,  while  full  explanations  are  presented  in  Section 
3.1.2.  Each  operation  corresponds  to  a  state  of  the  same  label  in  Figure  1.4,  where 
“network  partitioning”  and  “waveform”  map  to  implementing  “clustering”  and  “hopset,” 
respectively,  in  this  research.  Additionally,  each  operation  maps  to  the  color-corresponding 
OODA  function  in  Figure  1.1.  This  coloring  scheme  is  continued  through  the  remainder  of 
the  document. 
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Figure  1 .4:  System  function. 


•  0-1.  Sense  RF  spectrum.  Detect  the  presence  of  other  spectrum  users. 

•  0-2.  Form  network  partitions.  Partition  the  network  into  sub-networks  based  on 
geographic  location. 

•  0-3.  Share  spectrum  maps.  Distribute  spectrum  knowledge  throughout  the 
network. 

•  0-4.  Merge  spectrum  maps.  Fuse  the  sensed  data  for  optimal  waveform  selection. 
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•  0-5.  Select  waveform.  Generate  a  waveform  common  to  all  nodes  within  a  sub¬ 


network. 

•  0-6.  Send  waveform.  Distribute  waveform  information  to  all  nodes  within  network. 

•  0-7.  Use  waveform.  Engage  in  normal  communication  using  the  adaptive 
waveform. 

The  issues  identified  in  Section  1.1  represent  a  sampling  of  the  issues  currently 
recognized  within  the  DSA  community.  It  is  expected  a  framework  that  accomplishes  the 
preceding  seven  steps  is  a  step  toward  answering  some  of  these  questions.  This  research 
proposes  such  a  framework  in  the  form  of  a  middleware  architecture  and  implements 
several  key  components.  Succinctly,  the  goal  is  to  answer  the  following  question:  How  can 
we  implement  a  frequency  hopping  cognitive  radio  network  for  coexistence  with  multiple 
users? 

Within  the  above  problem  statements,  several  terms  are  defined: 

•  Frequency  hopping.  Devices  must  be  able  to  use  a  frequency  hopping  waveform 
that  utilizes  multiple  frequencies  for  discrete  periods  of  time. 

•  Cognitive  radio  network.  Two  or  more  radios  which  can  operate  according  to  the 
OODA  loop,  specifically  the  architecture  proposed  in  this  work,  form  a  cognitive 
radio  network. 

•  Coexist.  Cognitive  radios  are  designed  for  the  purpose  of  avoiding  other  energy 
sources  in  the  RF  spectrum. 

•  Spectrum  users.  It  is  highly  likely  that  other  devices  similar  to  the  one  in  this 
research  also  are  capable  of  using  the  same  RF  spectrum  usable  by  the  proposed 
system. 
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1.3  Research  Contributions 


This  research  makes  six  contributions  to  the  field  of  cognitive  radio: 

1.  The  framework  for  a  new  hybrid  hardware/software  middleware  architecture  (see 
Figure  3.1); 

2.  A  framework  for  testing  and  evaluating  clustering  algorithms  in  the  context  of 
cognitive  radio  networks  (see  Figure  3.1,  0-2  and  Figure  3.5); 

3.  A  new  RF  spectrum  map  representation  technique  (see  Figure  3.6); 

4.  A  new  RF  spectrum  map  merging  technique  (see  Figure  3.1,  0-4  and  see  Figure 
3.22); 

5.  A  new  method  for  generating  a  random,  key-based  adaptive  hopset  frequency 
hopping  waveform  (see  Figure  3.1,  0-5  and  see  Figure  3.23);  and 

6.  Initial  integration  testing  toward  implementing  the  proposed  system  on  a  field- 
programmable  gate  array  (FPGA)  (see  Figure  4.14). 

It  is  expected  that  these  contributions  form  the  foundation  for  fulfilling  the  original 
architecture,  where  the  eventual  result  is  a  deployable  prototype  ad  hoc  network  of 
cognitive  radios  which  operate  in  the  presence  of  additional  spectrum  users. 

1.4  Thesis  Organization 

Related  existing  work  and  background  information  on  implemented  components  is 
summarized  in  the  next  chapter.  Components  implemented  in  this  project  are  validated  in 
the  Methodology  chapter  and  analyzed  in  the  Results  chapter.  Conclusions  are  presented 
in  the  fifth  chapter,  and  future  work  is  detailed  in  the  last  chapter.  Appendices  containing 
additional  figures  and  software  code  are  at  the  end  of  this  document. 
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2  Related  Work 


This  work  combines  the  results  of  previous  related  but  independent  endeavors,  specif¬ 
ically  spectrum  mapping  techniques,  network  clustering,  FPGA-based  cognitive  ra¬ 
dio  implementations,  and  adaptive  frequency  hopping  (AFH).  This  chapter  addresses  orig¬ 
inal  work  in  these  areas  and  this  work’s  contributions  to  each  area. 

2.1  Cognitive  Radio 

The  ever-crowded  RF  spectrum  demands  a  utility  by  which  spectrum  may  be  used 
more  efficiently  by  more  users  for  less  cost.  Primary  users  purchase  spectrum  bands  for 
unrestricted  usage.  Such  users  do  not  always  use  their  allotted  spectrum.  Secondary  users 
(SUs),  on  the  other  hand,  do  not  possess  ownership  of  any  RF  spectrum  bands  and  are 
relegated  to  the  unlicensed  bands  (Industry-Science-Medicine  (ISM),  Unlicensed  National 
Information  Infrastructure  (UNII),  etc.).  Given  the  permeation  of  radio-centric  devices 
(cellular  telephones,  wireless  internet,  Bluetooth,  Zigbee,  etc.)  into  typical  daily  functions, 
secondary  users  are  highly  prone  to  inadvertent  interference  from  other  secondary  users. 
This  realization  begs  the  question:  How  can  secondary  users  with  their  limited  spectrum 
assignments  efficiently  use  underutilized  spectrum  allocated  to  primary  users? 

Past  systems  designed  for  resilience  in  an  interference-  or  jamming-prone  environment 
paved  the  way  for  today’s  advanced  radio  systems.  For  example,  radios  such  as  the 
military’s  SINCGARS  and  HAVE  QUICK  hop  frequencies  based  on  a  user-supplied  key, 
but  the  hopping  sequence  is  independent  of  the  surrounding  spectrum.  Legacy  adaptive 
frequency  hopping  radios  date  back  to  the  early  1990s,  but  only  in  the  past  decade  has  the 
requisite  commercial  technology  made  such  systems  feasible  in  non-military  applications. 
The  need  for  pervasive  adaptable  radios  has  long  been  recognized,  and  the  supporting 
research  and  competencies  are  now  coming  to  fruition. 
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Enter  the  cognitive  radio.  A  progression  of  software-defined  radio  (SDR),  a  cognitive 
radio  boasts  the  technologies  necessary  to  use  large  portions  of  RF  spectrum  designated 
for  primary  users  as  the  transmission  medium  for  secondary  users.  One  of  the  largest 
attractions  to  cognitive  radios  is  that  they  are  designed  to  adapt  to  the  needs  of  spectrum 
regulators,  network  operators,  and  user  objectives  [10].  This  promise  has  driven  explosive 
growth  in  just  the  past  five  years.  For  example,  in  [10],  it  is  mentioned  that  an  internet 
search  on  the  term  “cognitive  radio”  in  2009  yields  138,000  hits-triple  the  number  in  2006. 
A  similar  search  today,  however,  gives  6,820,000  hits -nearly  50  times  the  number  four 
years  ago.  Clearly,  cognitive  radio  has  taken  hold. 

Cognitive  radio  is  defined  by  the  SDR  Forum  [11]  and  the  Institute  of  Electrical  and 
Electronics  Engineers  (IEEE)  1900.1-2008  [12]  working  group  as  follows: 

1 .  “Radio  in  which  communications  systems  are  aware  of  their  environment  and  internal 
state,  and  can  make  decisions  about  their  radio  operating  behavior  based  on  that 
information  and  predefined  objectives.  The  environmental  information  may  or  may 
not  include  location  information  related  to  communication  systems. 

2.  “Cognitive  radio  (as  defined  in  1)  that  uses  SDR,  adaptive  radio,  and  other 
technologies  to  automatically  adjust  its  behavior  or  operations  to  achieve  desired 
objectives.” 

When  fueled  by  popularity  among  engineers  and  businesses,  formal  adoption  within 
the  scientific  community,  and  maturation  of  key  technologies,  the  cognitive  radio  is  quickly 
becoming  the  means  to  achieve  what  is  needed:  a  smarter  radio.  This  work  follows  the 
above  definitions  for  building  and  implementing  a  cognitive  radio 

2.1.1  An  Increasingly  Smarter  Radio. 

A  popular  and  growing  field,  cognitive  radio  has  seen  steady  advancement,  most 
notably  in  the  past  five  years.  Efficient  spectrum  usage,  to  include  dynamic  allocation 
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techniques  and  sharing  schemes,  drives  most  of  this  research  for  the  purpose  of  allowing 
secondary  users  to  avoid  primary  users  [3].  In  their  journal  article,  the  authors  of  [3]  cover 
the  technological  advances  in  cognitive  radio.  These  advances  are  addressed  in  Tables 
2.1,  2.2,  2.3,  and  2.4  along  with  the  projected  level  of  contribution  of  this  research  to 
each  area  (Primary,  Secondary,  or  Not  Applicable  [N/A]).  Primary  areas  receive  a  direct 
contribution  from  this  work,  secondary  areas  are  those  on  which  this  work  is  formed  but  do 
not  receive  any  new  contribution  from  this  work,  and  irrelevant  areas  are  not  the  focus  of 
any  contribution  in  this  work.  An  argument  is  presented  for  each  for  each  irrelevant  item 
as  relevant  items  are  further  discussed  throughout  the  document. 

Each  fundamental  advance  is  relevant  to  this  research,  due  in  large  part  to  the 
broadness  of  each  area.  However,  spectrum  sensing  is  assumed  to  be  an  existing  function 
prior  to  implementing  the  proposed  architecture.  Cognitive  capability  (i.e.,  adaptive  hopset 
selection)  and  dynamic  network  reconfiguration  comprise  the  core  contributions  of  this 
research. 

Interference  temperature  is  a  non-issue  in  this  research  as  it  is  assumed  each  radio  can 
innately  and  accurately  identify  spectrum  whitespace  and  extract  signal  contents.  Likewise, 


Table  2.1:  Recent  advances  in  CR  fundamentals. 


Area 

Advances 

Contribution 

CR  characteristics 

New  communications  and  networking 

Primary 

paradigms,  cognitive  capability,  etc. 

CR  functions 

Spectrum  sensing/analysis,  manage- 

Secondary 

ment/handoff,  and  allocation/sharing 

Network  architecture  and  ap¬ 

Spectrum  broker  entities,  dynamic  net¬ 

Primary 

plications 

work  reconfiguration,  etc. 
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Table  2.2:  Recent  advances  in  spectrum  sensing  and  analysis. 


Area 

Advances 

Contribution 

Interference  temperature 

Tolerable  interference  levels 

N/A 

Spectrum  sensing 

Energy  detection,  feature  detection, 

match  filtering/coherent  detection,  etc. 

Secondary 

Cooperative  sensing 

User  selection,  decision  fusion,  and  ef¬ 
ficient  information  sharing,  and  dis¬ 
tributed  cooperative  sensing 

Primary 

hardware-based  sensing  and  energy  detection  are  not  direct  contributions,  but  determining 
how  to  use  said  whitespace  is  a  large  part  of  forming  the  adaptive  hopset.  Further,  doing 
so  within  the  CRN  and  in  a  distributed,  collaborative  fashion  is  key  to  the  architecture’s 
implementation. 

Adaptive  clustering  and  hopset  selection  are  the  central  thrusts  of  this  research, 
so  licensed  spectrum  sharing,  power  control,  game  theory  in  spectrum  sharing,  and 
cooperation  enforcement  receive  no  contribution.  Spectrum  handoff,  cognitive  relay,  and 
routing  are  all  targets  of  continued  development  but  are  not  immediately  addressed  by  this 
research.  Similarly,  this  work  makes  use  of  advances  in  the  CR  MAC  layer  and  lays  the 
groundwork  for  a  control  channel  and  related  management  schemes,  but  does  not  contribute 
to  either  area.  Security  in  this  network  is  a  noble  consideration  and  one  that  must  be 
investigated.  This  research  does  not,  however,  address  security  in  order  to  maintain  public 
distribution  standards. 

While  this  research  does  use  existing  standards  and  regulations  as  guides,  contribu¬ 
tions  to  standards  are  not  central  to  the  work.  Research  implementations,  however,  are 
highly  relevant.  In  addition  to  being  a  research  implementation  in  itself,  this  work  builds 
on  several  existing  implementations. 
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Table  2.3:  Recent  advances  in  dynamic  spectrum  allocation  and  sharing. 


Area 

Advances 

Contribution 

Licensed  spectrum  sharing 

Spectrum  underlay  and  overlay 

N/A 

Media  Access  Control 

(MAC)  in  CRNs 

MultiMAC  protocols  and  applicability 

of  the  MAC  layer  to  spectrum  access 

schemes 

Secondary 

Spectrum  handoff 

Suspending  transmission  when  a  PU 

reappears  and  resuming  operation  via 

contingency  planning 

Secondary 

Cognitive  relaying 

Spectrum-opportunistic  packet  for¬ 
warding  within  a  CRN 

Secondary 

Spectrum  sensing  and  access 

Partially  observable  Markov  decision 

process  (POMDP),  optimality  of  my¬ 
opic  policies,  and  inclusion  of  SU 

residual  energy  and  buffer  state  in 

POMDP 

Primary 

Power  control  in  a  CRN 

Impact  of  transmission  power  on  avail¬ 
able  spectrum,  dynamic  programming 

for  optimal  power  and  rate  control,  col¬ 
laborative  power  sensing,  etc. 

N/A 

Control  channel  manage¬ 
ment 

Cluster-specific  control  channel  based 

on  common  channels,  “swarm-based” 

selection  according  to  whitespace  anal¬ 
ysis,  etc. 

Secondary 

Distributed  spectrum  sharing 

Independent  action  for  fair  resource 

sharing,  time-spectrum  blocking 

Primary 
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(similar  to  GSM  cellular  technology), 

and  multi-agent  learning  for  resource 

management 

Spectrum  sharing  game 

Game  theory  in  spectrum  resource 

management,  no-regret  learning 

while  considering  cooperative/non¬ 
cooperative  users,  auction  mecha¬ 
nisms,  etc. 

N/A 

Routing  in  a  CR  network 

Accounting  for  lack  of  channel  com¬ 
monality,  multiple  channel  switches 

along  a  path,  variable  channel  switch¬ 
ing  delay,  etc. 

Secondary 

Cooperation  stimulation  and 

enforcement 

Assuming  not  all  nodes  are  uncondi¬ 
tionally  cooperative  in  system  design, 

cooperation  degree  metrics,  credit 

mechanisms,  etc. 

N/A 

Security  in  CRNs 

Localization-based  defense,  spectrum 

sensing  data  falsification  attacks,  in¬ 
duced  SU/PU  interference,  etc. 

N/A 
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Table  2.4:  Recent  advances  in  spectrum  sensing  and  analysis. 


Area 

Advances 

Contribution 

IEEE  802.22 

Relatively  new  standard  governing 

cognitive  radio  and  secondary  user  ac¬ 
cess  in  the  TV  spectrum  [13] 

Secondary 

IEEE  PI 900.1 

Standards  series  centered  around  spec¬ 
trum  management  and  next- generation 

radios  [12] 

Secondary 

Research  implementations 

Berkeley  Emulation  Engine  2  (BEE2) 

[14],  Public  Safety  Cognitive  Radio 

(PSCR)  at  Virginia  Tech  [15],  and 

Rugers  University’s  Open  Access  Re¬ 
search  Testbed  for  Next- Generation 

Wireless  Networks  (ORBIT),  etc. 

Primary 

This  research  takes  advantage  of  and  builds  on  numerous  recent  cognitive  radio  ad¬ 
vances  including  fundamental  characteristics,  dynamic  network  architectures,  cooperative 
sensing,  spectrum  access,  distributed  spectrum  sharing,  and  existing  research  implementa¬ 
tions.  These  advances  and  many  others  are  addressed  in  [3],  providing  a  broad  survey  from 
which  to  launch  into  the  actual  research  process.  Contributions  to  relevant  areas  are  further 
validated  in  this  chapter  and  the  Methdology. 

2.1.2  Governing  Standards. 

As  a  technology  growing  in  both  conceptual  popularity  and  practical  implementation, 
new  standards  are  needed  to  support  CR  use  [16].  National  governing  organizations 
including  the  American  Federal  Communications  Commission  (FCC)  and  the  British 
Office  of  Communications  (Ofcom)  are  already  considering  adopting  laws  to  govern 
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cognitive  radio  usage-specifically,  the  ability  of  SUs  to  use  fallow  RF  spectrum  allocated  to 
PUs.  Additionally,  multiple  organizations  including  the  International  Telecommunications 
Union  -  Radio  Sector  (ITU-R),  SDR  Forum,  and  the  IEEE  have  begun  work  on  CR-centric 
standards.  Namely,  the  IEEE  has  already  published  two  standards  as  they  relate  to  CR, 
SCC41  (formerly  P1900)  and  IEEE  802.22. 

The  authors  of  [16]  also  note  that  the  emergence  of  such  standards  is  particularly 
useful  when  multiple  users  are  competing  for  the  same  RF  spectrum  whitespace.  When 
two  users  compete  in  this  fashion,  they  often  are  forced  by  such  standards  to  share  channels 
with  applicable  users.  This  is  known  as  self-c oexistence  and  is  implicitly  applicable  to  this 
research  as  the  proposed  system  blindly  adapts  to  unavailable  spectrum.  In  this  way,  this 
work’s  CRN  performs  self-coexistence  naturally  as  new  users  make  their  presence  in  the 
spectrum. 

A  detailed  examination  of  the  IEEE  802.22  standard  is  presented  in  [13].  The  sticking 
point  of  this  article  is  that  IEEE  802.22  is  the  first  worldwide  standard  based  on  cognitive 
radio.  Although  the  article  is  from  2005,  it  serves  to  highlight  the  growing  prominence 
of  CR  as  the  solution  to  spectrum  overcrowding.  According  to  the  article,  in  2005  it  was 
estimated  that  at  any  given  time,  only  5.2%  of  the  RF  spectrum  is  usable  by  SUs-even 
though  PUs  often  leave  their  frequencies  dormant  and  unused.  The  authors  also  define  the 
standard’s  relation  to  existing  network  standards  (see  Figure  2.1),  where  the  outermost  ring 
is  the  proposed  CR  domain.  Although  not  designed  to  operate  at  higher  throughput  rates 
(i.e.,  18-24  Mbps  compared  to  IEEE  802.1  In  at  over  100  Mbps),  IEEE  802.22  provides  the 
key  functionality  of  allowing  SUs  to  opportunistically  operate  in  PU-designated  (but  often 
vacant)  spectrum. 
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Figure  2.1:  Relation  of  IEEE  802.22  to  other  IEEE  network  standards  [13]. 
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2.2  Adaptive  Frequency  Hopping  (AFH) 

2.2.1  Legacy  Frequency  Hopping  Systems. 

2.2. 1.1  Single  Channel  Ground  and  Airborne  Radio  System  (SINCGARS). 
SINCGARS  is  a  cryptographically-keyed  frequency  hopping  radio  originally  intended 
for  ground  forces.  The  radio  can  operate  over  2,320  channels  from  33  to  80  megahertz 
(MHz)  and  has  modular  construction  to  ease  the  upgrade  process.  The  operator  is  able 
to  adjust  the  radio’s  functionality  using  a  collection  of  controls  on  the  front  of  the  device, 
shown  in  Figure  2.2  [17,  18]. 


Figure  2.2:  Standalone  SINCGARS  radio  [19]. 


2.2. 1.2  HAVE  QUICK. 

HAVE  QUICK  is  the  U.S.  military’s  standard  ultra-high  frequency  (UHF)  radio 
standard.  Implemented  on  a  host  of  airborne  platforms,  HAVE  QUICK  is  a  sophisticated 
frequency  hopping  technology.  Similar  to  SINCGARS,  it  uses  a  cryptographic  key  to 
produce  a  frequency  hopping  sequence  that  is  chronologically  synchronized  among  all 
users  for  a  given  time  period.  Figure  2.3  shows  the  cockpit  panel  mount  in  a  NATO  F- 
16  aircraft  [18,  20] 
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Figure  2.3:  HAVE  QUICK  panel  mount  in  an  F-16  [20] 


2.2. 1.3  High  Frequency  (HF)  AFH. 

Adaptive  frequency  research  has  been  ongoing  since  the  early  1990s.  In  their 
1993  MILCOMM  paper,  the  authors  of  [21]  introduce  a  scheme  for  AFH  based  on 
observing  channel  states  and  generating  a  frequency  hopping  sequence  that  uses  “good” 
channels  more  often  than  “bad”  channels.  Also  mentioned  is  the  stipulation  that  available 
frequencies  are  not  selected  uniformly.  This  system  is  modeled  in  Figure  2.4.  Later,  in 
1996,  another  approach  for  AFH  involving  correlated  frequency  hopping.  While  not  overtly 
focused  on  AFH,  [22]  does  highlight  the  need  for  improved  frequency  hopping  measures 
for  the  purpose  of  avoiding  “limited  bandwidth  or  intentional  jamming”-both  of  which  are 
addressable  using  CR  technology. 

2.2. 1.4  Summary. 

SINCGARS  and  HAVE  QUICK  represent  modem  U.S.  and  NATO  frequency  hopping 
radio  standards.  However,  each  implementation  is  not  adaptable  to  the  current  RF  spectrum 
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Figure  2.4:  Early  AFH  system  [21]. 


and,  even  if  the  RF  environment  is  known,  requires  a  “human-in-the-loop”  to  manipulate 
the  radio  settings  to  fit  such  an  environment.  By  definition,  cognitive  radio  has  the  innate 
ability  to  sense  the  spectrum,  orient  its  operation,  decide  on  a  course  of  action,  and  perform 
the  appropriate  actions.  A  human  operator  must  engage  the  same  process,  but  because  a  CR 
can  autonomously  tweak  the  required  knobs  and  meters  to  fit  the  existing  RF  environment, 
it  can  operate  much  faster  and  at  the  expense  of  little  or  no  valuable  manpower.  AFH 
systems  began  to  emerge  in  the  early  1990s,  but  not  in  the  context  of  a  truly  cognitive 
radio.  As  compared  to  legacy  systems  such  as  SINCGARS  and  HAVE  QUICK,  CR  allows 
for  a  quicker,  more  precise  method  for  implementing  a  frequency  hopping  radio.  The  AFH 
iteration  presented  in  this  research  is  implemented  as  a  CR  with  all  necessary  automation. 

2.2.2  Bluetooth  (802.15.1 )  &  WLAN  (802.11 ). 

Adaptive  frequency  hopping  is  an  emerging  trend  in  CR  networks.  As  they  are  already 
able  to  sense  spectrum  availability,  CR  networks  are  natural  candidates  for  expanding  the 
adaptive  envelope  to  include  entire  hopsets  in  addition  to  single  frequencies.  An  example  of 
AFH  in  practice  is  the  current  Bluetooth  implementation.  Whereas  the  early  standard  was 
spread  out  across  79  of  the  possible  83.5  ISM  channels,  Bluetooth  AFH  versions  need  only 
hop  over  15  channels.  This  allows  for  much  less  interference  to  other  devices  operating  in 
the  same  band,  including  WEAN  (802.11)  and  ZigBee  (802.15.4)  [23]. 
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Adaptive  frequency  hopping  is  an  emerging  trend  in  CR  networks.  As  they  are  already 
able  to  sense  spectrum  availability,  CR  networks  are  natural  candidates  for  expanding  the 
adaptive  envelope  to  include  entire  hopsets  in  addition  to  single  frequencies.  An  example  of 
AFH  in  practice  is  the  current  Bluetooth  implementation.  Whereas  the  early  standard  was 
spread  out  across  79  of  the  possible  83.5  ISM  channels,  Bluetooth  AFH  versions  need  only 
hop  over  15  channels.  This  allows  for  much  less  interference  to  other  devices  operating  in 
the  same  band,  including  WLAN  (802.11)  and  ZigBee  (802.15.4)  [23,  24]. 

2.2.3  Dynamic  Adaptive  Frequency  Hopping. 

The  authors  of  [25]  propose  the  dynamic  adaptive  frequency  hopping  (DAFH) 
algorithm  in  wireless  personal  area  networks  (WPANs).  They  define  DAFH  as  a  general 
method  that  enables  the  FH-networks  in  the  unlicensed  band  to  avoid  mutual  interference 
in  a  distributed  manner.  Their  approach  is  similar  to  ours  in  that  they  determine  available 
frequencies  and  select  a  hopset  using  a  common  seed  for  pseudorandom  hopset  generation. 
The  authors  also  implement  per-subnet  band  division,  i.e.  the  overall  band  is  partitioned 
such  that  every  neighboring  subnet  uses  a  separate  frequency  sample  space  to  avoid 
collisions  between  hopsets.  We  currently  intend  to  simply  ensure  no  two  subnets  will  share 
the  same  hop  at  the  same  time,  but  the  technique  in  [25]  is  certainly  worthwhile  to  explore 
for  future  development. 

2.2.4  Summary. 

Adaptive  frequency  hopping  (AFH)  forms  the  core  basis  by  which  the  proposed 
system  reliably  and  unpredictably  avoids  interference.  Frequency  hopping  applications 
currently  drive  military  radio  communication  systems,  but  none  hops  according  to 
adaptation  to  the  RF  environment.  Commercial  applications  such  as  Bluetooth  have  begun 
implementing  simple  packet  loss-induced  AFH  schemes,  but  such  methods  continue  to 
use  occupied  frequencies  to  avoid  changing  the  predefined  hopset.  Older  AFH  schemes 
exist  primarily  in  the  HF  bands,  but  newer  AFH  systems  operate  in  the  widely-used  ISM 
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and  UNII  bands  and  are  focused  on  hopset  collision  avoidance.  This  research  also  avoids 
hopset  collisions.  Additionally,  the  proposed  system  uses  a  combination  of  random  keying 
and  spectrum  sensing  to  accurately  avoid  occupied  frequencies  using  frequency  hopping. 

2.3  RF  Spectrum  Sensing  &  Mapping 

2.3.1  Cooperative  Sensing. 

Cooperative  sensing  consists  of  multiple  CRs  within  a  CRN  working  together  to 
formulate  an  accurate  picture  of  the  RF  environment.  This  topic  is  of  great  interest  in 
distributed  CR  systems  research  as  cooperative  sensing  both  minimizes  overall  workload 
for  the  individual  node  and  yields  a  more  reliable  depiction  of  environmental  conditions. 
However,  there  are  tradeoffs  to  cooperative  sensing:  the  CRN  must  both  communicate 
valuable  RF  spectrum  information  within  the  network  and  minimize  the  amount  of  network 
flow  over  low-bandwidth  and  potentially  vulnerable  control  channels. 

In  [26],  the  authors  examine  the  physical  aspects  of  the  motivation  for  cooperative 
sensing.  The  authors  recognize  bandwidth  as  a  confounding  factor  of  the  amount  of 
information  detectable  by  CRs:  if  a  radio  has  a  low-bandwidth  control  channel,  it  is 
likely  only  able  to  perform  either  energy  detection  or  signal  statistics;  a  CR  with  a  high- 
bandwidth  control  channel,  however,  may  implement  all  possible  detectors.  The  type  of 
control  channel  is  largely  dependent  on  the  hardware  used  to  implement  the  board. 

In  line  with  opinions  presented  in  [26],  the  authors  of  [27]  also  examine  the  issue 
of  bandwidth  limitations  in  cooperative  sensing.  Nodes  must  report  within  the  network 
their  observations,  and  the  granularity  of  these  observations  is  limited  by  the  bandwidth 
delivered  in  hardware.  The  implementation  studied  in  [27]  uses  sensing  bits  to  quantize  a 
likelihood  ratio  (LR)  for  determining  which  nodes  have  enough  worthwhile  information  to 
disseminate  through  the  network.  Using  LRs  is  thereby  a  means  for  reducing  meaningless 
observation  traffic  within  the  network  via  low-bandwidth  control  channels.  Their  tests  show 
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that  with  only  a  minor  cost  to  sensing  performance,  the  number  of  sensing  bits  decreases 
greatly. 

Cooperative  sensing  is  also  very  applicable  to  cluster-based  networks,  as  noted  in  [28]. 
Sensing  in  this  manner  minimizes  overall  traffic  throughout  the  overall  CRN  and  allows  for 
sensing  to  be  geographically  localized.  This  is  important  because  geographically  distant 
nodes  are  unlikely  to  sense  and  attempt  to  use  the  same  RF  spectrum  environment  as  nodes 
within  close  proximity  of  one  another. 

Cooperative  sensing  plays  a  large  part  in  this  network  in  terms  of  data  fusion  (RF 
map  merging),  and  in  order  to  minimize  overhead  time  associated  with  adapting  to  the  RF 
environment,  limiting  extraneous  information  over  control  channels  is  vital.  This  research 
proposes  minimizing  network  traffic  not  via  minimal  data  representation  as  in  [27],  but 
instead  by  partitioning  the  network  in  to  sub-networks  (“subnets”)  using  clustering.  There 
are  two  distinct  benefits  to  partitioning  the  network  into  clusters  that  are  likely  to  operate 
in  the  same  RF  environment. 

2.3.2  Spectrum  Map  Storage  &  Usage. 

An  RF  spectrum  can  be  represented  in  many  different  ways,  such  as  real+complex 
number  pairs,  a  histogram-like  FFT  model,  and  so  forth.  Prior  to  conducting  this 
experiment,  it  is  necessary  to  choose  a  method  such  that  it  is  both  portable  between 
components  and  extensible  for  future  research.  One  method  generalizes  the  RF  spectrum  as 
part  of  a  larger  set  of  data,  the  radio  environment.  In  their  pioneering  Radio  Environment 
Map  (REM)  research,  Zhao,  Le,  and  Reed  present  the  RF  spectrum  as  part  of  a  larger 
model  containing  multiple  quantitative  and  qualitative  radio  environment  characteristics  to 
include  terrain,  historical  or  time-based  usage  data,  and  others  in  addition  to  the  obvious 
inclusion  of  present  RF  spectrum  usage  [10]. 

This  research  implements  the  REM  solely  as  an  avenue  for  storing  spectrum 
availability  data  where  frequency-specific  data  is  contained  in  a  binary  vector.  This  concept 
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is  validated  in  section  three.  Therefore,  and  for  the  rest  of  this  paper,  the  spectrum  map  is 
referred  to  simply  as  the  REM  to  allow  for  expansion  of  this  research  to  include  more 
aspects  of  the  radio  environment. 

2.3.2. 1  Radio  Environment  Map. 

Zhao,  Le,  and  Reed  introduce  the  radio  environment  map  (REM)  as  an  integrated 
database  consisting  of  multi-domain  information  such  as  geographical  features,  available 
services,  spectral  regulations,  locations  and  activities  of  radios,  policies  of  the  user  and/or 
service  provider,  and  past  experience.”  Figure  2.5  illustrates  this  concept. 

A  cognitive  radio  (CR)  uses  this  REM  to  improve  its  situational  awareness  (SA)  and 
reasoning  process,  and  it  disseminates/receives  REM  information  to/from  other  CRs  in 
order  to  form  a  highly  situation-aware  CR  network.  According  to  the  authors,  the  formation 
of  a  REM  is  both  a  system-level  solution  to  cognitive  networking”  and  a  natural  step 
from  legacy  radios  to  CR;  more  succinctly,  REM  usage  in  CR  networks  leverages  prior 
knowledge  and  collective  intelligence.”  Finally,  the  authors  of  [10]  recognize  five  open 
questions  related  to  the  REM,  two  of  which  apply  specifically  to  this  thesis: 

1.  How  can  we  assure  the  REM’s  integrity,  security,  privacy,  and  reliability? 

2.  In  order  to  provide  desired  performance,  how  current  and  what  level  of  granularity 
does  the  information  contained  in  the  REMs  need  to  be? 

Regarding  these  questions,  this  thesis  concerns  the  REM  dissemination  integrity,  REM 
reliability,  and  REM  currency. 

2. 3. 2. 2  Overhead  Analysis  for  REM-enabled  Cognitive  Radio  Networks. 

An  important  aspect  of  using  a  REM  in  conjunction  with  a  CR  is  the  overhead 
associated  with  such  action.  Zhao,  Reed,  Mao,  and  Bae  [29]  present  REM  dissemination 
schemes  and  issues,  analyze  REM  dissemination  overhead,  and  compare  various  overhead 
scenarios  with  network  simulation  results.  This  research  mainly  pertains  to  REM 
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Figure  2.5:  Possible  REM  characteristics  [29]. 


dissemination  schemes  and  issues  with  an  eventual  focus  on  simulating  the  network  and 
analyzing  subsequent  REM  overhead. 

Dissemination  schemes  include  periodic  REM  broadcasts  to  the  entire  network  (via 
plain  flooding),  an  CR-specific  extension  of  the  mobile  wireless  LAN  optimized  link 
state  routing  (OLSR)  protocol  (via  adapting  or  extending  the  OLSR  HELLO  or  Topology 
Control  (TC)  messages),  adapting  the  REM  dissemination  rate  based  on  existing  primary 
users  (PUs)  or  interference  in  an  application-specific  ad  hoc  network,  and  disseminating 
REM  information  in  an  on-demand”  scheme  (i.e.  only  transmitting  requested  REM 
information,  as  opposed  to  always  transmitting  the  entire  map)  [29]. 

When  choosing  the  REM  dissemination  scheme,  the  overarching  goals  are  minimal 
retransmissions,  minimal  transmission  sizes,  and  a  minimal  number  of  REM  sources.  REM 
implementations  can  be  any  data  structure  containing  the  desired  information.  In  ascending 
computational  overhead,  the  implementation  can  take  the  form  of  a  multi-dimensional 
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matrix,  a  C++  structure,  or  a  multi-dimensional  database.  REM  memory  footprint  depends 
on  both  the  implementation  and  the  amount  of  information  stored  in  the  REM.  The  authors 
suggest  using  a  common  control  channel  in  any  of  the  following  flavors:  Narrowband 
channel  in  a  licensed  band;  channel  in  license-free  ISM  or  UNII  band;  an  ultra-wideband 
(UWB)  channel;  or  sharing  with  the  traffic  channel  [29]. 

While  the  aforementioned  methods  are  valid,  this  research  proposes  REM  dissemina¬ 
tion  by  a  different  means.  Totally-ordered  multicast  (TOMC)  is  an  information  dissemina¬ 
tion  protocol  in  which  all  nodes  receive  messages  in  the  same  order  [30].  Such  a  protocol 
is  valuable  in  a  network  where  information  must  be  distributed  among  a  large  number  of 
nodes  and  network  structure  decisions  are  contingent  upon  that  information.  This  is  essen¬ 
tially  the  state  of  the  proposed  network  during  formation  and  prior  to  clustering.  Nodes 
must  receive  information  in  the  same  order  to  guarantee  that  each  node  receives  the  same 
data  in  the  same  order.  While  TOMC  is  not  explored  further  than  conceptual  integration  in 
this  document,  it  is  expected  to  become  a  vital  part  of  future  research  involving  inter-node 
communication. 

Table  2.5  describes  a  potential  implementation  for  REM  information,  ideally 
contained  in  a  C++  structure  or  multi-dimensional  database.  The  rightmost  column  shows 
how  each  field  applies  to  the  system  proposed  in  this  research,  if  applicable.  While  the 
IEEE  802.22  standard  is  designed  to  be  infrastructure-based,  our  CRN  is  is  designed  as  an 
infrastructure-less,  ad  hoc  network.  This  table  helps  define  where  the  proposed  system  is 
similar  to  an  infrastructure-based  CRN  like  one  governed  by  IEEE  802.22  [10]. 
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Table  2.5:  REM  information  element  for  IEEE  802.22  systems  [10]. 


Domain  and 

Application 

index  range 

Syntax  and  index 

to  research 

Application  type 

Voice  (701),  packet  data  (702),  vido  conference 

Packet  data 

(700-799) 

(703),  etc. 

Optimization 

Minimize  interference  to  PU  (600),  maximize 

Maximize  SU 

layer (600-699) 

SU  throughput  (601),  etc. 

goodput 

Topology  and  net¬ 

Infrastructure -based  network  {WCDMA  (500), 

ad  hoc  network 

work  type  (500- 

cdma2000  (501),  WRAN  (502),  etc.};  Ad  hoc 

599) 

network  (510),  mesh  network  (520),  etc. 

MAC  and  duplex 

TDMA  (400),  FDMA  (401),  CDMA  (402), 

Frequency  hop¬ 

(400-499) 

OFDMA  (403);  FDD  (410),  TDD  (411),  etc. 

ping 

Geography  and 

Indoor  {home  (300),  office  (301),  airport  (302), 

Various 

mobility  informa¬ 

factory  (303),  etc};  outdoor  {urban  (310), 

tion  (300-399) 

suburban  (311),  open  rural  (312),  highway 

(313),  etc.};  in-vehicle  {train  (320),  bus  (321), 

car(322),  plane  (323),  etc.};  etc. 

Modulation  type 

AM  (200);  FM  (210);  M-PSK  {BPSK  (220), 

N/A 

(200-299) 

QPSK  (221),  etc};  M-QAM  {16-QAM  (230), 

64-QAM  (231),  etc.};  etc. 

Radio  device 

Channel  coding  {Convolutional  coding  (100)}, 

Sensitivity,  opera¬ 

capability  (1  GO- 

Turbo  coding  (110),  etc.};  maximum  RF  trans¬ 

tional  frequencies 

199) 

mit  power  (120),  sensitivity  (130),  operational 

bands  (140);  antenna  type  (150),  etc. 

Experience  (0-99) 

Blind  zone  (10),  hot  spot  (20),  hidden  node 

Dynamic  RF  en¬ 

(30),  etc. 

vironment 
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2.3.23  Radio  Environment  Map-enabled  Learning  Algorithms. 

The  authors  of  [31]  propose  a  generic  top-down  approach  for  obtaining  situational 
awareness  (SA)  via  REM  exploitation  for  CRs  and  a  framework  for  CR  learning  loops. 
Additionally,  the  authors  conduct  simulations  to  evaluate  the  efficiency  and  effectiveness 
of  the  approach  and  framework  they  propose.  Metrics  include  adaptation  time,  average 
received  signal  to  interference  and  noise  ratio  (SINR),  average  throughput,  and  existing  PU 
average  packet  delay. 

The  authors  use  Table  2.5  to  describe  one  possible  method  for  digitizing  and  indexing 
radio  environment  information.  For  example,  suppose  that  the  REM  shows  that  the  CR  is 
within  the  service  area  of  a  wireless  regional  area  network  (WRAN),  and  that  the  REM 
also  lists  TV  channel  9  as  available  for  use.  Once  the  CR  senses  the  open  channel  and 
adjacent  channels  and  confirms  that  TV  channel  9  is  indeed  available  and  suitable  for  use, 
the  WRAN  determines  the  current  situation,  which  in  this  example  can  be  defined  by  the 
Table  2.5  as  {703,  600,  502,  403,  411,  312,  231,  etc.}.  Thus  different  REM/CR  scenarios 
can  be  modeled  across  different  CRs  [29,  31]. 

The  use  of  the  REM  has  also  been  explored  as  a  means  to  spatially  interpolate  radio 
locations.  The  authors  of  [32]  present  a  new  robust  method  for  using  inverse  distance 
weighting  to  determine  feature  locations.  The  REM  forms  the  backbone  of  several  other 
areas  of  research,  including  building  a  mobile  ad  hoc  network  architecture  around  the  REM 
[33]  and  investigating  the  effect  of  a  noise  threshold  on  REM  effectiveness  [34].  Both 
concepts  are  central  to  this  research-the  REM  is  the  primary  means  for  storing,  using,  and 
transferring  this  system’s  RF  environment  observations,  and  the  applicability  of  a  noise 
threshold  is  also  thoroughly  vetted  when  comparing  the  RF  spectrum  environments  of 
neighboring  CRs. 
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2.3.3  Threshold  Detection. 


Threshold  detection  is  an  important  part  of  extracting  valuable  information  from  a 
signal.  While  developing  a  threshold  detection  method  is  not  a  focus  of  this  research, 
implementing  an  existing  solution  plays  a  major  role  in  both  the  clustering  evaluation  and 
hopset  selection  steps  of  the  overall  cycle  of  operations.  The  authors  of  [35]  highlight  one 
such  method  using  the  noise  floor.  The  noise  floor  is  determined  to  be  the  mean  of  the  entire 
spectral  density  function,  and  the  threshold  is  the  noise  floor  (mean)  multiplied  by  some 
coefficient.  This  coefficient  serves  to  move  the  threshold  above  the  noise  floor  by  some 
amount.  In  this  research,  the  threshold  is  implemented  in  a  similar  manner  for  the  purpose 
of  converting  the  measured  spectrum  to  a  format  usable  for  comparison  and  lightweight 
hardware  manipulation.  The  conversion  method  and  the  usage  of  the  conversion  output  is 
discussed  in  Section  3.2.5. 1  and  validated  both  Sections  3. 2. 5. 3  and  3. 3. 4. 2. 

2.3.4  Summary. 

In  order  for  a  radio  to  accurately  sense  its  environment  and  correctly  adapt  its 
behavior,  it  must  be  able  to  “see”  its  surroundings  and  “remember”  how  they  affect 
its  operation.  Many  different  techniques  for  sensing  in  the  context  of  CRNs  currently 
exist  and  have  been  demonstrated  as  effective  purveyors  of  RF  environment  information. 
The  radio  environment  map  (REM)  provides  an  ideal  medium  for  efficiently  storing  and 
communicating  RF  environment  information.  While  this  research  does  not  implement  a 
spectrum  sensing  implementation,  a  method  for  mapping  the  spectrum  and  storing  it  in  a 
custom  REM  is  proposed  and  demonstrated. 

2.4  Clustering 

2.4.1  Overview. 

Clustering  is  the  manner  in  which  similar  objects  are  “clustered”  together  using  some 
distance  or  similarity  metric.  In  practice,  clustering  serves  a  number  of  purposes.  In 
statistical  analysis,  clustering  is  a  way  to  identify  groupings  of  data.  In  image  processing, 
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clustering  enables  similarly-colored  pixels  to  be  marked  as  such.  This  research  uses 
clustering  to  form  IP  subnets  from  an  initially-flat  radio  network. 

2.4.2  k-means  Clustering. 

As  Kanungo,  et  al  explain  in  [36],  although  there  are  a  variety  of  k-means  clustering 
implementations,  Lloyd’s  algorithm  is  known  effectively  as  the  defacto  k-means  algorithm. 
Before  defining  the  algorithm,  though,  it  is  necessary  to  define  terms.  First,  a  centroid  is 
simply  the  geometric  weighted  center  of  a  cluster.  For  example,  if  all  nodes  in  a  cluster 
are  of  equal  weight,  the  centroid  is  the  center  of  mass  of  the  shape  composed  by  the  nodes. 
Second,  a  center  is  the  point  location  of  the  centroid.  After  the  clustering  algorithm  runs, 
every  node  is  assigned  to  a  center.  Finally,  a  center’s  neighborhood  is  defined  as  those 
nodes  for  which  z  is  the  closest  neighbor.  The  focal  point  of  Lloyd’s  algorithm  is  observing 
the  optimal  placement  of  the  center  of  at  the  centroid  of  the  relevant  cluster.  The  actual 

algorithm  is  defined  below. 

1:  procedure  LLOYD(Some  set  of  k  centers  Z) 

2:  repeat 

3:  for  each  center  zeZdo 

4:  Let  V(z)  denote  z’s  neighborhood 

5:  Move  z  to  the  centroid  of  V(z) 

6:  Update  V(z)  with  the  distance  from  each  point  to  its  nearest  center 

7:  end  for 

8:  until  some  convergence  condition  is  met 

9:  end  procedure 

2.4.3  University  of  Maryland  Testbed. 

The  authors  of  [37]  present  a  framework  for  testing  the  k-means  clustering  algorithm. 
This  testbed  is  written  in  C++  and  forms  the  centerpiece  of  this  work’s  clustering 
experiment.  A  total  of  34  different  runtime  options  are  available.  For  simplicity,  and 
to  limit  the  number  of  experiments,  this  research  varies  several  and  maintains  the  rest 
as  default  values.  One  of  the  most  important  options  is  the  ability  to  select  among  four 
different  clustering  algorithms,  all  of  which  are  variations  on  the  core  Lloyd’s  algorithm. 
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The  testbed  uses  the  original  Lloyd’s  algorithm  and  three  more;  the  additional  three 
algorithms  involve  some  modification  or  addition(s)  to  the  Lloyd’s  algorithm.  All  four 
algorithms  are  listed  and  summarized  below.  Full  explanations  are  contained  in  [37]. 

•  Lloyd’s.  Original  centroid-based  algorithm  that  runs  until  convergence  of  cluster 
assignments.  Initial  centers  are  randomly  sampled. 

•  Swap.  Maintains  a  set  of  candidate  centers  and  swaps  between  this  list  and  the 
existing  centers. 

•  Hybrid.  Complex  hybrid  algorithm  comprised  of  both  Lloyd’s  and  Swap  in  which  a 
number  of  Swap  iterations  are  performed  followed  by  several  iterations  of  Lloyd’s. 
The  Hybrid  algorithm  uses  an  approach  akin  to  simulated  annealing  to  avoid  local 
minima  entrapment. 

•  EZ-Hybrid.  Simplified  Hybrid  algorithm.  Performs  one  swap  and  a  number  of 
subsequent  iterations  of  Lloyd’s. 

2.4.4  Cluster  Visualization. 

Every  node  is  associated  with  a  particular  RF  spectrum  measurement.  These 
measurements  can  be  compared  between  two  nodes  to  yield  a  similarity  ratio  ranging  from 
totally  dissimilar  (zero)  to  totally  similar  (one).  While  several  cluster  visualizations  are 
informally  presented  in  [36],  this  research  shows  both  cluster  assignments  via  projected 
communication  links  and  the  RF  spectrum  similarity  between  nodes. 

2.4.5  Summary. 

A- means  clustering  is  the  basis  of  the  middleware’s  ability  to  reduce  a  flat  network 
into  a  hierarchical  set  of  sub-networks.  In  the  context  of  this  research,  an  IP-based  network 
composed  of  cognitive  radios  is  partitioned  into  subnets.  While  the  A- means  algorithm 
itself  is  not  new,  implementing  the  algorithm  in  the  realm  of  congitive  radio  networks  and 
the  manner  in  which  experimental  data  is  visualized  is  a  novel  contribution. 
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2.5  FPGA-Based  Cognitive  Radio 

2.5.1  Kansas  University  Agile  Radio  (KUAR). 

The  University  of  Kansas  uses  a  Virtex-II  Pro  FPGA  for  the  digital  signal  processing 
operations  and  digital  communication  components  of  their  software-defined  radio  platform 
KUAR.  It  performs  communications  processing  in  VHDL  on  the  FPGA,  does  signal 
processing,  radio  control,  and  RF  environment  sensing  with  both  VHDL  on  the  FPGA 
and  C  code  on  the  FPGAs  embedded  processor,  and  interfaces  with  a  transceiver  and  a 
Linux  control  processor  [38]. 

2.5.2  Wireless  Open-Access  Research  Platform  (WARP). 

Rice  University  uses  the  PowerPC  processor  on  the  Virtex-II  Pro  FPGA  for 
communication  processing  on  their  Wireless  Open-Access  Research  Platform.  WARP 
is  used  to  prototype  wireless  networks.  They  have  tested  and  verified  the  FPGAs  on 
the  WARP  board  in  over-the-air  tests.  The  single-input  single-output  (SISO)  Orthogonal 
Frequency  Division  Multiplexing  (OFDM)  transceiver  uses  the  FPGA  for  all  the  baseband 
processing.  It  uses  the  radio  daughtercards  on  the  WARP  board  to  convert  to  the  RF  band 
[39], 

2.5.3  Trinity  College’s  Cognitive  Radio  Framework. 

The  authors  of  [40]  at  Trinity  College  performed  a  case  study  that  demonstrated  the 
implementation  of  cognitive  radio  on  an  FPGA.  They  created  a  set  of  tools  for  radio 
designers  to  be  able  to  implement  cognitive  radio  on  FPGAs.  This  research  used  a  partial 
reconfiguration  of  the  FPGAs  during  runtime,  allowing  a  cognitive  engine  to  reconfigure 
both  the  software  and  hardware. 

2.5.4  Berkeley  Emulation  Engine  2  (BEE2). 

The  BEE2  is  a  modular  radio  prototyping  testbed  for  examining  both  narrowband 
and  wideband  approaches,  and  uses  a  combination  of  Simulink  and  Linux  alongside  the 
BORPH  operating  system  in  implementation.  In  their  experimental  design,  the  authors  of 
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[14]  use  the  BEE2  as  the  basis  for  two  separate  examples,  a  wideband  configurable  testbed 
in  the  400-500  MHz  range,  and  a  narrowband  multiple-input/multiple-output  CR  example 
in  the  20  MHz  range.  This  example  uses  a  Virtex  II  Pro  FPGA  as  the  central  part  of  the 
BEE2’s  CR  extension,  much  like  the  implementation  proposed  in  this  research. 

2.5.5  Virginia  Tech  Public  Safety  Cognitive  Radio  (PSCR). 

The  authors  of  [15]  present  a  CR  implementation  specifically  for  public  safety  (police, 
fire,  etc.)  personnel,  Public  Safety  Cognitive  Radio  (PSCR).  The  research  implements 
the  traditional  operation  cycle  of  a  cognitive  radio  in  very  fine  detail.  Included  in 
the  implementation  is  the  radio  environment  map  (REM)  much  like  in  this  document’s 
proposed  system.  In  addition  to  the  radio’s  “cognitive  cycle”  (as  also  defined  with  this 
research),  the  authors  present  a  fully-described  product  from  physical-layer  hardware  up 
through  the  graphical  user  interface  (GUI).  One  important  feature  of  the  PSCR  is  its 
dependency  on  policy-based  verification  such  that  the  radio  fits  both  the  physical  RF 
spectrum  environment  and  the  public  policies  governing  spectrum  usage  [15]. 

2.5.6  Summary. 

Like  those  in  previous  work,  this  research  leverages  the  FPGA  features  of  flexibility 
and  compact  area.  In  using  the  base  WARP  system,  a  custom  IP  core  written  in  VHDL,  and 
the  Virtex  IV’s  embedded  PowerPC  core  with  C  code  to  implement  the  radio  and  hopset 
selector. 

2.6  Background  Summary 

Many  underlying  technologies  compose  what  we  envision  to  be  a  viable,  field- 
deployable  cognitive  radio  network  and  the  supporting  architecture.  While  some  existing 
systems  are  indeed  similar  to  the  one  proposed  in  this  research,  this  work  employs  a  wholly 
new  architecture  and  implements  such  components  in  a  new  and  novel  manner.  Dynamic 
adaptive  frequency  hopping,  network  clustering,  and  packaging  the  system  into  a  hybrid 
hardware/software  FPGA  are  the  core  pieces  of  this  work  and  have  already  been  explored 
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or  accomplished  individually.  This  research  shows  that  such  sub-systems  can  be  integrated 
into  a  novel  cognitive  radio  middleware  architecture  as  part  of  an  ongoing  effort  to  develop 
a  “smarter  radio.” 
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3  Methodology 


This  research  proposes  the  architecture  for  a  frequency  hopping  cognitive  radio 
network  which  coexists  with  other  RF  spectrum  users,  and  implements  several  of 
the  necessary  features.  Specifically,  the  goals  of  this  research  are  to  answer  the  following 
questions: 

1.  Is  it  possible  and  feasible  to  implement  such  a  system? 

2.  What  are  the  properties  of  such  a  system? 

3.  What  are  the  pieces  for  implementing  such  a  system? 

4.  How  is  the  system  implemented? 

5.  What  is  the  operation  of  the  system? 

Similarly,  the  following  hypotheses  are  investigated: 

1.  Feasibility  of  clustering  an  ad  hoc  CR  network  based  on  spectrum  similarity  and 
physical  distance.  The  network  is  partitioned,  or  clustered,  according  to  the 
similarities  between  spectra  at  individual  nodes.  It  is  expected  clustering  a  CR 
network  according  to  local  RF  spectrum  similarity  and  physical  distance  produces 
clusters  with  nodes  whose  sensed  spectra  are  more  similar  than  clustering  on  physical 
distance  alone. 

2.  Feasibility  ofAFH  in  an  interference  environment.  Since  the  core  focus  of  this  system 
is  to  adaptively  select  and  use  a  frequency  hopping  sequence,  or  hopset,  based  on 
the  existing  RF  spectrum,  it  entails  that  this  research  must  prove  such  operation  as 
tractable.  It  is  expected  the  system  can  recognize  available  and  unavailable  spectrum, 
react  with  an  appropriate  AFH  sequence,  and  continue  normal  operation. 
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3.  Feasibility  of  communicating  information  between  nodes  in  an  orderly,  organized 
fashion.  Nodes  must  have  the  ability  to  form  and  maintain  a  dynamic  network,  share 
spectrum  data,  and  achieve  synchronized  frequency  hopping  activity.  Since  multiple 
nodes  must  communicate  on  the  same  hopping  pattern,  they  must  chronologically 
share  spectrum  resources.  Further,  the  spectrum  may  change  independently 
at  different  network  locations.  To  maintain  adaptive  operations,  nodes  must 
communicate  spectrum  data,  engage  their  cognitive  operations,  and  (re)distribute  an 
adaptive  hopset.  It  is  expected  all  nodes  can  maintain  effective  communications  while 
conducting  AFH  operations. 

4.  Ability  to  successfully  transmit  in  the  presence  of  a  dynamic  RF  spectrum.  In  addition 
to  deciding  an  adaptive  hopset,  the  system  must  transmit  and  receive  between  radios. 
This  is  necessary  to  measure  typical  traffic-based  network  metrics  such  as  goodput 
and  latency.  It  is  expected  the  system  can  combine  the  first  three  hypotheses  to  form 
a  novel  AFH  CR  network  system. 

3.1  Whole  System 

This  section  presents  the  entire  adaptive  frequency  hopping  (AFH)  architecture.  When 
in  use,  the  system  is  able  to  coexist  with  other  spectrum  users.  This  section  explains  how 
the  system  functions  and  its  internal  structure  and  lists  which  components  are  implemented 
in  this  research. 

3.1.1  Assumptions. 

In  order  to  properly  scope  this  research  and  therefore  limit  the  design  decisions  for 
unimplemented  components,  the  following  assumptions  are  made: 

1.  Rendezvous  and  routing  are  complete  prior  to  middleware  actions. 

2.  Every  radio  is  within  transmission  range  of  all  other  radios. 
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3.  Every  radio  acts  as  a  node  in  an  IP  network  such  that  each  radio  has  both  an  IP  and 
MAC  address. 

4.  Radios  have  the  ability  to  communicate  with  one  another  prior  to  network  clustering. 

5.  The  network  is  immobile  and  the  spectrum  is  fixed,  i.e.,  nodes  cannot  enter  or  leave 
the  network  and  the  RF  spectrum  remains  constant  once  the  experiment  begins. 

6.  Nodes  are  uniquely  identified  based  on  network  arrival  time,  i.e.,  the  first  radio  to 
join  the  network  has  the  first  identifier. 

7.  The  radio  with  the  first  identifier  (e.g.  “Node  0”)  serves  as  the  network  leader. 

8.  Each  radio  can  detect  “quiet”  periods  and  engage  in  spectrum  sensing  during  these 
periods. 

9.  Nodes  have  the  ability  to  synchronize  their  time  upon  rendezvous  and  remain 
synchronized  for  the  duration  of  operations. 

10.  The  observable  spectrum  is  contiguous  and  spans  the  same  frequency  range  across 
the  network  such  that  all  REMs  are  of  equivalent  structure  during  map  merging. 

11.  Operations  occur  precisely  in  the  order  defined  in  Figure  1.4  such  that  transmission 
and  sensing  cannot  occur  simultaneously. 

12.  A  reliable  transfer  protocol  is  used  during  data  transmission  in  0-3  and  0-6  in  Figure 
1.4. 

3.1.2  Whole  System  Function. 

System  functionality  describes  the  “what”  of  the  proposed  architecture.  In  other 
words,  the  system  function  is  an  direct  answer  to  the  requirements  posed  at  the  beginning 
of  this  chapter  and  accounts  for  the  assumptions  given  above.  The  following  list  details 
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each  operation  in  chronological  order,  and  Figure  3.1  shows  how  these  operations  relate  to 
each  other  and  the  conditions  that  must  be  met  to  transition  between  operations. 

0-1.  Sense  RF  spectrum.  All  nodes  scan  the  radio  frequency  environment  to 
obtain  a  snapshot  of  occupied  frequency  bands  and  available  whitespaces.  Effectively, 
this  snapshot  amounts  to  a  2,048-bin  Fast  Fourier  Transform  (FFT)  of  the  analog  spectrum. 
A  threshold  is  applied  to  the  FFT  output  to  convert  the  spectrum  to  a  REM  representation. 


Figure  3.1:  Whole  system  function. 
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Each  node  applies  this  threshold  individually  and  stores  this  information  into  its  local  REM. 
This  threshold  is  discussed  in  Section  3.2.5. 1. 

0-2.  Form  clusters.  During  ad-hoc  network  formation,  each  node  must  perform 
rendezvous/neighbor  discovery.  Per  assumption  one  above,  this  is  already  accomplished 
prior  to  0-1.  For  the  purpose  of  this  research,  it  is  assumed  nodes  begin  in  a  non-hopping 
mode  of  operation,  where  discovery  occurs  on  a  fixed  control  channel.  Per  assumptions  two 
and  three  above,  the  radios  initially  form  one  large,  flat-hierarchy  IP  network.  Nodes  then 
form  geographic  clusters,  meaning  nodes  within  close  geographic  proximity  are  assigned  a 
common  cluster  identifier.  These  clusters  serve  as  IP  subnets,  and  each  cluster/subnet  has 
a  leader  given  by  assumption  six.  Such  subnets  serve  to  reduce  the  total  amount  of  traffic 
within  the  overarching  network.  It  is  hypothesized  nodes  within  a  cluster  have  similar 
REMs  and  are  able  to  form  a  common  REM  with  limited  data  exchange. 

0-3.  Share  REMs.  A  radio  may  freely  transmit  their  REM  information  at  any 
time,  but  will  only  do  so  if  its  observed  RF  spectrum  environment  has  changed.  This  is 
accomplished  via  TOMC.  In  TOMC,  many  nodes  can  attempt  to  “speak”  at  a  given  time, 
but  all  nodes  will  only  “hear”  from  one  node  at  a  time  due  to  the  nature  of  the  protocol.  In 
this  protocol,  all  intended  transmissions  are  eventually  received  in  a  totally-ordered  fashion, 
meaning  each  node  will  receive  all  data  in  the  same  order  as  all  other  nodes.  This  operation 
occurs  after  0-2  because  doing  so  in  a  flat,  complex  network  floods  the  network  with 
REM  data  traffic,  thereby  voiding  the  concept  of  using  clustering  to  minimize  unnecessary 
transmissions. 

0-4.  Merge  spectrum  maps.  The  leader  of  a  given  cluster  merges  the  maps  of  the 
cluster’s  nodes  to  form  a  common  REM.  Because  REMs  are  represented  as  binary  data, 
the  common  REM  is  simply  the  result  of  a  logical  AND  operation  between  all  REMs.  This 
operation  is  performed  within  a  circuit  on  an  FPGA.  Section  3. 3.4. 2  validates  this  concept. 
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0-5.  Select  hopset.  Using  a  random  key,  an  FPGA  circuit  randomly  selects  available 
channels  for  use  as  a  hopset.  This  hopset  is  2,048  hops  in  length.  Section  3. 3. 4. 3  validates 
this  concept. 

0-6.  Send  hopset.  As  in  operation  0-3,  the  radios  use  TOMC  to  distribute  an 
adaptive  hopset  through  the  network.  Hopsets  are  only  valid  within  a  subnet,  i.e.,  each 
subnet  generates  and  uses  a  unique  hopset. 

0-7.  Use  hopset.  Once  all  nodes  within  a  cluster  have  knowledge  of  the  adaptive 
hopset,  the  radios  commence  data  transmission  and  reception.  This  operation  is  the  result 
of  the  preceding  six,  and  functions  just  as  any  typical  frequency  hopping  radio. 

The  system  returns  from  0-7  to  0-1  when  one  of  two  conditions  are  met:  (1)  the 
hopset  expires,  or  (2)  the  spectrum  has  potentially  changed.  A  hopset  expires  if  all  hops 
have  been  used.  While  individual  frequencies  may  be  reused  in  order  to  generate  a  hopset 
of  desired  length,  the  hopset  itself  is  generally  for  one-time  use  only.  However,  if  for  any 
number  of  reasons  (i.e.,  present  environment  is  deemed  benign,  a  new  key  is  unavailable, 
etc.)  a  new  hopset  cannot  be  generated,  the  hopset  is  reusable.  The  spectrum  is  known  to 
have  potentially  changed  if  any  radio  sends  new  REM  data  during  any  operation  other  than 
0-3. 

Of  the  seven  operations  outlined  above,  operations  0-2 , 0-4.  and  0-5  are  implemented 
in  this  research.  These  operations  are  shaded  green  in  Figure  3.1.  Operations  0-1 ,  0-3.  O- 
6,  and  0-7  are  only  explained  in  this  section  and  are  not  validated  any  further  in  this  work. 

3.1.3  Whole  System  Structure. 

The  solution  proposed  in  this  research  is  a  framework  for  accomplishing  the  operations 
described  above.  Specifically,  a  middleware  architecture  for  adaptive  frequency  hopping 
within  cognitive  radio  network  is  presented.  This  architecture  exists  in  a  hybrid  hardware- 
software  system  on  a  WARP  II  FPGA-centric  radio  board.  This  middleware  architecture  is 
shown  in  Figure  3.2. 
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In  Figure  3.2,  there  are  four  primary  components:  the  RF  spectrum  emulator,  a 
database  for  storing  spectrum  map  data,  the  radio  hardware  with  which  the  radio  transmits 
and  receives,  and  the  middleware  architecture  itself.  The  architecture  itself  consists  of  eight 
sub-components:  database  interface,  spectrum  map  parser,  map  distributor,  clustering,  map 
merger/hopset  selector,  secure  hash  chainer,  multicast  communication  layer,  and  hardware 
interface  layer.  Of  these  sub-components,  the  clustering  and  map  merger/hopset  selector 
are  constructed  and  tested  for  this  work  (shown  here  in  green).  The  remaining  pieces  (i.e., 
on-board  spectrum  sensing,  radio  card  operation,  and  network  communication)  are  to  be 
demonstrated  in  future  iterations  of  this  project.  The  overall  system  and  its  operation  is 
proposed  in  Section  A. 

The  clustering  implementation  is  integrated  into  the  architecture  purely  in  software. 
In  order  to  minimize  traffic  within  the  overall  network,  the  radios  exist  within  a  packet- 
switched  network.  Geographically-close  nodes  belong  to  the  same  subnet.  By  using  a 
clustering  algorithm  to  partition  a  collection  of  many  nodes,  it  is  expected  nodes  contained 
within  the  same  subnet  will  have  highly-similar  RF  spectra.  In  Section  B,  an  experiment  is 
presented  for  selecting  the  best  combination  of  number  of  clusters  (equivalent  to  number 
of  subnets)  coupled  with  the  ideal  clustering  algorithm  heuristic. 

The  map  merger/hopset  selector  component  exists  primarily  in  hardware  with  limited 
software  interfacing.  After  accepting  an  indefinite  number  of  maps  from  radios  within  the 
network  and  a  random  key,  it  generates  an  adaptive  hopset  which  maps  to  the  frequencies 
(or  channels)  usable  by  radios  within  the  network.  This  capability  allows  the  system  to  fit 
within  the  available  spectrum  “whitespace”  in  true  DSA  fashion.  This  design  is  validated 
in  Section  C. 
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Figure  3.2:  Proposed  middleware  architecture. 
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3.2  Network  Clustering 

3.2.1  Problem  Definition. 

3.2. 1.1  Goals  and  Hypothesis. 

The  goal  of  this  experiment  is  to  demonstrate  a  method  for  selecting  a  clustering 
algorithm  for  integration  with  the  rest  of  the  middleware  architecture.  The  “defacto”  k- 
means  algorithm  is  Lloyd’s  algorithm  [36,  37],  so  Lloyd’s  algorithm  is  the  default  choice. 
Therefore,  the  purpose  of  this  experiment  is  to  investigate  which  clustering  heuristic 
described  in  Section  2.4.3,  if  any,  performs  best  over  a  range  of  RF  spectrum  environments 
and  center  counts. 

3.2. 1.2  Approach. 

We  introduce  a  method  by  which  the  four  clustering  heuristics  in  Section  2.4.3  can  be 
evaluated  for  expected  best  performance  in  a  cognitive  radio  network.  We  first  establish  a 
baseline  for  clustering  performance  using  “canned”  maps  and  then  use  the  point  distribution 
generator  within  the  kmltest .  exe  software  to  form  “real”  maps  of  four  distribution  types: 
uniform,  Gauss,  and  multi-cluster.  The  points  used  for  these  maps  represent  receivers  in  a 
network,  as  opposed  to  the  transmitters  we  intend  to  avoid. 

Apply  threshold.  In  the  actual  system,  operation  0-1  (see  Figure  3.1)  converts  the 
observed  spectrum  to  a  format  usable  by  the  map  merging  and  hopset  selection  component 
described  in  Section  3.3.  This  conversion  is  also  necessary  to  validate  the  clustering  test 
framework  proposed  in  this  section.  Because  the  spectrum  is  input  as  an  unprocessed, 
emulated  entity  for  this  experiment,  the  threshold  conversion  is  completed  within  the  test 
framework. 

Establish  baseline.  In  order  to  validate  the  performance  of  kmltest .  exe,  we  form 
several  trivial  maps  to  verify  that  the  four  clustering  heuristics  perform  as  expected.  For 
example,  if  there  are  n  discrete  clusters  as  seen  by  the  human  eye,  it  is  expected  n  clusters 
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Figure  3.3:  Baseline  receiver  configurations. 


will  be  formed  as  they  are  geographically  positioned.  The  baseline  receiver  configurations 
(i.e.,  the  canned  maps)  are  shown  in  Figure  3.3. 

Clustering  algorithm  testing.  The  clustering  algorithm  partitions  a  flat  cognitive 
radio  network  into  sub-networks  (analogous  to  IP  subnets).  This  partitioning  decomposes 
from  a  larger  problem  in  which  every  node  can  “talk”  to  every  other  node  (notionally 
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Figure  3.4:  Network  clustering  test  framework  function. 


of  0(n 2)  complexity)  into  a  simpler  problem  (notionally  toward  0(n  login ))  complexity). 
(Note:  Network  complexity  is  not  proven  as  part  of  this  research,  and  is  used  in  this  research 
as  a  method  for  generalizing  expected  trends.)  The  most  desirable  performance  yields  the 
highest  spectrum  similarity  within  a  cluster  such  that  extraneous  traffic  is  minimized  and 
the  adaptive  hopset  is  chosen  from  the  broadest  possible  space. 

Results  evaluation.  Using  the  data  from  this  experiment,  we  implement  one  or 
more  heuristics  for  use  in  the  overall  cognitive  radio  network.  Selection  is  based  on  best 
performance  across  a  range  of  RF  spectrum  environments  and  node  configurations. 

3.2.2  System  Services. 

S-l.  Generate  node  maps.  Given  the  number  of  nodes,  the  number  of  maps  needed, 
and  map  dimensions,  component  C-l  generates  the  desired  number  of  maps  with  the 
desired  quantity  and  size  using  the  specified  distribution.  There  are  no  bounds  on  the 
inputs.  The  generated  maps  can  be  reused  in  future  tests  for  a  consistent  dataset. 
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5-2.  Decode  spectrum.  Spectrum  data  is  originally  formatted  according  to  the 
virtual  device  that  recorded  the  measurement.  Given  the  original  data  and  the  map 
dimensions,  component  C-2  extracts  that  data  to  produce  a  grid  representation  of  the  RF 
spectrum  of  the  same  dimensions  as  the  maps  generated  during  5-7.  The  output  of  this 
component  takes  the  form  of  multiple  discrete  FFTs. 

5-3.  Generate  REMs.  Using  the  data  generated  by  5-2,  component  C-3  transforms 
the  spectrum  data  into  binary  REMs.  These  REMs  are  used  for  calculation  of  the 
performance  metric. 

S-4.  Generate  kmltest.exe  commands.  Clustering  is  performed  in  an  external 
program,  kmltest.exe,  which  implements  the  KMlocal  testbed  described  in  [37]. 
Normally,  the  clustering  heuristic  used  and  the  number  of  centers  are  input  by  the  user 
at  command  line.  To  automate  the  testing  process,  commands  must  be  prefabricated  and 
piped  in  using  system  calls.  This  service  generates  a  text  file  containing  all  commands 
needed  to  run  the  appropriate  test. 

5-5.  Run  map  through  kmltest.exe.  This  service  is  the  core  of  the  experimental 
design,  and  is  therefore  the  Component  Under  Test  (CUT).  The  program  is  run  using 
system  calls  and  input  is  given  via  the  text  file  generated  during  S-4.  The  program  generates 
output  to  a  predefined  text  file. 

5-6.  Parse  cluster  assignments.  When  kmltest.exe  runs,  it  generates  a  report 
containing  the  assignments  of  all  points.  Because  this  information  cannot  be  returned 
directly  to  the  test  framework  without  extensive  modification,  it  is  more  viable  simply  to 
parse  the  output  file  generated  by  the  7- means  software.  Along  with  the  REMs  derived  in 
C-3,  cluster  assignments  are  used  for  calculation  of  the  performance  metric. 

5-7.  Calculate  spectrum  similarity.  The  similarity  calculation  component  computes 
the  performance  metric  described  in  Section  3. 2. 5. 4.  This  component  accepts  the  list  of  all 
REMs  and  the  list  of  cluster  assignments  as  inputs  and  returns  a  single  number. 
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System  Parameters 


Figure  3.5:  Network  clustering  test  framework  structure. 


3.2.3  System  Boundaries. 

The  system  under  test  (SUT)  is  comprised  of  six  components: 

C-l.  Node  map  generator.  This  component  exists  within  the  kmltest.exe  program 
and  can  generate  multidimensional  point  collections  of  various  distribution  types.  This 
research  uses  two-dimensional  maps  over  uniform,  Gaussian,  and  bimodal  Gaussian 
distributions. 

C-2.  Spectrum  decoder.  Data  from  the  is  parsed  from  a  structure  formed  by 
the  DYSE  into  a  structure  more  easily  used  in  determining  spectrum  measurements  to 
individual  locations. 

C-3.  REM  encoder.  Applies  the  threshold  value  and  forms  binary  vectors  for  use  in 
comparing  RF  spectra. 
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C-4.  k-means  clusterer.  Performs  the  clustering  operation  using  one  of  four 
heuristics  and  outputs  the  results  to  a  file. 

C-5.  Assignments  parser.  Reads  the  cluster  assignments  from  the  output  file  and 
generates  a  data  structure  containing  all  nodes  and  their  assignments. 

C-6.  Similarity  calculator.  Calculates  the  similarity  metric  for  all  clusters. 

3.2.4  Workload. 

The  component  under  test  is  affected  by  four  workload  parameters:  the  number  of 
nodes,  the  number  of  maps,  each  map’s  dimensions,  and  the  distribution  type  used  to 
create  the  map.  All  four  parameters  are  input  to  the  map  generator  component  in  the 
lcmltest .  exe  software. 


Table  3.1:  Workload  parameters  and  descriptions. 


Workload  Parameter 

Description 

Number  of  nodes 

Defines  how  many  nodes  should  be  created  for  each 

map.  The  number  of  nodes  in  a  map  remains  constant 

for  all  maps  generated  during  a  given  test. 

Number  of  maps 

Defines  the  number  of  maps  generated  for  a  test. 

Map  dimensions 

Specifies  the  x  and  y  dimensions  of  a  map  such  that 

the  minimal  coordinates  will  be  (1, 1)  and  the  maximal 

coordinates  will  be  (x,y).  When  nodes  are  generated, 

they  are  placed  within  these  bounds. 

Distribution  type 

Selects  the  type  of  distribution  used  to  generate  the 

map  of  nodes. 
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3.2.5  Performance  Metrics. 

This  system’s  performance  metric  is  determined  by  the  similarity  between  observed 
spectra  at  nodes  within  clusters.  This  metric  is  determined  using  the  comparisons  of 
multiple  spectrum  maps.  The  method  for  determining  spectrum  maps  is  presented  first, 
followed  by  an  clustering  affects  a  cognitive  radio  network.  Last,  a  node-wise  similarity 
metric  (the  main  system  performance  metric)  is  examined. 

3.2.5. 1  Spectrum  Representation  and  Threshold  Determination. 

In  this  research,  the  REM  is  represented  as  a  binary  vector.  This  vector  is  derived  from 
an  analog-to-digital  conversion  of  the  analog  spectrum.  In  the  case  of  using  a  discrete  FFT 
to  analyze  the  spectrum,  the  number  of  FFT  bins  that  corresponds  to  one  discrete  channel 
depends  on  the  channel  bandwidth  required  by  the  user.  For  example,  if  there  are  N  FFT 
bins  and  M  discrete  channels  span  the  given  frequency  band,  [  jj]  bins  map  to  each  channel. 
We  assume  that  the  mapping  of  bins  to  channels  and/or  center  frequencies  has  already  been 
established. 


H(co) 


co 
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Figure  3.6:  Threshold  application  example. 
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Using  the  discrete  FFT  output,  we  apply  a  threshold  value  to  each  bin.  In  [35], 
researchers  present  a  simple  method  for  applying  a  threshold  to  a  signal  to  extract 
information.  The  average  of  each  discrete  FFT  bin’s  value  is  multiplied  by  some  coefficient 
to  increase  the  likelihood  that  the  threshold  is  applied  above  the  noise  floor.  In  this  research, 
the  coefficient  is  0.2.  This  translates  to  moving  the  threshold  20%  higher  than  the  average 
FFT  bin  value  when  applied  to  a  linear  scale  (see  Figure  3.13).  We  use  this  threshold 
method  to  evaluate  each  bin  for  energy  presence.  Pilot  experiments  show  that  in  order  to 
eliminate  as  much  noise  as  possible  without  overlooking  actual  spectrum  usage,  this  value 
should  be  held  constant  at  0.2. 

Calculation  and  application  of  the  threshold  is  performed  in  the  spectrum  sensing 
operation  (O-l  in  Figure  3.1)  by  marking  each  bin  as  a  ’1’  if  the  spectrum  is  unused  above 
the  threshold.  Otherwise  the  bin  is  marked  as  a  ’O’.  The  result  of  this  process  is  the  binary 
vector  spectrum  representation  where  ones  represent  available  frequencies  and  zeros  show 
unavailable  spectrum.  An  example  of  the  threshold  conversion  method  is  shown  in  Figure 
3.6,  and  the  threshold  is  applied  to  an  actual  RF  map  in  Figure  3.13. 

3.2.5. 2  REM  Scenarios. 

A  geographically-separated  network  does  not  share  a  common  REM.  A  single 
frequency  hopping  sequence  cannot  be  selected  which  every  node  can  concurrently  use. 
Individual  clusters,  however,  are  more  likely  to  have  common  unoccupied  spectrum.  To 
achieve  this  goal  and  build  the  aforementioned  hierarchy,  the  strategy  of  clustering  is  break 
the  network  into  smaller  clusters,  which  do  share  common  REM.  In  this  manner,  Scenario  C 
distills  into  separate  instances  of  Scenarios  A  and/or  B.  As  shown  in  the  notional  examples 
of  Figures  3.7,  3.8,  and  3.9,  the  greatest  amount  usable  spectrum  is  available  when  clustered 
nodes  have  similar  observed  spectra,  likely  the  result  of  geographical  closeness. 

S-A,  Similar.  Nodes  are  geographically  close  to  each  other  and  have  nearly  identical 
REMs.  A  fused  common  REM  is  approximately  equivalent  to  any  individual  REM  within 
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Figure  3.7:  Cluster  scenario  S-A. 


the  cluster.  This  is  the  simplest  scenario  as  the  aggregate  REM  is  easiest  to  compute,  i.e., 
the  cluster  leader’s  REM.  Any  given  node  can  communicate  with  any  other  node  at  will. 
An  example  of  individual  and  merged  spectra  in  S-A  is  shown  in  Figure  3.7  where  the 
spectrum  is  partitioned  into  32  bins.  Dark  areas  represent  occupied  spectrum  and  light 
areas  represent  usable  “white  space.” 

S-B,  Overlapping.  Nodes  with  moderate  geographical  separation  likely  have  non¬ 
identical,  but  overlapping,  REMs.  A  fused  common  REM  is  equivalent  to  the  intersection 
of  the  individual  REMs.  This  scenario  is  identical  to  A  with  the  caveat  that  the  space  of 
possible  hop  destinations  is  more  limited.  An  example  of  individual  and  merged  spectra 
in  S-B  is  shown  in  Figure  3.8  where  the  spectrum  is  partitioned  into  32  bins.  Dark  areas 
represent  occupied  spectrum  and  light  areas  represent  usable  “white  space.” 


REMl  RE  M2  REM3 


Merge(REM  1 ,  REM2 ,  REM3)  Valid  common  REM 



Figure  3.8:  Cluster  scenario  S-B. 


S-C,  Disjoint. 
overlapping  REMs. 


Geographically-far  nodes  have  non-identical  and  potentially  non- 
In  this  scenario,  no  more  than  two  nodes  share  common  available 
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spectrum  such  that  the  cluster’s  REM  intersection  is  disjoint.  Therefore,  a  node  has  a 
mutually-exclusive  REM  intersection  with  only  one  other  node.  In  terms  of  ad  hoc  network 
communications,  this  is  the  worst  scenario  as  a  gateway  node  would  need  to  relay  packets 
between  adjacent  nodes.  An  example  of  individual  and  merged  spectra  in  S-C  is  shown  in 
Figure  3.9  where  the  spectrum  is  partitioned  into  32  bins.  Dark  areas  represent  occupied 
spectrum  and  light  areas  represent  usable  “white  space.” 
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Figure  3.9:  Cluster  scenario  S-C. 


Therefore,  the  goal  of  clustering  is  to  geographically  partition  a  set  of  nodes  with 
disjoint  REMs  (Scenario  C)  into  clusters  which  have  similar  or  overlapping  (Scenario  A/B) 
REMs  with  the  intent  of  high  intra-cluster  spectrum  similarity.  High  intra-cluster  spectrum 
similarity  implies  a  broader  sample  space  from  which  to  select  an  adaptive  hopset.  A 
broader  sample  space,  then,  is  expected  to  decrease  the  possibility  of  interference. 

3.2. 5. 3  Spectrum  Map  Comparison. 

In  this  research,  REMs  are  represented  as  binary  vectors.  As  such,  it  a  trivial  operation 
to  compare  two  REMs.  Hamming  distance,  a  commonly  used  metric  in  coding  theory, 
counts  the  number  of  coefficients  by  which  two  strings  differ.  If  bits  are  considered  as 
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coefficients  and  REMs  as  strings,  the  Hamming  distance  between  two  REMs  of  equal 
length  is  the  number  of  bits  which  are  different  for  all  equivalent  positions.  An  graphical 
example  of  Hamming  distance  is  shown  in  Figure  3.10.  Positions  at  which  bits  are  equal 
are  shaded  green;  likewise,  unequal  bits  are  shaded  red. 
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Figure  3.10:  Hamming  distance  example. 


Hamming  distance  is  directly  applicable  to  determining  the  effectiveness  of  clustering. 
While  the  original  metric  determines  the  bit-wise  difference  between  two  REMs, 
subtracting  that  number  from  unity  and  dividing  the  difference  by  the  number  of  bits  yields 
the  percent  similarity  between  two  REMs.  The  latter  is  used  to  evaluate  how  effectively 
clustering  groups  together  nodes  with  like  spectra  and  will  be  referred  to  as  “REM  distance” 
from  this  point. 

3.2. 5.4  Intra-Cluster  Spectrum  Similarity. 

Using  the  cluster  scenarios  of  Section  3. 2. 5. 2  and  the  distance  metric  of  Section 
3. 2. 5. 3,  the  test  system’s  sole  metric  is  intra-cluster  spectrum  similarity  (ICSS).  To 
determine  the  overall  spectrum  similarity  within  a  cluster,  the  REM  distance  is  computed 
for  every  pair  of  nodes  within  the  cluster.  No  pair  carries  a  particular  “weight”  over  any 
other  pair,  so  intra-cluster  spectrum  similarity  is  defined  simply  as  the  average  of  all  REM 
distances  within  a  cluster.  An  example  of  ICSS  calculation  is  shown  in  Figure  3.11  where 
ICSS  =  84.5%. 
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ICSS  =  (0.93  +  0.89  +  0.85  +  0.82  +  0.80  +  0.78)/ 6  =  0.845 
Figure  3.11:  ICSS  calculation  example. 


We  do  not  use  Hamming  distance  as  our  experimental  metric  because  it  is  dependent 
on  the  number  of  bits  (which  map  to  the  bins,  channels,  frequencies,  etc.  in  the  observed 
spectrum)  in  the  vectors  being  compared.  ICSS  is  only  a  slight  modification  of  Hamming 
distance  and,  as  a  ratio,  it  is  independent  of  vector  size.  Therefore,  ICSS  can  also  be  used  to 
compare  experiments  using  arbitrary-length  spectrum  vectors.  (Note:  We  do  not  use  ICSS 
for  this  purpose  in  this  paper,  but  instead  build  our  experiment  for  future  expansion.) 

3.2.6  System  Parameters. 

The  system  under  test  accepts  five  system  parameters:  the  set  of  RF  spectrum  maps, 
the  number  of  bins  used  in  the  FFT  (corresponding  to  the  number  of  bits  in  a  REM),  the 
REM  threshold  coefficient,  k- means  heuristic,  and  the  number  of  C means  centers.  These 
parameters  and  the  reasons  they  are  chosen  are  shown  in  Table  3.2. 

The  reason(s)  for  choosing  each  system  parameter  follows: 

•  RF  spectrum  maps.  The  set  of  RF  spectrum  maps  is  a  parameter  because  varying 
the  RF  environment  presents  different  REMs  for  comparison  using  ICSS.  An  RF 
spectrum  map  is  unique  to  each  grid  location  and  contains  the  FFT  of  the  spectrum 
observed  at  that  location.  All  spectrum  map  inputs  are  shown  in  Appendix  A. 
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•  Number  of  FFT  bins.  This  parameter  affects  determination  of  the  REMs  as  the 
number  of  FFT  bins  is  also  the  number  of  bits  in  a  REM. 


•  REM  threshold  coefficient  (a).  Given  that  each  map  contains  a  unique  energy  density 
at  each  bin,  the  REM  threshold  coefficient  sets  the  REM  threshold  (see  Figure  3.6) 
a  percent  above  the  mean  of  the  energy  levels  in  that  spectrum.  The  equation  for 
computing  the  REM  threshold  using  the  threshold  coefficient  a  is  shown  in  3.1, 
where  N  is  the  number  of  bins,  and  S  represents  the  vector  containing  all  N  energy 
density  values.  For  this  research,  spectral  density  values  are  recorded  in  decibels  and 
are  negative  in  sign.  The  value  for  a  is  limited  to  the  range  [-0.5, 0.5]. 


(3.1) 


This  research  uses  a  coefficient  of  0.2  as  in  Figure  3.13,  which  is  similar  to  the 
original  threshold  coefficient  example  in  [35].  If,  in  a  given  bin  of  the  FFT,  the  signal 
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Figure  3.13:  Applied  threshold  coefficient  example. 
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exists  above  the  threshold,  the  corresponding  bit  of  the  resulting  REM  is  marked  as  a 
‘1’  vice  a  ‘0’  otherwise.  This  technique  is  applied  in  Section  3. 2. 5.4  when  computing 
the  performance  metric  and  in  Section  3. 3. 4. 2  for  map  merging. 

•  k-means  heuristic.  Exposing  the  different  heuristics  allows  for  a  decision  on  which 
heuristic,  if  any,  outperforms  the  others  and  should  therefore  be  integrated  with  the 
rest  of  the  middleware  architecture. 

•  Number  of  k-means  centers.  It  is  expected  the  number  of  k-means  centers  affects 
how  many  nodes  will  belong  to  a  cluster,  thereby  influencing  the  commonly-available 
spectrum  throughout. 

•  Node  distribution  standard  deviation  (cr).  The  clustering  software  accepts  a 
number  of  parameters,  including  the  standard  deviation  used  in  generating  several 
distributions  (including  Gauss).  In  general,  increasing  this  value  increases  the  area 
covered  by  the  Gauss  distribution. 
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Table  3.2:  System  parameters. 


System  Parameter 

Description 

RF  spectrum  maps 

Represents  the  maps  used  to  simulate  varying  spec¬ 
trum  environments.  Figures  A.l  and  A. 2  shows  two 

examples  of  ten-by-ten  spectrum  measurement  grids 

where  the  top  left  plot  corresponds  to  position  (1,1) 

and  the  bottom  right  plot  to  (10,10).  Thresholds  are 

denoted  with  a  red  dashed  line. 

Number  of  FFT  bins 

The  number  of  bins  in  an  FFT  of  the  spectrum.  The 

value  used  for  this  experiment  is  2,048  due  to  the 

equipment  configuration. 

REM  threshold  coefficient 

Determines  how  much  the  mean  energy  level  is 

adjusted  when  evaluating  the  spectrum  for  availability. 

A-means  heuristic 

The  different  heuristics  by  which  A- means  clustering 

is  performed.  These  heuristics  are  listed  and  briefly 

described  in  Section  2.4.3. 

Number  of  A'- means  cen¬ 

Defines  how  many  centers  should  be  used  by  the 

ters 

specified  A-means  heuristic. 

Node  distribution  cr 

Provides  a  non-default  standard  deviation  for  deter¬ 
mining  several  distributions.  Pilot  experiments  show 

that  cr  =  0.3  yields  Gauss  distributions  that  consume 

roughly  the  same  area  as  the  other  distribution  types. 
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3.2. 7  Factors. 


Arguably  the  single  most  important  parameter  is  the  RF  spectrum  map  under  which 
a  node  map  exists.  If  the  spectrum  map  is  held  constant,  it  is  useless  to  assert  clustering 
produces  expected  ICSS  results.  Thirteen  total  RF  spectrum  maps  are  used  to  validate 
clustering  algorithm  performance.  Additionally,  four  clustering  heuristics  and  seven  center 
counts  are  used. 


Table  3.3:  System  factors. 


Factor 

Levels 

RF  spectrum  maps 

(see  Figures  3.14  and  3.15) 

Number  of  maps 

10 

Node  distributions 

Uniform,  Gaussian  (//  =  0),  Multi-Cluster  ( k  =  2). 

A- means  heuristic 

Lloyd’s,  Swap,  EZ-Hybrid,  Hybrid 

Number  of  centers 

2  to  98,  in  intervals  of  2  (49  total) 

Ten  RF  spectrum  measurement  grids  are  used  to  validate  the  system’s  performance 
across  a  range  of  transmitter  configurations.  These  spectra  are  emulated  and  their  FFTs 
taken  (“snapshots”  of  the  spectrum)  at  every  whole-number  coordinate  within  the  specified 
dimensions  (i.e.,  (4, 5)).  The  REMs  formed  from  these  snapshots  are  used  in  determining 
ICSS.  All  ten  spectra  are  shown  in  Appendix  A. 

For  the  non-baseline  data,  the  spectrum  at  each  grid  location  is  generated  via  spectrum 
emulation  equipment  where  a  custom  transmitter  map  may  be  specified.  The  transmitter 
maps  used  to  generate  the  RF  spectra  for  this  experiment  are  shown  in  Figures  3.14  and 
3.15.  In  these  diagrams,  the  map  has  dimensions  10  x  10  such  that  each  has  a  maximal 
coordinate  of  (10, 10).  Triangles  represent  transmitters,  large  circles  represent  high-power 
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Figure  3.14:  Transmitter  configurations  1-6. 
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Figure  3.15:  Transmitter  configurations  7-10. 


transmitters  (as  opposed  small  circles  for  low-power),  and  dashed  lines  represent  wide¬ 
band  transmitters  (as  opposed  to  dotted  lines  for  narrow-band). 

We  use  kmltest.exe  to  generate  ten  different  receiver  maps  of  three  distribution 
types.  Each  map  is  pseudo-randomly  generated  using  integer  seed.  Because  the  program  is 
known  to  crash  when  using  seeds  of  three  and  seven,  the  seed  number  is  incremented  twice 
upon  reaching  those  numbers  such  that  they  are  not  used.  One  map  of  each  distribution 
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Figure  3.16:  Sample  uniform  (left),  Gauss  (center),  and  multi-cluster  (right)  node 
distributions  with  seed  =  1 . 


type  (seed  =  1)  is  shown  in  Figure  3.16  where  the  x-  and  y-axcs  correspond  to  the  nodes’ 
respective  physical  locations.  All  ten  maps  (using  seeds  {1,  2,  4,  5,  6,  8,  9,  10,  11,  12})  are 
shown  in  Appendix  B. 

The  /.'-means  testbed  described  in  Section  2.4.3  features  four  clustering  heuristics. 
Therefore,  all  four  are  used  for  this  experiment  for  exploring  any  variability  and/or 
improvement  among  heuristics  prior  to  incorporation  with  the  system. 

All  center  counts  are  powers  of  two  for  binary  division  of  the  map.  We  use  this  rule  to 
facilitate  simpler  creation  of  a  more  vertical  hierarchy  in  future  work.  Additionally,  pilot 
experiments  revealed  that  kmltest .  exe  crashed  when  the  number  of  centers  was  equal  to 
the  number  of  nodes.  Therefore,  the  maximum  number  of  centers  is  the  highest  power  of 
two  less  than  the  number  of  nodes;  in  this  case,  100. 

To  emulate  the  ability  of  nodes  to  organically  sense  the  RF  spectrum,  nodes  are 
mapped  to  DYSE  spectrum  measurements.  When  nodes  are  assigned  spectra  for  clustering 
heuristic  evaluation,  they  receive  the  spectra  assigned  to  the  grid  cell  in  which  they  are 
located,  per  Figure  3.16  and  the  rest  of  the  maps  located  in  Appendix  B.  The  grid  cell 
is  the  whole  number  of  the  node’s  coordinates.  For  example,  a  node  with  a  location 
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Figure  3.17:  Spectrum-to-node  mapping  example. 


(4.23288, 7.31245)  will  be  assigned  the  spectra  of  grid  cell  (4, 7).  A  demonstration  of  this 
mapping  of  spectrum  measurements  to  nodes  is  shown  in  Figure  3.17.  Black  lines  show 
how  spectrum  measurements  map  to  individual  grid  cells,  and  red  arrows  show  how  those 
grid  cells  with  spectrum  measurements  are  mapped  to  the  nodes  within  each  grid  cell. 

3.2.8  Evaluation  Technique. 

3.2.8. 1  Technique. 

Evaluation  is  via  simulation  followed  by  analysis.  It  is  not  practical  to  evaluate  this 
system  purely  by  simulation  because  the  ICSS  metric  must  be  analyzed  using  simulation 
output. 

3.2.8. 2  Experimental  configuration. 

Host  machine.  This  work  is  performed  on  a  laptop  containing  an  Intel  Core  i7- 
2720QM  containing  eight  cores  running  at  2.20  GHz.  The  machine  has  16.0  GB  of  RAM 
and  runs  Windows  7  with  Service  Pack  2. 
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Software  and  code.  MATLAB  is  used  to  build  this  test  framework,  run  the 
simulation,  and  analyze  output.  All  associated  code  is  included  in  Appendix  F. 

3.2.83  Results  validation. 

Results  are  validated  by  examining  (a)  whether  clustering  does  indeed  behave  as 
expected  using  the  baseline  maps,  (b)  whether  the  ICSS  value  does  asymptotically  approach 
unity  as  the  number  of  clusters  approaches  the  number  of  nodes,  and  (c)  whether  the  ICSS 
value  is  consistent  between  distribution  types  and  across  different  node  distributions  and  RF 
maps.  If  the  two  latter  observations  hold  throughout  the  experimental  results,  the  A'- means 
algorithm  is  a  solid  fit  for  use  in  a  network  where  the  network’s  configuration  depends  on 
the  RF  spectrum  environment.  Additionally,  if  one  heuristic  consistently  performs  better 
than  Lloyd’s  algorithm,  then  that  heuristic  is  deemed  the  primary  choice  for  use  in  such  a 
network.  If  no  heuristic  does  so  consistently  or  at  all,  then  Lloyd’s  algorithm  is  chosen  for 
implementation.  This  validation  approach  is  visual  confirmation  based  on  a  large  dataset, 
so  it  is  expected  any  trend  is  distinct  and  readily  identifiable  by  the  human  eye. 

3.2.9  Experimental  Design. 

This  experiment  is  full  factorial  as  the  number  of  possible  input  combinations  is 
feasible  to  implement.  There  are  58,800  total  simulations  (10  RF  spectrum  maps  x  10 
node  distributions  x  3  distribution  types  x  4  clustering  heuristics  x  49  center  counts). 

3.2.10  Methodology  Summary. 

This  experiment  is  designed  to  evaluate  the  effectiveness  of  several  different  clustering 
algorithms  for  incorporation  into  a  cognitive  radio  system.  Experimental  data  and  analysis 
is  generated  using  MATLAB  code,  and  results  are  compared  using  the  proposed  new 
metric,  ICSS. 
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3.3  Adaptive  Hopset  Selection 

3.3.1  Spectrum  Input. 

Spectrum  sensing  is  the  first  step  in  forming  the  adaptive  hopset.  First,  a  radio  front- 
end  receiver  (Rx)  senses  the  spectrum.  An  analog-to-digital  (ADC)  converter  converts  this 
signal  to  its  digital  equivalent.  Next,  a  Fast  Fourier  transform  (FFT)  is  performed  on  the 
digitized  signal  to  generate  an  N-bin  frequency-domain  model  of  the  signal.  The  FFT 
output  is  then  converted  to  a  binary  vector  corresponding  to  the  spectrum  representation 
definition  presented  in  Section  3.2.5. 1.  The  entire  spectrum  input  process  is  shown  in 
Figure  3.18,  where  components  designed  as  part  of  this  section  are  shown  in  red. 

The  portion  of  Figure  3. 1 8  in  dashed  lines  shows  the  REM  distribution  process  and  the 
centralized  REM-based  network  clustering  process.  Once  clusters  exist,  common  REMs 
can  be  formed  for  individual  clusters  and  the  network  as  a  whole  using  the  same  merging 
process.  Clustering  is  not  a  focus  of  this  paper,  however,  and  will  not  be  explained  further 
here. 

In  this  section,  we  assume  spectrum  sensing  is  complete  for  the  network,  that  every 
node  can  represent  its  REM  in  this  manner,  and  that  all  nodes  have  received  the  REMs  of 
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Figure  3.18:  Spectrum  input  diagram. 
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all  nodes  in  the  network.  In  terms  of  Figure  3. 18,  we  assume  that  all  steps  prior  to  the  REM 
Merger  block  are  complete  and  REM  vectors  are  ready  for  merging. 

3.3.2  FPGA  Internal  Structure. 

We  implement  our  design  on  a  Xilinx  Virtex  IV  FPGA  as  part  of  the  WARP  II  radio 
transceiver  board.  Using  the  Xilinx  Embedded  Development  Kit  (EDK),  we  assemble  a 
collection  of  custom  and  existing  intellectual  property  (IP)  cores  using  a  32-bit  bus  and 
the  FPGA’s  embedded  PowerPC  processor.  In  Figure  4,  green  cores  are  affiliated  with  the 
PowerPC,  blue  cores  are  WARP-specific,  tan  cores  are  those  needed  to  support  our  IP  core, 
and  our  core  is  shown  in  red.  Cores  with  an  asterisk  (*)  are  not  used  in  this  experiment,  but 
will  be  used  in  future  development. 

The  IP  cores  supporting  our  core  are  RS-232  (for  output  to  PC  serial  port  terminal)  and 
a  removable  CompactFlash  memory.  The  latter  two  IP  cores  are  grouped  together  by  a  gray 
dashed  fine  because  they  are  interchangeable  as  method  for  exchanging  REMs.  Initially, 
we  simply  write  each  node’s  REM  to  the  removable  memory  and  transfer  REMs  between 
nodes  in  this  fashion. 

Future  iterations  of  this  design  will  include  transmit  REM  data  via  wired  network  for 
system  testing  purposes.  Finally,  radios  will  operate  in  a  completely  wireless  fashion  using 
their  radio  cards.  The  bus  structure  of  our  design  is  shown  in  Figure  3.19. 

3.3.3  IP  Core  Internal  Structure. 

Our  IP  core  consists  of  six  components:  REM  merger,  key  loader,  two  64x32  register 
files  for  storing  the  aggregate  REM  and  the  key,  the  hopset  selector,  and  a  BRAM  section 
in  which  open  channels  are  stored.  Their  organization  is  shown  in  Figure  3.20.  Buses  are 
shown  in  bold,  and  bus  widths  are  in  square  brackets.  A  larger  version  of  3.20  is  shown  in 
Appendix  C  as  Figure  C.2.  All  VHDL  code  necessary  to  implement  this  IP  core  is  included 
in  Appendix  G. 
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Figure  3.19:  FPGA  bus  structure  diagram. 
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Figure  3.20:  AHS  structural  diagram. 


3.3.4  IP  Core  Internal  Function. 

3.3.4. 1  Bandwidth  Masking. 

Different  radio  protocols  have  different  bandwidth  requirements.  In  order  to  provide 
flexibility  with  regard  to  these  requirements,  our  device  accepts  two  32-bit  bandwidth 
masks.  When  the  hardware  scans  the  aggregate  REM  for  available  channels,  the  two  masks 
are  applied  to  determine  whether  the  required  band  is  available.  For  example,  if  a  given 
protocol  requires  B  bins  to  be  available  around  a  center  frequency,  the  user  specifies  the 
corresponding  bits  in  the  bandwidth  mask  input  vectors  as  ’high’.  The  construction  of  the 
bandwidth  mask  is  shown  in  Figure  3.21. 
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Figure  3.21:  Example  bandwidth  mask  vector. 
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3. 3.4.2  REM  Merging  Components. 

REM  merging  is  the  first  core  function  of  this  design.  In  this  paper,  we  assume  that 
each  REM  is  a  2048-bit  binary  vector,  the  direct  result  of  a  2048-bin  FFT.  REMs  can  be 
composed  of  a  variety  of  types  of  spectrum  data;  we  use  only  spectral  data  represented  as 
binary  vectors  per  Section  3.2.5. 1. 

Further,  REM  merging  is  the  process  of  successive  bitwise  AND  operations.  In  this 
manner,  the  AND  of  any  number  of  REMs  produces  a  new  REM  for  which  available 
channels  are  the  intersection  of  all  available  channels  across  all  input  REMs.  Each  vector 
is  loaded  successively  and  ANDed  with  the  previous  result  per  the  merging  technique 
presented  in  Figure  3.22.  There  is  no  physical  limit  to  the  number  of  input  REMs,  although 
the  number  of  possible  REMs  is  ultimately  limited  to  the  network  size.  All  vectors  are 
serially  loaded  at  a  rate  of  one  word  per  clock  cycle.  The  number  of  clock  cycles  needed  to 
load  a  map  is  modeled  by  Equation  3.2. 


#REMs  Bbits  1  cycle 
Cycles  =  - x _ x 


(3.2) 


map  REM  Buswidth 
For  example,  consider  the  scenario  in  which  64  REMs  are  used.  A  64-bit  map  with 

4, 096 cycles  to  load  all  maps.  If  the  number  of 
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Figure  3.22:  Map  merging  example. 
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vectors  is  less  than  the  total  number  of  possible  input  gates,  unused  input  gates  are  assigned 
all  ones  such  that  the  bitwise  AND  of  any  vector  with  that  of  all  ones  results  in  the  original 
vector.  All  vectors  are  serially  loaded  at  a  rate  of  one  word  per  clock  cycle. 

The  REM  merger  can  accept  an  arbitrary  number  of  maps  as  the  hardware  component 
logically  ANDs  each  32-bit  input  with  its  respective  section  in  the  aggregate  REM.  For 
example,  if  the  nth  32-bit  section  of  a  REM  is  loaded  into  the  REM  merger,  it  will  be 
ANDed  with  the  nth  section  of  the  REM  stored  by  the  merger.  The  result  is  then  written  to 
the  register  storing  that  section  of  the  aggregate  REM.  In  this  way,  any  number  of  devices 
in  a  network  can  share  REMs  for  the  purpose  of  creating  an  adaptive  hopset. 

33.4.3  Adaptive  Hopset  Selection. 

Adaptive  hopset  selection,  or  AHS,  represents  the  second  core  function  of  our  design. 
Map  merging  generates  the  aggregate  REM  as  a  binary  vector  in  which  each  bit  maps 
to  the  center  frequency  of  an  FFT  bin  within  the  device’s  usable  spectrum.  Because  all 
radios  within  the  network  know  this  mapping,  adaptive  hopset  selection  becomes  a  matter 
of  randomizing  these  channels.  We  assume  the  hopset  S  is  H  hops  in  length,  bit  position- 
to-channel  mappings  are  stored  in  the  device  (i.e.,  in  a  look-up  table),  and  an  A-bit  random 
key  K  exists  prior  to  forming  the  hopset.  Under  these  assumptions,  our  circuit  performs  the 
following  steps  on  an  A-bit  REM  to  convert  the  map  into  a  corresponding  adaptive  hopset: 
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Next_hop  flag 
asserted 


1.  Starting  at  the  lowest-number  channel,  begin  iterating  through  the  binary  vector. 
Initialize  the  available  channel  count  C  =  0. 

2.  For  every  bit,  if  the  bit  is  a  ‘1’,  look  up  the  corresponding  channel  number  and  store 
it  in  a  list  L  containing  available  channels  and  increment  C. 

3.  IfC  <  H,  replicate  the  first  H  -  C  at  the  end  of  the  list  L  such  that  L  has  H  elements. 

4.  For  i  =  0  to  i  =  H  -  1,  look  up  the  available  channel  in  L  at  the  index  indicated  by 
the  key  K,  scaled  to  the  number  of  open  channels.  Output  the  channel  in  L  as  the  ith 
hop  in  hopset  S .  Rotate  K  one  bit  to  the  right. 
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5.  Repeat  step  four  H  times,  i.e.,  until  all  H  hops  have  been  assigned  within  S . 

The  hopset  S  of  length  N  is  this  system’s  output.  Because  each  index  is  derived  from 
a  random  number,  every  index  is  inherently  random.  This  process  is  graphically  illustrated 
as  a  finite  state  machine  in  Figure  3.23. 

3.3.5  System  Testing. 

VHDL  code  was  developed  in  ModelSim  PE.  Once  we  developed  a  consistent  simula¬ 
tion  design,  we  synthesized  the  hopset  selector  core  in  the  Xilinx  ISE  implementation  tool. 
We  used  the  mapped  version  of  the  circuit  to  confirm  the  design’s  functionality  with  the 
effects  of  realistic  delays  and  other  hardware  characteristics  prior  to  loading  the  integrating 
the  circuit  with  the  rest  of  the  WARP  system.  We  show  the  results  of  this  simulation  in  the 
Results  section. 

We  tested  the  hardware  system  using  C  code  implemented  on  the  FPGA’s  embedded 
processor  core  by  reading  in  the  spectrum  maps  from  a  compact  flash  file,  submitting  the 
maps  to  the  VHDL  code  described  above,  and  writing  the  resulting  hopset(s)  to  a  new 
compact  flash  file  to  simulate  REM  reception  and  hopset  transmission  over  the  network. 
We  use  random  maps  and  a  random  key  to  ensure  a  random  hopset  is  indeed  generated 
from  a  realistic  input. 

3.3.6  Optimization  Goals. 

Virtex  IV  FPGA  resources  are  already  limited,  and  incorporating  the  existing  WARP 
system  causes  the  design  to  be  even  more  crowded.  Therefore,  this  design  is  optimized 
for  area.  Because  the  state  machines  and  other  logic  in  our  IP  core  require  relatively  little 
area  (on  the  order  of  5%  total  resources),  our  primary  optimization  was  to  minimize  the 
effort  needed  to  synthesize  and  implement  memory.  This  strategy  implicitly  attempts  to  use 
memory  components  built  into  the  board.  We  used  a  generic  VHDL  block  RAM  (BRAM) 
structure  given  in  the  XST  User  Guide  for  storing  available  channels  [41]. 
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4  Results 


4.1  Network  Clustering 

4.1.1  Clustering  Baseline. 

In  order  to  verify  that  clustering  does  behave  as  expected,  four  canned  node 
distributions  were  generated  such  that  clusters  were  readily  apparent.  The  A;- means 
algorithm  was  run  on  example  distributions  containing  two,  four,  eight,  and  16  clusters. 
Sixteen  clusters  was  the  upper  limit  as  it  was  assumed  the  trend  would  continue  in  the 
same  fashion  as  the  first  four  greater-than-two  binary  counts.  As  shown  in  Figure  4.1,  each 
heuristic  clustered  the  distributions  as  expected  (i.e.,  two  clusters  were  formed  in  the  map 


2  centers  4  centers  8  centers  16  centers 


Figure  4.1:  Canned  cluster  verification  test. 
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containing  two  distinct  clusters,  etc.).  Different  line  colors  are  irrelevant  to  this  portion  of 
the  experiment. 

This  proof-of-concept  validates  the  simple  (but  very  important)  hypothesis  that  the 
clustering  testbed  software  does  perform  as  expected.  As  a  result,  it  is  assumed  all  results 
obtained  using  the  clustering  software  are  accurate  representations  such  that  any  obscurities 
in  the  results  pertain  solely  to  the  spectrum-based  clustering  metric. 

4.1.2  Clustering  Visualization. 

Lines  are  drawn  between  all  nodes  within  each  cluster  to  represent  potential  (pre¬ 
routing)  communication  links.  When  nodes  within  a  cluster  have  highly  similar  REMs 
(i.e.,  a  high  ICSS  value),  the  connecting  links  are  colored  red.  Likewise,  links  between 
nodes  with  dissimilar  REMs  are  colored  blue.  This  is  accomplished  by  selecting  colors 
from  the  MATLAB  colormapO  function  according  to  the  inter-node  ICSS  value.  The 
clustering  visualization  method  show  only  center  counts  of  two,  four,  eight,  16,  32,  and  64 
in  the  interest  of  preserving  space.  Figures  4.2,  4.2,  and  4.2  show  examples  of  this  cluster 
visualization  method  using  the  third  DYSE  map  on  all  three  distribution  types. 
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Figure  4.2:  Cluster  visualization  of  a  uniform  distribution  using  DYSE  map  #3. 
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EZ-Hybrid  Hybrid  Swap  Lloyd  EZ-Hybrid  Hybrid  Swap  Lloyd 
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Figure  4.3:  Cluster  visualization  of  a  Gauss  distribution  using  DYSE  map  #3. 
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Figure  4.4:  Cluster  visualization  of  a  multi-cluster  distribution  using  DYSE  map  #3. 


Multiple  cluster  visualizations  are  presented  in  Appendix  D.  Each  distribution  type 
is  used  ten  times,  and  each  pseudo-random  distribution  (based  on  a  random  seed)  is  used 
three  times.  This  configuration  yields  a  total  of  thirty  different  cluster  visualizations. 
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4.1.3  ICSS  Evaluation. 


ICSS  plots  for  the  first  and  fifth  DYSE  maps  are  shown  for  all  three  distributions  in 
Figures  4.5,  4.6,  and  4.7.  ICSS  is  found  to  vary  between  approximately  60%  and  100% 
across  the  range  of  clustering  heuristics  and  distribution  types.  Each  datapoint  represents 
running  the  relevant  clustering  heuristic  on  the  ten  distributions  shown  in  Appendix  B  and 
averaging  the  ICSS  result.  It  should  be  noted  that  while  the  trends  do  look  nearly  identical, 
close  inspection  reveals  a  slight  yet  distinct  difference  between  plots. 

Five  characteristics  are  readily  apparent  in  all  three  ICSS  plots.  An  asymptotic 
relationship  to  unity  at  high  cluster  counts,  a  dip  at  lower  counts  in  uniform  and  Gauss 
distributions,  the  inconsistency  of  the  EZ-hybrid  heuristic  at  cluster  counts  over  50  in  a 
100-node  distribution,  the  relative  consistency  of  both  the  Lloyd  and  Hybrid  heuristics 
over  (nearly)  the  entire  range,  and  the  immediate  drop-off  in  ICSS  for  the  Swap  heuristic 
at  high  cluster  counts  are  all  plainly  visible  on  first  inspection.  These  characteristics  are 
examined  below. 

The  asymptotic  relationship  over  the  range  of  cluster  values  is  as  expected  and 
confirms  the  hypothesis  that  spectrum  similarity  between  nodes  within  the  same  cluster 
asymptotically  approaches  unity  as  the  mapping  of  nodes  to  clusters  becomes  one-to- 
one.  In  addition  to  confirming  this  trend,  the  radio  can  also  be  pre-loaded  with  a  range 
of  clustering  results  such  that  different  cluster  counts  are  available  for  selection  based 
on  desired  ICSS.  For  example,  if  the  RF  environment  contains  a  significant  amount  of 
interference,  nodes  likely  need  a  higher  ICSS  value  in  order  to  increase  the  likelihood  that 
a  common  hopset  can  be  found.  Per  the  results  of  evaluating  the  ICSS  metric,  this  need 
translates  to  more  clusters.  Conversely,  an  RF  environment  with  little  interference  requires 
less  commonality  (i.e.,  lower  ICSS)  for  approximately  the  same  number  of  common 
frequencies. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim1_Grid.mat,  uniform  distribution) 
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Figure  4.5:  ICSS  for  uniform  distributions  using  DYSE  maps  #1  and  #5. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim1_Grid  mat.  gauss  distribution) 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim5_Grid.mat,  gauss  distribution) 


Figure  4.6:  ICSS  for  Gauss  distributions  using  DYSE  maps  #1  and  #5. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim1_Grid.mat,  multi_clus  distribution) 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim5_Grid.mat,  multi_clus  distribution) 


Figure  4.7:  ICSS  for  uniform  distributions  using  DYSE  maps  #1  and  #5. 
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There  exists  a  dip  in  ICSS  occurs  around  16  clusters,  with  the  ICSS  value  remaining 
below  that  of  two  clusters  until  the  number  of  clusters  reaches  approximately  25.  More 
testing  is  needed  to  explore  whether  ICSS  “breaks  even”  at  25  (in  this  case,  25%  of  the 
total  number  of  nodes),  and  why  this  occurs  in  the  uniform  and  Gauss  distributions,  but  not 
in  the  multi-cluster  distribution. 

The  EZ-hybrid  heuristic  experiences  a  series  of  oscillations  at  higher  cluster  counts 
(i.e.,  those  counts  greater  than  half  the  total  number  of  nodes).  This  is  contrary  to  the 
monotonic  trend  in  the  three  other  heuristics  (aside  from  Swap  at  high  counts).  This  is 
because  EZ-hybrid  is  a  simple  hybrid  of  the  Lloyd  and  Swap  heuristics  and  does  not  make 
any  effort  to  optimize  clustering.  The  Hybrid  heuristic,  on  the  other  hand,  does  not  show 
any  such  dropoff  except  at  the  highest  cluster  count.  This  result,  however,  is  because  the 
test  program  crashes  with  regularity  at  that  cluster  count. 

Similarly,  the  Swap  heuristic  crashes  for  each  of  the  last  four  cluster  counts.  Because 
of  this  instability,  the  Swap  heuristic  appears  to  be  a  poor  choice  for  implementation  with 
the  remainder  of  the  system.  If  it  not  for  the  consistent  crash  tendency  at  high  cluster  counts, 
Swap  performs  evenly  with  (and  often  better  than)  the  three  other  heuristics.  However, 
if  the  desired  ICSS  is  lower  than  approximately  97%  (a  very  likely  scenario),  the  Swap 
heuristic  is  an  ideal  choice. 

4.1.4  Application  to  System  Implementation. 

This  experiment  highlights  the  tradeoffs  present  in  selecting  the  number  of  clusters 
chosen  for  a  CRN  implementation.  From  Figures  4.5,  4.6,  4.7,  and  those  in  Appendix 
E,  it  is  apparent  that  increasing  the  number  of  clusters  yields  higher  bandwidth.  Further, 
partitioning  the  nodes  into  many  clusters  decreases  network  complexity  and  limits  excess 
traffic  loads  on  the  overall  network  as  intra-cluster  communication  is  independent  of 
other  clusters.  However,  while  high  bandwidth  potentially  means  less  interference,  it 
is  not  always  necessary.  Networks  experiencing  few  or  no  occupied  channels  require 
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less  commonality  between  nodes,  which  leads  to  a  lower  ICSS  value  and  fewer  clusters. 
Therefore,  these  results  can  be  used  to  retrieve  in  real-time  a  desirable  cluster  count  given 
a  minimum  or  maximum  ICSS  value. 

The  cluster  count  used  as  input  to  the  clustering  algorithm  in  actual  operation  is  driven 
solely  by  the  desired  ICSS  value,  but  ICSS  is  affected  by  the  radio  environment  itself.  We 
intend  to  incorporate  the  results  of  this  research  as  a  cache  of  experimentally-obtained 
data.  A  cache  reduces  computing  time,  and  it  is  expected  that  the  data  obtained  using 
100  nodes  can  be  extrapolated  to  both  smaller  and  larger  networks  with  the  same  result. 
Ultimately,  using  a  table-lookup  scheme  balances  the  tradeoffs  anticipated  in  different  RF 
environments  where  greater  interference  can  be  offset  by  increasing  the  number  of  clusters 
in  order  to  improve  the  similarity  of  spectrum  observations  between  nodes. 
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4.2  Adaptive  Hopset  Selection 

4.2.1  Mapped  Simulation. 

In  order  to  validate  our  design  before  running  the  time-intensive  place-and-route 
process  and  to  facilitate  quicker  debugging,  we  simulated  the  synthesized/mapped  version 
of  our  circuit.  Figures  4.8-4.13  detail  this  simulation  and  reinforce  our  assertion  that  the 
circuit  performs  as  intended  given  realistic  hardware  constraints.  Signals  of  interest  are 
described  in  captions.  The  synthesized  component  was  instantiated  in  the  same  VHDL  test 
bench  used  to  evaluate  our  code  base.  Xilinx  automatically  generated  the  VHDL  file  using 
Unisim  components  which  accurately  model  hardware  constraints  in  simulation. 

Initially,  the  new_rems  and  new_key  flags  are  toggled  to  clear  the  existing  aggregate 
REM  and  key.  The  REMs  are  loaded  by  writing  the  32  bits  of  REM  data  and  toggling  the 
load_rem  flag  (Figure  4.8). 
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Figure  4.8:  REM  loading 


The  key  is  loaded  in  a  similar  fashion  to  the  REM(s)  with  the  load_key  flag  being 
toggled  (Figure  4.9). 
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Figure  4.9:  Key  loading 


Once  the  software  has  loaded  the  REMs  and  key,  it  toggles  the  start  flag.  The 
circuit  immediately  begins  counting  and  tracking  the  open  channels  by  first  deasserting  the 
finished  flag  (Figure  4.10). 
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Figure  4.10:  Channel  counting 
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After  the  open  channels  have  been  counted,  the  hopset_ready  flag  is  asserted.  The 
software  then  toggles  the  next  Jhop  flag  in  order  to  retrieve  the  hopset  contents  (Figure 
4.11). 
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Figure  4.1 1:  Hopset  retrieval 


The  finished  flag  is  asserted  after  all  hops  have  been  extracted.  The  hopset_ready 
flag  is  deasserted  (Figure  4.12). 


Ijf.  next  hop 

0 

1-ft  hopset  ready 

1 

^finished 

0 

hop_number[10:0) 
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hop_channel[10:0] 

673 
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10000  ps 

...  X  2044  X  2045  x  2046X20 

97"  ::1915  600  :  2045  ;  6 


Figure  4.12:  Hopset  generation  finished 
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Figure  4.13  shows  full  system  simulation.  Hopset  output  occupies  the  most  time 


because  its  speed  is  governed  by  that  of  the  software  running  on  the  processor. 


loadkey  0 

newrems  0 

newkey  0 

inputl  3 1 : 0  ]  7fffffff 

key 131:0)  lce54ebl 

lower_bw_mask|31:0)  80O0O0OO 

upper_bw_mask(31:0)  00000001 

start  0 

next_hop  0 

hopsetready  0 

finished  1 

hop_number[10:0)  0 

hop  channel[10:0]  2029 


Figure  4.13:  Full  operation 


4.2.2  Standalone  Device  Usage. 

When  placed  and  routed,  the  AHS  (including  a  bus  communication  wrapper)  occupies 
approximately  10%  of  the  board’s  slice  resources.  Because  BRAM  is  used,  the  associated 
slice  usage  is  limited  to  a  small  footprint.  Two  BRAMs  out  of  376  total  units  are  used 
to  create  the  open  channels  table.  Table  4.1  provides  a  comparison  between  device  usage 
summaries  for  both  the  individual  IP  core  and  the  WARP-based  system  as  a  whole. 

4.2.3  Timing  Analysis. 

While  the  maximum  frequency  given  by  XST  is  only  70.442  MHz,  the  actual  system 
can  be  run  at  a  much  higher  frequency  because  of  the  automatic  placement  of  clock  dividers 
around  the  chip.  For  example,  we  run  the  PowerPC  clock  at  125  MHz  and  the  bus  clock  at 
100  MHz  with  no  adverse  effects.  This  is  because  the  relatively  low  frequency  applies  only 
to  a  subset  of  clock  nets. 
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Table  4.1:  Device  Resource  Usage  Summary. 


Resource 

Available 

AHS 

AHS + WARP 

Number  of  BUFGs 

32 

1  (3%) 

11  (34%) 

Number  of  External  IOBs 

768 

256  (33%) 

462  (69%) 

Number  of  RAMB  16s 

376 

2  (1%) 

59  (15%) 

Number  of  Slices 

42,176 

4,614(10%) 

15,032  (35%) 

Number  of  SLICEMs 

21,088 

768  (3%) 

1,119(5%) 

4.2.4  Hopset  Selection  Demonstration. 

In  addition  to  actually  building  the  AHS  system  and  implementing  on  an  FPGA,  C 
code  was  written  for  testing  and  verifying  the  system’s  functionality.  Figure  4.14  shows 
the  output  of  the  test.  The  example  shown  is  the  result  of  a  64-bit  implementation.  The 
smaller  size  (as  opposed  to  the  2,048-bit  full  implementation)  is  used  simply  to  aid  in 
demonstration.  Using  all  2,048  bits  would  simply  create  a  large,  complicated  output.  The 
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Figure  4.14:  Hopset  output  example. 
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Figure  4.15:  AHS+WARP  device  usage  diagram. 
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hopset  shown  is  the  result  of  a  64-bit  random  key.  The  code  used  to  produce  this  output 
forms  the  skeleton  onto  which  the  rest  of  the  system  is  to  be  constructed. 

4.2.5  Device  Usage  Floorplans. 

The  PlanAhead  diagrams  in  Figures  4.15  and  C.l  show  how  the  slices  (shown  as  light 
blue  dots)  fill  the  FPGA  fabric  with  the  existing  WARP  system  and  the  AHS  IP  core  and 
IP  core  alone,  respectively.  The  most  notable  difference  is  the  addition  of  routed  PowerPC 
cores  in  the  larger  design  (shown  as  two  dim  blue  blocks). 

4.2.6  Optimization  Achievements. 

Directing  the  synthesis  and  implementation  tools  to  build  the  memory  structure  out  of 
logic  slices  ran  and  did  not  complete,  even  when  given  over  12  hours  to  run.  Purposefully 
using  BRAM  and  using  the  XST  option  for  automatically  extracting  BRAM  was  quick  and 
resulted  in  virtually  no  synthesis  and  implementation  overhead. 
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5  Conclusions 


5.1  Research  Contributions 

This  research  makes  five  contributions  to  the  field  of  cognitive  radio: 

1.  The  framework  for  a  new  hybrid  hardware/software  middleware  architecture  (see 
Figure  3.1); 

2.  A  framework  for  testing  and  evaluating  clustering  algorithms  in  the  context  of 
cognitive  radio  networks  (see  Figure  3.1,  0-2  and  Figure  3.5); 

3.  A  new  RF  spectrum  map  representation  technique  (see  Figure  3.6); 

4.  A  new  RF  spectrum  map  merging  technique  (see  Figure  3.1,  0-4  and  see  Figure 
3.22); 

5.  A  new  method  for  generating  a  random,  key-based  adaptive  hopset  frequency 
hopping  waveform  (see  Figure  3.1,  0-5  and  see  Figure  3.23);  and 

6.  Initial  integration  testing  toward  implementing  the  proposed  system  on  a  field- 
programmable  gate  array  (FPGA)  (see  Figure  4.14). 

5.2  Whole  System 

We  believe  it  is  both  possible  and  feasible  to  implement  an  adaptive  frequency  hopping 
cognitive  radio  for  use  among  existing  primary  and  secondary  users.  To  support  this 
assertion,  we  propose  a  new  middleware  architecture  for  use  in  such  a  system.  We  also 
implement  several  components  as  the  basis  for  proving  feasibility.  Through  defining  the 
architecture  and  implementing  several  operations  within  the  system,  we  conclude  such  a 
system  is  both  possible  and  feasible  to  build. 
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5.3  Network  Clustering 

We  presented  a  methodology  for  evaluating  the  k- means  clustering  algorithm  in 
a  cognitive  radio  network  over  a  range  of  node  distributions,  cluster  counts,  and  RF 
spectrum  maps.  We  also  introduced  a  new  metric,  intra-cluster  spectrum  similarity  (ICSS), 
for  comparing  the  effectiveness  of  clustering  in  the  context  of  a  dynamic  RF  spectrum 
environment.  Our  experiment  showed  that  ICSS  asymptotically  approaches  unity  as  the 
number  of  clusters  in  a  notional  cognitive  radio  network  of  100  nodes  approaches  the 
number  of  nodes.  This  trend  was  as  expected  and  remained  consistent  across  all  ten  RF 
spectrum  maps,  distributions,  distribution  types,  and  clustering  heuristics  for  a  range  of 
cluster  counts.  Because  of  its  consistent  performance,  Lloyd’s  algorithm  (the  de  facto 
heuristic)  was  selected  for  integration  with  the  whole  system. 

In  addition  to  selecting  an  algorithm,  the  consistency  of  results  allows  the  design  to 
incorporate  the  ICSS  results  in  a  table-lookup  format.  For  example,  if  the  radio  determines 
that  spectra  must  be  K-pcrccnt  similar,  the  radio  runs  the  clustering  algorithm  using  X 
clusters.  Given  the  requirement  for  Y,  X  is  easily  derived  from  the  plotted  ICSS  data. 
Further,  if  it  is  determined  a  cluster  must  have  high  bandwidth,  then  a  higher  ICSS  value 
(X)  is  chosen,  and  therefore  a  large  number  of  clusters  (I7)  is  chosen.  Likewise,  a  low- 
bandwidth  connection  requires  a  lower  ICSS  value  and  therefore  fewer  clusters.  Within 
this  decision,  there  exists  a  tradeoff.  High  bandwidth  leads  to  more  clusters  which  leads  to 
decreased  complexity  (toward  0(n  login )),  whereas  low  bandwidth  leads  to  fewer  clusters 
which  leads  to  increased  complexity  (toward  O(n)).  We  conclude  that  the  decision  for  Y 
is  contingent  upon  a  tradeoff  between  ICSS  (representing  similarity  between  REMs)  and 
required  bandwidth,  and  that  the  number  of  clusters  is  not  dependent  on  ICSS  alone,  but 
also  network  requirements. 
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5.4  Adaptive  Hopset  Selection 

The  last  three  contributions  were  made  in  adaptive  hopset  selection.  First,  we 
implemented  a  cognitive  radio  with  adaptive  hopset  selection  on  an  FPGA.  Second, 
we  proposed  and  demonstrated  a  new  technique  for  merging  RF  spectrum  maps  when 
represented  as  binary  vectors.  Finally,  we  proposed  and  demonstrated  a  new  adaptive 
hopset  selection  technique.  The  whole  system  (composed  of  the  WARP  architecture  and 
our  custom  core)  fit  within  the  resources  available  on  the  Virtex  IV  FPGA.  Because  our 
code  simulation,  mapped  hardware  simulation,  and  software -based  testing  yielded  the 
desired  results,  this  portion  of  the  whole  system  functioned  as  expected. 

5.5  Final  Remarks 

The  architecture  proposed  in  this  research  is  a  novel  method  by  which  to  implement 
a  frequency  hopping  cognitive  radio  network  for  coexistence  with  other  RF  spectrum 
users.  By  demonstrating  the  practicality  of  using  the  /.-means  clustering  algorithm  and  RF 
spectrum  measurements  as  an  effective  method  for  partitioning  an  otherwise-fiat  network 
into  sub-networks,  we  provide  an  efficient  way  to  decrease  network  complexity  in  a 
way  that  accounts  for  RF  spectrum  differences.  Further,  by  demonstrating  a  hardware- 
based  implementation  for  selecting  a  random  frequency  hopping  hopset,  we  introduce  an 
avenue  for  adaptively  hopping  frequencies  as  a  means  to  quickly  and  continually  avoid 
RF  interference.  To  support  these  two  contributions,  we  propose  a  lightweight,  hardware- 
portable  spectrum  representation  technique.  It  is  expected  these  contributions  forge  a  path 
for  building  on  the  proposed  architecture  and  ultimately  developing,  testing,  and  fielding 
an  operational  frequency  hopping  cognitive  radio  network. 


94 


6  Future  Work 


As  stated  in  the  Introduction,  several  pieces  within  this  research  remain  to  be 
implemented  before  a  full-scale,  standalone  prototype  network  can  exist.  At  the  time 
of  publication,  one  post-doctoral  student,  two  masters  students,  and  several  interns  are 
currently  working  on  expanding  the  envelope  of  this  research. 

1.  Multicast  communication  layer.  Multicast  communication,  specifically  totally- 
ordered  multicast  (TOMC),  enables  a  conceivably  large  and  complex  network  to 
communicate  in  a  structured  fashion.  Integrating  the  TOMC  layer  onto  the  board 
means  tying  in  existing  software.  Currently,  using  the  Spread  API  is  yielding 
promising  results.  This  layer  will  be  fully  validated  once  nodes  can  communicate 
over  either  wired  or  wireless  links.  This  layer  requires  some  sort  of  MAC  layer  to  also 
be  in  place,  so  it  is  possible  small  and/or  embedded  operating  system  (i.e.,  TinyCore 
Linux,  Windows  Embedded,  etc.)  may  need  to  run  on  the  embedded  processor.  At 
time  of  print,  this  step  is  being  completed  as  part  of  a  follow-on  student’s  thesis. 

2.  Ethernet  core.  An  Ethernet  core  is  the  next  step  toward  inter-node  communication. 
When  implemented,  Ethernet  allows  two  or  more  boards  to  be  a  part  of  the  same 
wired  network  in  much  the  same  way  as  the  eventual  wireless  network.  At  time  of 
print,  wireless  transmission  via  Ethernet  has  been  demonstrated  but  has  not  yet  been 
integrated  with  the  middleware  architecture. 

3.  Frequency  hopping  operation.  The  WARP  board  supports  up  to  four  radio  cards, 
and  these  cards  must  be  integrated  with  the  hopset  selection  method  described  in 
Section  3.3.  Current  work  on  this  step  (outside  the  scope  of  this  document)  is 
demonstrating  the  WARP’s  wireless  functionality,  so  integration  with  the  proposed 
middleware  architecture  is  pending.  In  order  for  the  radio  to  operate  as  intended, 
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though,  this  task  is  extremely  critical.  In  addition  to  implementing  frequency 
hopping-spread  spectrum  communication,  it  is  desirable  to  include  non-contiguous 
OFDM  where  carriers  are  adaptively  selected  using  REM  data. 

4.  Secure  hash  chaining.  To  realistically  generate  new  hopsets  periodically,  keys  need 
to  be  distributed  securely  and  quickly.  Secure  hash  chaining  provides  a  method 
by  which  such  values  can  be  distributed.  This  task  is  intended  for  software-only 
implementation,  but  it  is  feasible  the  operation  could  be  completed  much  quicker  in 
a  hardware -based  FPGA  solution. 

5.  On-board  spectrum  sensing.  Once  the  board  can  sense  its  own  spectrum  and 
produce  a  model  satisfying  the  previously  requirements  laid  out  and  implemented, 
each  board  can  begin  to  function  as  a  fully  autonomous  device.  Spectrum  sensing 
is  currently  performed  in  an  emulation  environment;  this  is  for  prototyping  only. 
Simply  put,  on-board  sensing  completes  the  OODA  loop. 

6.  Field  testing.  Lab  testing  and  emulation  environments  only  validate  this  research 
to  the  accuracy  of  such  tests.  In  order  to  fully  prove  this  system’s  worth,  it  must  be 
tested  in  an  outdoor  environment  where  the  spectrum  is  truly  unknown  and  the  radio 
must  truly  adapt.  This  also  allows  the  radio  to  grow  as  necessary  without  the  physical 
limitations  of  doing  so  indoors. 
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Appendix  A:  DYSE-Generated  RF  Spectrum  Maps 


The  following  plots  represent  “snapshots”  of  various  RF  spectrum  environments. 
These  plots  were  generated  using  AFRL’s  DYSE  system.  Dashed  lines  represent  the 
threshold  applied  to  each  plot. 
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Figure  A.l:  DYSE-generated  RF  spectrum  map  #1. 
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Figure  A. 2:  DYSE-generated  RF  spectrum  map  #2. 


Figure  A. 3:  DYSE-generated  RF  spectrum  map  #3. 
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Figure  A.4:  DYSE-generated  RF  spectrum  map  #4. 


Figure  A. 5:  DYSE-generated  RF  spectrum  map  #5. 
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Figure  A. 6:  DYSE-generated  RF  spectrum  map  #6. 
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Figure  A. 8:  DYSE-generated  RF  spectrum  map  #8. 
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Figure  A. 9:  DYSE-generated  RF  spectrum  map  #9. 


107 


1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047 


ii'^'P*Ynhriw 


rrw^^iw 


•nmrw1^ 

fewiiw 

Tfrwpr^^ 

’1nyFTH|rv'  '  HrpiF'Fj'ii 

'  ’Wifwii 

"fwi* 


1023  20 

[ 

7wT"^,j  ",w"rr 

NWW*1^^ 

1023  2047 


^1/y^rrf  ■ 


. . ' 


1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047  0  1023  2047 


Figure  A.  10:  DYSE-generated  RE  spectrum  map  #10. 
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Appendix  B:  Node  Distributions 


Figure  B.l:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  1). 


Figure  B.2:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  2). 
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Figure  B.3:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  4). 


Figure  B.4:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  5). 


Figure  B.5:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  6). 
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Figure  B.6:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  7). 


Figure  B.7:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  8). 


Figure  B.8:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  9). 
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Figure  B.9:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  10). 


Figure  B.10:  Uniform  (1),  Gauss  (c),  and  multi-cluster  (r)  distributions  (seed  =  11). 
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Appendix  C:  Additional  FPGA  Design  Figures 


Figure  C.l:  AHS  device  usage  diagram. 
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Figure  C.2:  AHS  internal  structure. 
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Appendix  D:  Additional  Clustering  Visualization  Plots 
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Figure  D.l:  Baseline  node  distribution  with  two  apparent  clusters. 
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Figure  D.2:  Baseline  node  distribution  with  four  apparent  clusters. 
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Figure  D.3:  Baseline  node  distribution  with  eight  apparent  clusters. 
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Figure  D.4:  Baseline  node  distribution  with  16  apparent  clusters. 
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Figure  D.5:  Uniform  node  distribution  using  DYSE  map  #1. 
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Figure  D.6:  Gauss  node  distribution  using  DYSE  map  #E 
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Figure  D.7:  Multi-cluster  node  distribution  using  DYSE  map  #1. 
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Figure  D.8:  Uniform  node  distribution  using  DYSE  map  #2. 
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Figure  D.9:  Gauss  node  distribution  using  DYSE  map  #2. 
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Figure  D.10:  Multi-cluster  node  distribution  using  DYSE  map  #2. 


119 


EZ-Hybrid  Hybrid  Swap  Lloyd  EZ-Hybrid  Hybrid  Swap  Lloyd 


2  centers  4  centers  8  centers  16  centers  32  centers  64  centers 


IP 

'Ups#* 

-  •  ■^>  N. 

x* : ; 

■  S  V  -"V  v 

4§? 

Jpt 

K, 

.  ■  ’  ^ 

•  v  y.  'vT; 

•  ’  _  •  *'  i  • 

3P 

jm  Ms- . 

A  ^  ^ 

-A  — 

/ ;  i  v,.  r 
-  -  :■ 
A  ' SJ 

.1  /  •  4  ;  •  ‘ 

I  ' 

' ; 

-MsC4, 

Cv  ^  - 

-.•  A  ^  -=^3 

V  ■  \ 

’  ■  ’  ’  >  -  . 

*  \  y  - 

■  ■■  A»  $ 

\  • . .  •'  ' 

Figure  D.ll:  Uniform  node  distribution  using  DYSE  map  #3. 
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Figure  D.12:  Gauss  node  distribution  using  DYSE  map  #3. 
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Figure  D.13:  Multi-cluster  node  distribution  using  DYSE  map  #3. 
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Figure  D.14:  Uniform  node  distribution  using  DYSE  map  #4. 
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Figure  D.15:  Gauss  node  distribution  using  DYSE  map  #4. 
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Figure  D.16:  Multi-cluster  node  distribution  using  DYSE  map  #4. 
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Figure  D.17:  Uniform  node  distribution  using  DYSE  map  #5. 
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Figure  D.18:  Gauss  node  distribution  using  DYSE  map  #5. 
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Figure  D.19:  Multi-cluster  node  distribution  using  DYSE  map  #5. 
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Figure  D.20:  Uniform  node  distribution  using  DYSE  map  #6. 
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Figure  D.21:  Gauss  node  distribution  using  DYSE  map  #6. 
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Figure  D.22:  Multi-cluster  node  distribution  using  DYSE  map  #6. 


125 


EZ-Hybrid  Hybrid  Swap  Lloyd  EZ-Hybrid  Hybrid  Swap  Lloyd 


2  centers  4  centers  8  centers  16  centers  32  centers  64  centers 


-mm 

m ^  /lot 

^  -X>"  K 

L  'k-'  ^ 

k  ^ 

-=s3l7 

.  A",) . 7  & 

•  ■  >  •■  . 
‘<\  -.  ._•• 

v.v  .  —  A  . 
•  s  ".  . 

‘  7  '  *  ;.  ‘  •• 

nn 

<r-  ^*fe 

.  A>-  .,  T  ■’  A  . 

<k  h  .S'  :  ,  ' 

"A?  \  V 

^  -  5  > 

l> 

‘  V  :.  ^  f.\  ~.' 

iSlS? 

/7L 

-s^lfeu  ■<sSfe: 

*7  % 

0>  'V 

^  ^ 

<r  A  $  -Ote 

/  V  '  h  \ 

\  >  ■  :• 

4  l> 

l  ’  ;.'.• 

A\ 

i§5 

^  ^ 

<r  A  ^ 

■  7  A1'  «  f  ■’  A  ■ 

•  E  ,**■  X  ;  v  " 

V  A  V 

a'  -  ^  " 

<r  A  ^  '='  • 

". 1  •  ■4  .■  ~  ,' 

7  \' -7  'V" 

Figure  D.23:  Uniform  node  distribution  using  DYSE  map  #7. 
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Figure  D.24:  Gauss  node  distribution  using  DYSE  map  #7. 
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Figure  D.25:  Multi-cluster  node  distribution  using  DYSE  map  #7. 
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Figure  D.26:  Uniform  node  distribution  using  DYSE  map  #8. 
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Figure  D.27:  Gauss  node  distribution  using  DYSE  map  #8. 
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Figure  D.28:  Multi-cluster  node  distribution  using  DYSE  map  #8. 
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Figure  D.29:  Uniform  node  distribution  using  DYSE  map  #9. 
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Figure  D.30:  Gauss  node  distribution  using  DYSE  map  #9. 
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Figure  D.31:  Multi-cluster  node  distribution  using  DYSE  map  #9. 
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Figure  D.32:  Uniform  node  distribution  using  DYSE  map  #10. 
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Figure  D.33:  Gauss  node  distribution  using  DYSE  map  #10. 
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Figure  D.34:  Multi-cluster  node  distribution  using  DYSE  map  #10. 
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Appendix  E:  Additional  ICSS  Plots 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  SimIJSrid.mat,  uniform  distribution) 


Figure  E.l:  ICSS  for  uniform  distributions  using  DYSE  map  #1. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim1_Grid.mat,  gauss  distribution) 


Figure  E.2:  ICSS  for  Gauss  distributions  using  DYSE  map  #1. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim1_Grid.mat,  multi_clus  distribution) 


Figure  E.3:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #1. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim2_Grid.mat,  uniform  distribution) 


Figure  E.4:  ICSS  for  uniform  distributions  using  DYSE  map  #2. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim2_Grid  mat,  gauss  distribution) 


Figure  E.5:  ICSS  for  Gauss  distributions  using  DYSE  map  #2. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim2_Grid  mat,  multi_clus  distribution) 


Figure  E.6:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #2. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim3_Grid.mat,  uniform  distribution) 


Figure  E.7:  ICSS  for  uniform  distributions  using  DYSE  map  #3. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim3_Grid.mat,  gauss  distribution) 


Figure  E.8:  ICSS  for  Gauss  distributions  using  DYSE  map  #3. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim3_Grid.mat,  multi_clus  distribution) 


Figure  E.9:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #3. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim4_Grid.mat,  uniform  distribution) 


Figure  E.10:  ICSS  for  uniform  distributions  using  DYSE  map  #4. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim4_Grid.mat,  gauss  distribution) 


Figure  E.ll:  ICSS  for  Gauss  distributions  using  DYSE  map  #4. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim4_Grid  mat,  multi_clus  distribution) 


Figure  E.12:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #4. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim5_Grid  mat,  uniform  distribution) 


Figure  E.13:  ICSS  for  uniform  distributions  using  DYSE  map  #5. 


0  95 


0.9 


co 

to 

o 

0) 

cn 

2 

<u 

> 

< 


085 


0.8 


0.75 


0.7 


0.65 


40  60 

Number  of  Clusters 


100 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim5_Grid.mat,  gauss  distribution) 


Figure  E.14:  ICSS  for  Gauss  distributions  using  DYSE  map  #5. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim5_Grid.mat,  multi_clus  distribution) 


Figure  E.15:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #5. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim6_Grid.mat,  uniform  distribution) 


Figure  E.16:  ICSS  for  uniform  distributions  using  DYSE  map  #6. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim6_Grid.mat,  gauss  distribution) 


Figure  E.17:  ICSS  for  Gauss  distributions  using  DYSE  map  #6. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim6_Grid  mat,  multi_clus  distribution) 


Figure  E.18:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #6. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim7_Grid.mat,  uniform  distribution) 


Figure  E.19:  ICSS  for  uniform  distributions  using  DYSE  map  #7. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim7_Grid.mat,  gauss  distribution) 


Figure  E.20:  ICSS  for  Gauss  distributions  using  DYSE  map  #7. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim7_Grid.mat,  multi_clus  distribution) 


Figure  E.21:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #7. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim8_Grid.mat,  uniform  distribution) 


Figure  E.22:  ICSS  for  uniform  distributions  using  DYSE  map  #8. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim8_Grid.mat,  gauss  distribution) 


Figure  E.23:  ICSS  for  Gauss  distributions  using  DYSE  map  #8. 


ICSS  vs.  Cluster  Count  for  100  Nodes 


(using  Sim8_Grid  mat,  multi_clus  distribution) 


Figure  E.24:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #8. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim9_Grid.mat,  uniform  distribution) 


Figure  E.25:  ICSS  for  uniform  distributions  using  DYSE  map  #9. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim9_Grid.mat,  gauss  distribution) 


Figure  E.26:  ICSS  for  Gauss  distributions  using  DYSE  map  #9. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim9_Grid.mat,  multi_clus  distribution) 


Figure  E.27:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #9. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim10_Grid.mat,  uniform  distribution) 


Figure  E.28:  ICSS  for  uniform  distributions  using  DYSE  map  #10. 
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ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim10_Grid  mat,  gauss  distribution) 


Figure  E.29:  ICSS  for  Gauss  distributions  using  DYSE  map  #10. 


ICSS  vs.  Cluster  Count  for  100  Nodes 
(using  Sim10_Grid.mat,  multi_clus  distribution) 


Figure  E.30:  ICSS  for  multi-cluster  distributions  using  DYSE  map  #10. 
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Appendix  F:  MATLAB  Code 


MATLAB  code  was  used  to  test  the  fc-means  clustering  algorithm  and  its  heuristics 
against  various  parameter  configurations.  In  addition  to  the  following  files,  parseArgs  .m 
and  subaxis. m  (both  available  online)  are  also  used.  To  load  the  spectrum,  run 
load_spectrum.m.  To  run  the  test  itself,  run  kmeans_experiment  .m.  For  reference, 
the  total  MATLAB  code  base  developed  for  this  experiment  is  863  lines,  including  vertical 
whitespace. 
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function  spectrum  =  build_baseline_spect ra  (type,  dim,  bins) 
spectrum  =  zeros (dim,  dim,  bins); 

if  strcmp(type,  'empty’)  ==  1 
for  i  =  1 : dim 

for  j  =  1 : dim 

spectrum(i,  j,  :)  =  ones (bins,  1); 

end 

end 

elseif  strcmp(type,  'disjoint')  ==  1 
interval  =  bins  /  (dim*dim) 
counter  =  1; 

if  interval  <  1 

fprintf (' Cannot  build  disjoint  baseline — not  enough  ... 

bins . \n\n ' ) ; 
spectrum  =  []; 
return; 

end 

for  i  =  1 : dim 

for  j  =  1 : dim 

temp  =  zeros (bins,  1); 
temp (counter*interval )  =  1; 
spectrum (i,  j,  :)  =  temp; 
counter  =  counter  +  1; 

end 

end 

elseif  strcmp(type,  'full')  ==  1 
for  i  =  1 : dim 

for  j  =  1 : dim 

spectrum(i,  j,  :)  =  zeros  (bins,  1); 

end 

end 

end 
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35  end 


i  function  clean_data() 

2  dos ( ' del  /Q 

. \commands\commands* ' ) ; 

3  dos ( ' del  /Q 

. \output\output* ' ) ; 

4  dos ( ' del  /Q 

. \input\input* ' ) ; 

5  end 
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function  spectrum  =  decode_spectrum  (mag,  n_devices,  info_grid,  ... 
bins,  dim) 

spectrum  =  zeros (dim,  dim,  bins); 

for  x  =  1 : dim 

for  y  =  1 : dim 

for  d  =  l:n_devices 

if  info_grid(x,  y)  *  0 

if  mag(d,  info_grid(x,  y)  ,  :)  *  zeros  (1,  1,  bins) 

spectrum (x,  y,  :)  =  mag(d,  info_grid(x,  ... 

y)  r  ■)  ; 

break; 

end 

end 

end 

end 

end 

end 


1  function  intra_cluster_similarity  =  draw_cluster  s  (map_number ,  ... 

distro,  v_n_centers,  map_assigned,  n.nodes,  rems,  baseline,  ... 
plot-dusters ) 

2  heuristics  =  char (' Lloyd ' ,  'Swap',  'Hybrid',  ' EZ— Hybrid ') ; 

3  n.heuristics  =  4; 

4  plot_number  =  1; 

5  intra_cluster_similarity  =  zeros  (numel  (v_n_centers )  ,  ... 

n_heuristics ) ; 

6 

7  distros  =  {'Uniform'  'Gauss'  ' Multi-Cluster '};% ;  'co_gauss'; 

8  %  '  co_laplace  '  ;  ' clus_gauss ' ;  ' clus_orth_f lats ' ; 

9  % ' clus_ellipsoids ' ;  ' multi.clus ' ] ; 

10 

11  for  j  =  1  :  n_heuristics 

12  fprintf ( ' \n\tProcessing  heuristic:  %s\n',  heuristics ( j ,:)) ; 

13  fprintf  (  ' \t  Cluster  counts  drawn:  '); 

14 

15  for  i  =  1 :  numel  (v_n_centers  ) 

16  cluster_similarities  =  zeros  (v_n_centers  (i)  ,  2); 

n  cluster_similarities  (  :  ,  2)  =  ones  (v_n_centers  (i)  ,  1); 

18 
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33 

34 

35 

36 

37 
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if  plot_clusters  == 

set (gca,  'Units',  'normal'); 

subaxis  (n_heurist  ics ,  numel  (v_n_centers)  ,  ... 

plot_number ,  . .  . 

'Spacing',  0.01,  'MT',  0.05,  'MR',  0.02,  .. 

'MB',  0.05,  'ML',  0.02); 

if  j  ==  1 

title  (sprintf (' %d  centers',  v_n_centers  (i) ) ) ; 

end 

if  i  ==  1 

y label (heuristics ( j , : ) ) ; 

end 

hold  on 

%center  s_x_coords  =  centers  (i,  j,  ... 

1  :v_n_centers  (i)  ,  1); 

%center  s_y_coords  =  centers  (i,  j,  ... 

1 :  v_n_centers  ( i ) ,  2); 

%scatter ( centers_x_coords  ( :  ) ,  ... 

center s_y_coords  (  : ) ,  '.r'); 

nodes_x-Coords  =  map_assigned (i,  j,  1,  :); 

nodes_y_coords  =  map_assigned  (i,  j,  2,  :); 

scatter  (nodes_x_coords  (:)  ,  nodes_y_coords  ( : )  ,  '  .k' 

end 

for  k  =  l:n_nodes 

current-center  =  map_assigned (i,  j,  3,  k) ; 
current.rem  =  rems (k,  : ) ; 

cluster-similarities (current_center+l,  1)  =  ... 

cluster -similar it ies ( current_center+l , 
1)  +  1; 

for  m  =  k+l:n_nodes 

compare.center  =  map.assigned (i,  j,  3,  m)  ; 
compare.rem  =  rems (m,  : ) ; 

if  compare.center  ==  current-center 

similarity  =  1 . 0— pdist ( [ current.rem;  ... 
compare.rem] ,  'Hamming'); 

cluster-similarities  ( currentxenter+1 ,  2) 


58  cluster-similarities  ( current_center+l , 

2)  +  similarity; 

59 

60  if  plot-duster s  ==  1 

61  xl  =  map_assigned  (i,  j,  1,  k)  ; 
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62 

yi 

=  map_assigned (i. 

ir 

2, 

k>  ; 

63 

x2 

=  map_assigned (i. 

ir 

i. 

m)  ; 

64 

y2 

=  map_assigned (i. 

ir 

2, 

m)  ; 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 


colors  =  colormapO; 

plot ( [xl  x2],  [yl  y2],  'Color',  ... 

colors (ceil (( 63*similarity) )  +1,  :)); 


end 


end 


end 


end 


temp_similar ity  =  cluster_similarities ( : ,  1); 

for  k  =  1 :  v_n_centers  (i) 

if  cluster_similarities (k,  1)  4  0 

cluster_similarities (k,  1)  =  ... 

ceil ( cluster_similarities (k,  1)~2  /  2); 

temp_similarity (k)  =  cluster_similarities (k, 

2)  /  cluster_similarities (k,  1); 

end 

end 

%dot (cluster_similarities  (  :  ,  1),  ... 

cluster_similarities  (  :  ,  2))  /  ... 

sum (cluster.similarities ( : ,  1)); 
%intra_cluster_similarity  (i,  j)  =  ... 
dot ( cluster_similarities  (  : ,  1),  ... 

1 . /cluster_similarities  (  :  ,  2)); 
intra_cluster_similarity  (i,  j)  =  mean (temp_similarity ) ; 

if  plot-duster s  ==  1 
hold  off 

set (gcf , 'Color', 'white') 

set (gcf ,  'Position',  [1  1  1920  978]) 

axis ( [0  10  0  10] ) ; 

set  (gca,  'xtick',  []); 

set (gca,  'ytick',  []); 

%set (gcf ,  'Units',  'normal'); 

%set (gca,  'Position',  get (gca,  ' OuterPosition ' ) ) ; 


plot_number  =  plot_number  +  1; 


end 


fprintf('%d  ',  v_n_centers  (i) ) ; 


end 

end 

if  baseline  ==  0 
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105 


export_f ig (gcf ,  ... 

sprint f ( '  . /figures/ Clustering/ clustering%d_%s . png '  ,  ... 

mapnumber,  char  (distros  (distro)  )  )  )  ; 

106  else 

107  export_fig  (gcf ,  ... 

sprint f ( ' ./ figures/ Clustering/baseline%d . png ' ,  baseline) ) ; 

108  end 

109 

no  close (gcf) ; 

in 

in  fprintf (  ' \n ' ) ; 

in  end 

114 

ns  function  e  =  n_edges  (n_vertices ) 

ii6  e  =  (n_vertices"2  /  2)  -  (n_vertices  /  2); 

in  end 


1  function  draw-distribution  (map_number,  n_nodes,  centers,  distro,  ... 

center.count,  dim,  baseline) 

2  distros  =  {'Uniform'  'Gauss'  'Multi-Cluster' }; %;  'co.gauss'; 

3  % ' co.laplace  '  ;  ' clus.gauss ' ;  ' clus_orth_f lats ' ; 

4  % ' clus_ellipsoids ' ;  ' multi_clus ' ] ; 

5 

6  seed  =  1; 

7 


8 

for  m  = 

1 : map_number 

9 

if  i 

m  ==  3 

10 

seed  =  seed  +  1; 

11 

12 

end 

13 

for 

d  =  l:distro 

14 

fig  =  figure (1000  +  (m-l)*distro  +  d)  ; 

15 

set(fig,  'Position',  [100  100  500  500]) 

16 

gen_kmeans_commands  (m,  nmodes,  centers,  d,  ... 
center.count,  seed,  baseline); 

17 

run_kmeans (m,  center.count) ; 

18 

points  =  read_kmeans_points (m,  n.nodes,  dim,  baseline) 

19 

20 

hold  on 

21 

h  =  scatter (points (: ,  1),  points (:,  2),  '  . k ' )  ; 

22 

hChildren  =  get (h,  'Children'); 

23 

24 

set (hChildren,  ' Markersize '  ,  15) 

25 

o, 

o 

title (sprintf (' %s  Distribution\nseed  =  %d',  ... 

char  (distros  (distro)  )  ,  seedmumber)  ) 

26 

O, 

O 

xlabel('X  Location') 

27 

28 

O, 

O 

ylabel('Y  Location') 

29 

set (gcf , 'Units', ' normal ' ) 

30 

set (gcf, 'Color', 'white') 

31 

axis  (  [0  10  0  10] ) 
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set (gca, ' Position ',[ 0 . 05  0.05  0.91  0.91]) 
set (gca, 'Box', 'on') 


%  x— axis  parameters 

set  (gca,  ' XTickMode ' ,  'manual') 

set (gca, ' XMinorTick ' ,  'on') 

set (gca, ' XTick ' , 0 : 10) 

set (gca, ' XAxis Location ' , ' Top ' ) 


42 

%  y— axis 

parameters 

43 

set (gca, ' 

YTickMode ' ,  'manual' 

44 

set (gca, ' 

YMinorTick ' ,  ' 

on '  ) 

45 

set (gca, 1 2 3 4 5 6 7 8 9 10 * 12 13 

YTick ' ,0:10) 

46 

set (gca, ' 

YDir ' , ' reverse 

' ) ; 

47 

48 

grid  on 

49 

set (gca, ' 

GridLineStyle ' 

f  _  1 

r 

1 ) ; 

50 

%set (gca. 

' XColor ' , [0.5 

LO 

o 

0  . 

51 

%set (gca. 

' YColor ' , [0.5 

LO 

o 

0  . 

52 

53 

export_f i 

■  g (gcf,  . . . 

sprint f  (  '  . \\ f igures\\Node_maps\\% s _%d . png ' , char (distro4 (d) ) , 
m)  )  ; 

close (gcf ) ; 

hold  off 
fprintf ( ' | ' ) ; 


end 

seed  =  seed  +  1; 

end 

fprintf ( ' \n ' ) ; 


end 


1  function  gen_kmeans_commands  (map_number ,  n.nodes,  centers, 

distro_indices ,  center_count ,  seed-number,  baseline) 

2  heuristics  =  {'lloyd',  'swap',  'hybrid',  ' EZ-hybrid ' } ; 

3 

4  for  h  =  1 : numel (heuristics ) 

5  show.assignments  =  ' show_assignments  yes ' ; 

6  kcenters  =  'kcenters ' ; 

7 

8  %  Point  generation  commands 

9  print_points  =  ' print_points  yes'; 

10  distribution  =  'distribution'; 

n  data_size  =  'data_size'; 

12  gen_data_pts  =  '  gen_data_pt s  '  ; 

13  read.data.pts  =  '  read_data_pts  '  ; 
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14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


seed  =  ' seed 1 ; 

%  Point  distribution  types 

distros  =  {'uniform'  'gauss'  ' mult i_clus ' } ; %;  'co_gauss  '  ; 

% ' co.laplace  '  ;  ' clus.gauss ' ;  ' clus_orth_f lat s ' 

% ' clus .ellipsoids  '  ;  ' multi_clus ' ] ; 

%  Point  distribution  type  parameters 


stcLdev  =  'stcLdev'; 

o, 

o 

for 

gauss,  clustered;  . . 

def  =  1 

std_dev_lo  =  '  std_dev_lo  '  ; 

-  1 

o, 

o 

for 

clus.ellipsoids ;  def 

—  ± 

std_dev_hi  =  '  std_dev_hi  '  ; 

=  i 

o, 

o 

for 

clus.ellipsoids ;  def 

corrxoef  =  1  corr_coef'; 

o, 

o 

for 

co_gauss ,  ... 

co.laplace;  def  =  0.05 

colors  =  'colors'; 

o, 

o 

for 

clustered;  def  =  5 

max_clus_dim  =  ' max_clus_dim ' 

•  9- 

r  o 

for 

clus_orth_f lats ,  ... 

clus -ellipsoids; 

o, 

o 

def  =  1 

%  Runtime  command/parameters 
run.kmeans  =  '  run.kmeans  '  ; 

%  System  commands 
quit  =  ' quit ' ; 

%  Open  the  file  for  writing 

fid  =  fopen  (sprintf ( '  . /commands/commands%d_%d_%s ' ,  ... 

map_number,  center_count,  char (heuristics (h) )) ,  'w'); 

%  Write  commands  to  the  command  input  file 
fprintf (fid,  '%s\n',  show_assignment s ) ; 
fprintf (fid,  '%s\n',  print_points) ; 
fprintf (fid,  '%s  %d\n',  seed,  seed-number) ; 
fprintf (fid,  '%s  %d\n',  data_size,  n_nodes); 

for  i  =  dist ro.indices 
%  Set  gauss  std.dev 
if  i  ==  2  ||  i  ==  3 

fprintf (fid,  '%s  %f\n',  std_dev,  0.3); 

end 

if  i  ==  3 

fprintf (fid,  '%s  %d\n',  colors,  2); 

end 

fprintf (fid,  '%s\n',  sprintf ('%s  %s',  distribution, 
char (distros (i) ) ) ) ; 

if  baseline  ==  0 

fprintf (fid,  '%s\n',  gen_dat a_pt s )  ; 
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60 

61 

62 

63 

64 
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66 

67 

68 

69 

70 

71 

72 
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2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
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22 

23 

24 
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26 

27 

28 

29 

30 

31 

32 


else 


fprintf (fid,  !%s  . /input /baseline%d * ,  ... 

read_data_pt  s ,  baseline); 

end 

fprintf (fid,  ' \n'); 

fprintf (fid,  '%s\n',  sprintf('%s  %df,  kcenters, 
center_count) ) ; 

fprintf (fid,  ' %s\n',  sprintf('%s  %s',  run.kmeans, 
char (heuristics  (h) ) ) ) ; 

end 


fprintf (fid,  ,%s\n',  quit); 


f close (fid) ; 

end 

end 


function  rems  =  gen.rems (spectrum,  n_nodes,  dim,  thresh,  bins) 
rems  =  ones (n.nodes,  bins); 

%aggregate  =  ones  (1,  bins); 
rem  =  zeros (1,  bins); 

%fid  =  f open ( filename,  'w'); 
rem.count  =  1; 

%if  fid  ==  -1 

%  error (' could  not  open  file1); 

%end 


for  x  =  1 : dim 

for  y  =  1 : dim 

cell  =  spectrum (x,  y,  :); 

c  =  cell ( : ) ; 

thresh  =  mean (cell); 

if  c  ==  zeros (bins,  1) 
rem  =  ones (1,  bins); 

else 

for  i  =  1 : bins 

if  c ( i )  <  thresh 
rem ( i )  =  1 ; 

else 

rem ( i )  =  0 ; 

end 

end 

end 

rems ( rem_count ,  :)  =  rem; 

%  aggregate  =  aggregate  &  rem; 

%  fprintf (fid,  '%s\n',  dec2hex (rem) ' ) ; 
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46 

47 
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rem_count  =  rem.count  +  1; 

end 

end 

%fprintf ( f id,  '\n%s\n',  dec2hex (aggregate) ') ; 

%openChannels  =  0; 

%f or  i  =  1 : bins 
%  if  aggregate (i)  ==  1 

%  openChannels  =  openChannels  +  1; 

%  end 

%end 

%fprintf ( ' Open  Channels:  %d\n',  openChannels); 
%f close (fid) ; 

end 


1  function  maps  =  gen_test_maps (n_maps ,  size,  dim) 

2  maps  =  zeros  (njnaps,  3,  size); 

3 

4  for  i  =  l:n_maps 

5  maps (i,  1:2,  :)  =  gen_data_pts (size,  dim, 

sprintf ( ' . /input/input%d' ,  i) ) ' ; 

6  maps  (i,  3,  :  )  =0; 

7  end 

8  end 


1  function  kmeans_experiment  (n.trials,  n.nodes,  v_n_centers,  ... 

thresh_ad just ,  map-dimensions,  bins,  spectrummumber,  ... 
distro_indices ,  baseline,  plot_clusters ) 

2  lloyd  =  zeros (numel  (v_n_centers ) ,  numel (distro_indices )) ; 

3  swap  =  zeros (numel  (v_n_centers ) ,  numel (distro_indices )) ; 

4  hybrid  =  zeros (numel (v_n_centers ) ,  numel (distro_indices )) ; 

5  ez_hybrid  =  zeros (numel (v_n_centers ) ,  numel (distro_indices) ) ; 

6 

7  n_heuristics  =  4; 

8  n_centers  =  numel (v_n_centers )  ; 

9 

10  for  t  =  0.6:0. 1:1. 4 

n  thresh_ad just  =  t; 

12 

13  %  Initialize  spectrum  data  if  it  does  not  already  exist 

14  fprintf (' Generating  spectrum  data  using  Sim%d_Grid . mat  ...',  ... 

spectrum.number) ; 

15  spectrum  =  load_spectrum  (  spectrum.number )  ; 

16  fprintf  (  '  Done  .  \n  '  )  ; 

17 
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18 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30  % 

31 

32  % 

33 

34 

35 

36 

37 

38 

39 

40  % 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 


%  Point  distribution  types 

distros  =  {'uniform1  'gauss'  ' mult i_clus ' } ; %;  'co_gauss  '  ; 

% ' co_laplace  '  ;  ' clus_gauss ' ;  '  clus_orth_f  lat  s  '  ; 

% ' clus_ellipsoids  '  ;  ' mult i.clus ' ] ; 


tic 

rems  =  gen.rems  (spectrum,  n_nodes,  map-dimensions ,  ... 

thresh_ad just ,  bins)  ; 

fclose  all; 
clean.data  ( ) ; 

icss_f  iguremumber  =  numel  (dist ro_indices )  *n_t rials ; 

for  j  =  1 : numel (dist ro_indices ) 
for  j  =  distro_indices 
seed_number  =  0; 

int  ra_cluster_similarity  =  zeros  (n_trials,  ... 

numel  (v_n_centers)  ,  Of¬ 
fprint  f  (' Running  experiment  using  %s  distribution  ... \n ' , 
char (distros ( j ) ) )  ; 

for  i  =  l:n_trials 
for  i  =  n.trials 

seed.number  =  seed_number  +  1; 

%  When  seed  is  3,  kmltest.exe  crashes.  (??) 
if  seed_number  ==  3 
seed-number  =  4; 

end 

fprintf ( ' Generating  K— means  commands ...'); 
for  v  =  v_n_centers 

gen_kmeans_commands  (if  n.nodes,  v_n_centers,  j, 
Vf  seed-number,  baseline)  ; 

end 

fprintf  (  '  Done  .  \n  Of¬ 
fprint  f  (' Running  kmltest.exe  for  map  %d...',  i) ; 
for  v  =  v_n_centers 
run_kmeans  (i,  v) ; 

end 

fprintf  (  '  Done  .  \n  Of¬ 
fprint  f  (' Parsing  points  for  this  distribution...'); 
points  =  read_kmeans_points  (if  n_nodes,  ... 

map-dimensions ,  baseline); 
fprintf ( ' Done . \n ') ; 
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fprintf ( 1  Parsing  assignments  for  map  %d...',  i) ; 
assignments  =  read_kmeans_assignments  (i,  n.nodes,  ... 

numel  ( v_n_cent er  s  )  )  ; 
fprintf  (  ' Done . \n ' ) ; 

fprintf (' Parsing  centers  for  map  %d...’,  i) ; 
centers  =  read_kmeans_centers  (i,  n_nodes,  ... 

numel  ( v_n_cent er  s  )  )  ; 
fprintf  (  ' Done . \n ' ) ; 

fprintf (' Entering  assignments  for  map  %d...',  i) ; 
map_assigned  =  set_assignments  (n_nodes,  points,  ... 

assignments,  v_n_center s ) ; 
fprintf ( 1  Done . \n ' ) ; 

%  Could  use  FFT  algorithm  for  faster  similarity  . . . 
comparison 

fprintf (' Drawing  clusters  and  determining  ... 

intra— cluster  similarity. . . ' ) ; 

%figure (( j-1 ) *n_t  rials  +  i) 

%colorbar 

%colorbar ( ' YTickLabel ' ,  {'0%',  '25%',  '50%',  '75%',  .. 

' 100%  '  },  ... 

%  'location',  ' EastOutside ' ) 

temp  =  draw.clusters  (i,  j,  v_n_centers,  map_assigned, 
n.nodes,  rems,  baseline,  plot_cluster s ) ; 
intra_cluster_similarity (i,  :,  :)  =  temp; 

fprintf ( ' Done . \n ' ) ; 

fprintf ('Map  %d  processing  complete ! \n ' ,  i); 

fprintf ( ' ========================================\n ' )  ; 

end 

for  k  =  1 : numel  (v_n_centers ) 

lloyd(k,  j)  =  mean ( intra_cluster_similarity  (  : ,  k,  1)); 
swap(k,  j)  =  mean ( intra_cluster_similarity ( : ,  k,  2)); 
hybrid(k,  j)  =  mean (intra_cluster_similarity  (  : ,  k,  3)) 
ez_hybrid(k,  j)  =  mean  (  intra_cluster_similarity ( : ,  k, 
4)  )  ; 

end 


fprintf (' Plotting  intra— cluster  similarity  for  %s  .. 
distribution . . . ' ,  char (distros ( j ) ) ) ; 

%int ra_clus ter .similarity 

figure (n_trials*numel (distro_indices)  +  j) 
hold  on 

plot (log2  (v_n_centers) ,  intra_cluster_similarity  (  :  , 

:,  1),  '—  r+') 

plot (log2  ( v_n_center s ) ,  intra_cluster_similarity  (  :  , 

:,  2),  ’ — b+ ' ) 
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108  % 
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111  % 

112  % 
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123 
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127 

128 
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131 
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133 

134 
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plot (log2  ( v_n_center s ) ,  intra_cluster_similarity  (  :  , 

' — k+ ' ) 

plot (log2  ( v_n_center  s ) ,  intra_cluster_similarity  (  :  , 
g+') 

set (gca,  ' XTickLabel ' ,  v_n_centers) 
xlabel (' Number  of  Clusters') 
ylabel (' Average  ICSS') 

title  ( sprintf (' ICSS  vs.  Cluster  Count  for  %s  ... 
Distribution '  ,  upper (char (distros ( j ) ) ) ) )  ; 
hold  off 

fprintf ( ' Done . \n ' ) ; 

fprintf (' Experiment  complete  for  %s  ... 
distribution ! \n ' ,  char (distros ( j ) ) ) ; 

fprintf ( ' ========================================\n\n 

colors  =  ['r'  'b'  'g'  ' k ' ] ; 

Plot  ICSS  for  Lloyd 

fprintf (' Plotting  ICSS  for  all  heuristics...'); 
figure  (n_heuristics*n_centers  +  1) 
h  =  zeros ( 4 ,  1 ) ; 

labels  =  {'Lloyd',  'Swap',  'Hybrid',  ' EZ-Hybrid ' } ; 
hold  on 

n_icss  =  n_trials*numel  (v_n_centers  )  ; 
cluster  =  zeros (n.icss,  1); 
icss  =  zeros (n_icss ,  1); 

for  m  =  1 : numel  (v_n_centers ) 
for  n  =  l:n_trials 

cluster  (  (m-1 ) *n_trials  +  n)  =  v_n_centers(m); 
icss ( (m-1) *n_trials  +  n)  =  ... 

intra_cluster_similarity  (n,  m,  1); 

end 

end 

boxplot (icss,  log (cluster ) ) 

h(l)  =  plot ( v_n_centers ,  lloyd(:,  j),  ' — ro ' ) ; 

h(2)  =  plot ( v_n_center s ,  swap ( : ,  j),  '— bs'); 

h(3)  =  plot (v_n_centers ,  hybrid (: ,  j),  '— g*'); 

h(4)  =  plot  (v_n_centers ,  ez_hybrid(:,  j),  '— cv'); 

set(gcf, 'Color', 'white') 

set (gca,  'XTickMode',  'auto') 

set (gca,  ' XTickLabelMode ' ,  'auto') 

set (gca, ' Box ' , ' on ' ) 

set (gca, ' XMinorTick ' ,  'on') 

set (gca, ' YMinorTick ' ,  'on') 

axis([0  n_nodes  0.6  1]); 


,  3), 
,  4), 
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xlabel (' Number  of  Clusters') 
ylabel (' Average  ICSS ' ) 

title  (sprintf (' ICSS  vs.  Cluster  Count  for  %d  ... 

Nodes\n (using  Sim%d_Grid . mat ,  %s  distribution)',  .. 
n.nodes,  spectrum_number ,  char (distros ( j ) ) ) ,  ... 

'Interpreter',  'none'); 
legend(h,  labels,  'Location',  'Northwest'); 
hold  off 

fprintf ( ' Done . \n ' ) ; 
export_f ig (gef ,  ... 

sprintf ( ' . /f igures/ICSS/icss_sim%d_%s .png ' ,  ... 

spectrum_number ,  char (distros ( j ) ) ) ) ; 
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o_ 

o 

close (gef ) ; 

162 

163 

o. 

o 

icss_f igure.number  =  icss_f igure_number  +  1; 

164 

o_ 

o 

165 

o_ 

o 

Plot  ICSS  for  Swap 

166 

167 

o_ 

o 

fprintf ( 'Plotting  ICSS  for  Swap. . . ' ) ; 

168 

o_ 

o 

figure (icss_f igure.number) 

169 

o. 

o 

h_swap  =  zeros (numel (thresholds) , 

1)  ; 

170 

o_ 

o 

labels  =  cell (numel (thresholds )  , 

i)  ; 

171 

o_ 

o 

172 

o. 

o 

hold  on 

173 

o. 

o 

n_icss  =  n_trials*numel  (v_n_centers)  ; 

174 

o_ 

o 

cluster  =  zeros (n_icss ,  1); 

175 

o_ 

o 

icss  =  zeros  (n.icss,  1); 

176 

o_ 

o 

177 

o, 

o 

for  m  =  1 : numel (v_n_centers ) 

178 

o_ 

o 

for  n  =  l:n_trials 

179 

o_ 

o 

cluster ( (m—1) *n_trials  +  n)  = 

v_n_cent< 

180 

o_ 

o 

icss ( (m-1) *n_trials  +  n)  =  .. 

intra_cluster_similarity (n,  m,  2); 

181 

o. 

o 

end 

182 

o_ 

o 

end 

183 

184 

o_ 

o 

boxplot (icss,  log (cluster ) ) 

185 

186 

o_ 

o 

for  r  =  1 : numel (thresholds ) 

187 

o_ 

o 

h_swap  (r)  =  plot  ( v_n_center s  , 

swap (r. 

' -bo ' ) ; 

188 

o. 

o 

labels (r)  =  {char (sprintf (' Threshold  @ 

thresholds (r) ) ) }; 

189 

o_ 

o 

end 

190 

191 

o. 

o 

axis  ( [ 0 . 90*log2 (min ( v_n_centers ) ) 

1 . 05*log2 (max (v_n_centers ) )  ... 

192 

o_ 

o 

0 . 95*min (min (min (lloyd) ,  min (swap) 

)  ,  ... 

min (min (hybrid) ,  ... 

193 

o_ 

o 

min  (ez_hybrid)  ) )  ... 
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194 

o_ 

o 

1 . 05*max (max (max (lloyd) ,  max (swap)),  .. 

max (max (hybrid) ,  ... 

195 

o_ 

o 

max  (ez_hybrid) ) ) ] ) 

196 

197 

o_ 

o 

set (gca,  'XTickMode',  'auto') 

198 

o_ 

o 

set (gca,  ' XTickLabel ' ,  v_n_centers) 

199 

o_ 

o 

200 

o_ 

o 

xlabel (' Number  of  Clusters') 

201 

o_ 

o 

ylabel (' Average  ICSS') 

202 

o_ 

o 

title  ( sprintf (' ICSS  vs.  Cluster  Count 

for  Swap 

Heuristic\n  (%s  distribution)',  char (distros ( j ) ) ) 

) ; 

203 

o. 

o 

legend (inswap ,  labels,  'Location',  ' SouthEast ' ) 

204 

o. 

o 

hold  off 

205 

o_ 

o 

fprintf ( ' Done . \n ' ) ; 

206 

207 

o. 

o 

icss_f igure.number  =  icss_f igure.number  + 

i; 

208 

o_ 

o 

209 

o_ 

o 

Plot  ICSS  for  Hybrid 

210 

211 

o_ 

o 

fprintf (' Plotting  ICSS  for  Hybrid...') 

; 

212 

o_ 

o 

figure  (icss_f iguremumber) 

213 

o. 

o 

h_hybrid  =  zeros (numel (thresholds )  ,  1) 

; 

214 

o. 

o 

labels  =  cell (numel (thresholds ) ,  1); 

215 

216 

o. 

o 

hold  on 

217 

o. 

o 

n_icss  =  n_trials*numel  ( v_n_centers  )  ; 

218 

o_ 

o 

cluster  =  zeros (n.icss,  1); 

219 

o_ 

o 

icss  =  zeros (n.icss,  1); 

220 

o_ 

o 

221 

o_ 

o 

for  m  =  1 : numel (v_n_centers ) 

222 

o_ 

o 

for  n  =  l:n_trials 

223 

o_ 

o 

cluster ( (m-1) *n_trials  +  n)  = 

v_n_cent  e: 

224 

o_ 

o 

icss ( (m— 1) *n_trials  +  n)  =  ... 

intra_cluster_similarity  (n,  m,  3); 

225 

o. 

o 

end 

226 

o. 

o 

end 

227 

228 

o_ 

o 

boxplot (icss,  log (cluster ) ) 

229 

o_ 

o 

230 

o_ 

o 

h_hybrid  =  plot  ( v_n_centers  ,  hybrid(r. 

j), 

231 

o_ 

o 

labels (r)  =  {char ( sprintf (' Threshold  @ 

%d  dB',  . 

thresholds (r) ) ) }; 

232 

o. 

o 

233 

o. 

o 

axis  ( [ 0 . 90*log2 (min  ( v_n_centers )  )  ... 

1 . 05*log2 (max ( v_n_center s ) )  ... 

234 

o_ 

o 

0 . 95*min (min (min (lloyd) ,  min (swap) ) ,  . . 

.  . 

min (min (hybrid) ,  ... 

235 

o_ 

o 

min  (ez_hybrid) ) )  ... 

236 

o_ 

o 

1 . 05*max (max (max (lloyd) ,  max (swap)),  ., 

.  . 

max (max (hybrid) ,  ... 

237 

o_ 

o 

max (ez_hybrid)  ) ) ] ) 

238 
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239 

o_ 

o 

set (gca,  'XTickMode',  'auto') 

240 

o_ 

o 

set (gca,  ' XTickLabel ' ,  v_n_centers) 

241 

o. 

o 

242 

o. 

o 

xlabel (' Number  of  Clusters') 

243 

o_ 

o 

ylabel (' Average  ICSS') 

244 

o_ 

o 

title  ( sprintf (' ICSS  vs.  Cluster  Count 

for  Hybrid  . 

Heuristic\n (%s  distribution)',  char (distros ( j ) ) ) 

) ; 

245 

o_ 

o 

legend (h_hybrid,  labels,  'Location', 

' SouthEast ' ) ; 

246 

o_ 

o 

hold  off 

247 

o_ 

o 

fprintf ( ' Done . \n ' ) ; 

248 

249 

o. 

o 

icss_f igure.number  =  icss_f igure_number  + 

i; 

250 

o. 

o 

251 

o_ 

o 

Plot  ICSS  for  EZ-Hybrid 

252 

253 

o. 

o 

fprintf (' Plotting  ICSS  for  EZ— Hybrid. 

254 

o_ 

o 

figure (icss_f igure.number) 

255 

o_ 

o 

h_ez_hybrid  =  zeros (numel (thresholds ) 

,  1); 

256 

o_ 

o 

labels  =  cell (numel (thresholds ) ,  1) ; 

257 

o_ 

o 

258 

o, 

o 

hold  on 

259 

o. 

o 

n_icss  =  n_trials*numel  (v_n_centers)  ; 

260 

o, 

o 

cluster  =  zeros (n.icss,  1); 

261 

o_ 

o 

icss  =  zeros (n.icss,  1); 

262 

o. 

o 

263 

o_ 

o 

for  m  =  1 :  numel  (v_n_centers  ) 

264 

o_ 

o 

for  n  =  l:n_trials 

265 

o_ 

o 

cluster  (  (m-1)  *n_trials  +  n)  = 

v_n_centers  (] 

266 

o_ 

o 

icss ( (m-1) *n_trials  +  n)  =  .. 

. 

intra_cluster_similarity  (n,  m,  4); 

267 

o_ 

o 

end 

268 

o_ 

o 

end 

269 

o_ 

o 

270 

o. 

o 

%boxplot (icss,  log  (cluster) ) 

271 

o. 

o 

272 

o. 

o 

for  r  =  1 : numel (thresholds ) 

273 

o_ 

o 

h_ez_hybrid  (r)  =  plot  (v_n_centers 

,  ez_hybrid (r 

: ,  j ) ,  ' -ko ' ) ; 

274 

o_ 

o 

labels (r)  =  {char ( sprintf (' Threshold  @  %d  dB ' , 

thresholds (r) ) ) }; 

275 

o_ 

o 

end 

276 

o, 

o 

277 

o. 

o 

set (gca,  'XTickMode',  'auto') 

278 

o. 

o 

set (gca,  'XTickLabel',  v_n_centers) 

279 

o. 

o 

280 

o_ 

o 

xlabel (' Number  of  Clusters') 

281 

o. 

o 

ylabel (' Average  ICSS') 

282 

o. 

o 

title  ( sprintf (' ICSS  vs.  Cluster  Count 

for  EZ— Hybri' 

Heuristic\n (%s  distribution)',  char (distros ( j )) ) 

)  ; 

283 

o_ 

o 

legend (h_ez_hybrid,  labels,  'Location 

'  ,  '  SouthEast 

284 

o_ 

o 

hold  off 

285 

o, 

o 

fprintf ( ' Done . \n ' ) ; 
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286 

287  %  icss_f  igure_number  =  icss_f  igure_number  + 

288  end 

289  toe 

290 

291  %  Catch  any  file  handles  not  closed  due  to  an  error 

292  f close  all ; 

293 

294  end 

295  end 


i; 


1  function  spectrum  =  load_spectrum  (simmumber) 

2  load (sprintf (' Sim%d_Grid. mat ' ,  simmumber)  ,  ... 

3  ' device_gr id ' ,  'freq',  'info_grid',  'mag'); 

4  spectrum  =  decode.spectrum (mag,  7,  info_grid,  2048,  10); 

5  end 


1  function  plot.spectrum (map.number,  spectrum,  thresh_ad just ,  bins) 

2  figure  (map.number)  ; 

3 

4  title (sprintf (' RF  Spectrum  Maps  Generated  by  ... 

DYSE\nSim%d_Gr  id  .  mat  '  ,  mapmumber)  ) 

5 

6  plot.number  =  1; 

7  unknown  =  0; 


9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 


colors  =  ['r'  'b'  'g'  ' k ' ] ; 

for  y  =  1:10 
%for  y  =  4 

for  x  =  1:10 
%for  x  =  3 

set (gca,  'Units',  'normal'); 

subaxis(10,  10,  plot_number,  'Spacing',  0.015,  ' MT ' ,  ... 

0.03,  'MR',  0.02,  'MB',  0.035,  'ML',  0.02); 

title (sprintf (' %d,  %d',  x,  y) ) 

s  =  spectrum (x,  y,  :); 

threshold  =  thresh.ad just *mean ( s ( : ) )  ; 

plot ( s ( : ) ) ; 

axis ( [0  2047  -150  10] ) ; 
hold  on 

plot([0  bins],  [threshold  threshold],  ' — ' ,  'Color',  ... 

(1/255)*  [5,  136,  5],  'LineWidth',  2); 
plot  (  [0  bins] ,  [mean (s  (  : ) )  mean (s  (  : ) ) ] ,  '  '  ,  ... 

'Color',  [1,  0,  0],  'LineWidth',  2); 
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31 
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34 
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37 

38 

39 
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50 

51 
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55 

56 

57 

58 

59 

60 
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1 
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12 

13 
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15 


hold  off 

set  (gca,  'xtick',  []); 
set  (gca,  'ytick',  []); 


if  x  ==  1 

set (gca,  'YTickMode',  ’manual') 
set (gca,  ' YTick',  [-120,  -60,  0]) 
elseif  x  ==  10 

set (gca,  'YTickMode',  'manual') 
set (gca,  'YTick',  [-120,  -60,  0]) 
set (gca,  ' YAxisLocation ' ,  'right') 

end 

if  y  ==  1 

set (gca,  'XTickMode',  'manual') 
set (gca,  'XTick',  [0  1023  2047]) 
set (gca,  ' XAxisLocation ' ,  'top') 
elseif  y  ==  10 

set (gca,  'XTickMode',  'manual') 
set (gca,  'XTick',  [0  1023  2047]) 

end 

plot.number  =  plot_number  +  1; 

end 

end 

fprintf ( ' %d  total  cells  with  unknown  properties . \n\n ' ,  unknown); 

set (gcf ,  'Color',  'white') 

set (gcf ,  'Position',  [1  1  1920  978]) 

export_fig (gcf ,  sprintf  (  '  . \\ f igures\\RF_maps \\dyse%d . png ' ,  ... 

map_number)  )  ; 

end 


function  assignments  =  read_kmeans_assignment s  (map.number,  .  .  . 
node_count,  n.centers) 
lloyd  =  1; 
swap  =  2 ; 
hybrid  =  3; 

EZhybrid  =  4; 
type  =  0; 

n_heuristics  =  4; 

centers.count  =  — 1; 
centers_index  =  0; 

assignments  =  zeros  (n_centers  ,  n_heurist  ics ,  node_count); 
fid  =  f open ( sprintf ('. /output /output%d ' ,  map.number) ,  'r'); 
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0) 


16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


while  -ifeof(fid) 

%  Spin  until  cluster  assignments 
while  ( strcmp ( f get 1 ( f id) ,  ' [Run_k— means :' )  == 

if  feof (fid) 
break; 

end 

end 


%  Get  the  heuristic  type 
str.type  =  fscanf (fid,  '  k— means.alg 
if  strcmp ( st r_type ,  'lloyd')  *  0 
type  =  lloyd; 

elseif  strcmp ( str_type ,  'swap’)  *  0 
type  =  swap; 

elseif  strcmp ( str.type ,  'hybrid')  *  0 
type  =  hybrid; 

elseif  strcmp ( str.type ,  'EZ— hybrid')  5 
type  =  EZhybrid; 
elseif  strcmp ( st r_type , 


1)  ; 


0 


)  *  0 

%disp (' Reached  end  of  output  file.'); 

%  Reached  end  of  file  =>  no  more  results  to  read 
break; 


end 


%  Read  out  the  next  line 
fgetl (fid) ; 
fgetl (fid) ; 

%  Get  the  number  of  centers 

temp_centers  =  fscanf (fid,  '  kcenters  =  %d',  1)  ; 

if  centers-count  *  temp_centers 
centers.count  =  temp_centers ; 
centers.index  =  centers.index  +  1; 

end 


%  Spin  until  cluster  assignments 

while  (strcmp (fgetl ( fid) ,  '  (Cluster  assignments:')  ==  0) 

if  feof (fid) 
break; 

end 

end 


%  Read  out  two  following  lines 
fgetl (fid) ; 
fgetl (fid) ; 

for  i  =  l:node_count 

%  Read  the  cluster  for  the  point  number 
data  =  fscanf (fid,  ' %d  %d  %f',  3)'; 

if (isempty (data) ) 
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break; 


end 


assignments (centers_index,  type,  i)  =  data (2); 

end 

end 


f close (fid) ; 

end 


function  center.data  =  read_kmeans_center s  (map_number,  .  .  . 
node_count,  n_centers) 
lloyd  =  In¬ 
swap  =  2  ; 
hybrid  =  3; 

EZhybrid  =  4; 
type  =  0; 

n_heuristics  =  4; 

centers.count  =  — 1; 
centers.index  =  0; 

center_data  =  zeros  (n_centers,  n_heurist ics ,  node_count,  2); 
fid  =  f open ( sprintf ('. /output /output%d ' ,  map_number) ,  'r'); 

while  -ifeof(fid) 

%  Spin  until  cluster  assignments 

while  ( strcmp ( f get 1 ( fid) ,  ' [Run_k— means :' )  ==  0) 
if  feof (fid) 
break; 

end 

end 


%  Get  the  heuristic  type 

str.type  =  fscanf (fid,  1  k-means_alg  =  %s',  1); 

if  strcmp ( st r_type ,  'lloyd')  *  0 
type  =  lloyd; 

elseif  strcmp ( st r_type ,  'swap')  *  0 
type  =  swap; 

elseif  strcmp ( str.type ,  'hybrid')  *  0 
type  =  hybrid; 

elseif  strcmp ( str_type ,  'EZ— hybrid')  *  0 
type  =  EZhybrid; 

elseif  strcmp ( str.type ,  '')  *  0 

%disp (' Reached  end  of  output  file.'); 

%  Reached  end  of  file  =>  no  more  results  to  read 
break; 

end 
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%df,  1); 


%  Read  out  the  next  line 
fgetl (fid) ; 
fgetl (fid) ; 

%  Get  the  number  of  centers 
temp.centers  =  fscanf(fid,  '  kcenters 
if  centers_count  *  temp_centers 
centers.count  =  temp_centers  ; 
centers_index  =  center  s_index  +  1; 

end 

%  Spin  until  center  points 
while  ( strcmp ( fgetl ( fid) ,  '  (Final  Center  Points:')  ==  0) 

if  feof (fid) 
break; 

end 

end 

for  i  =  1 : center s_count 

temp_center  =  fscanf (fid,  '  %d  [  %f  %f  . . . 

]  dist  =  %f 1 ,  4) ' ; 

center  =  temp_center ( 1 )  +  1; 
center_x  =  temp.center  (2 )  ; 
center_y  =  temp_center  ( 3 )  ; 
center.sqdist  =  temp.center ( 4 ) ; 

center.data (centers.index,  type,  center,  1:2)  =  ... 
[center_x  center.y]  ; 

end 

end 

f close (fid) ; 

end 


function  points  =  read_kmeans_point  s  (map.number ,  node.count,  dim,  ... 
baseline) 

points  =  zeros (node.count,  3); 

fidl  =  0; 
if  baseline  ==  0 

fidl  =  fopen (sprintf ( ' . /output /output%d ' ,  map.number) ,  ’r'); 

else 

fidl  =  fopen ( sprint f ('. /input /baseline%d ' ,  baseline),  'r'); 

end 

fid2  =  fopen ( sprint f ('. /input /input%d ' ,  map_number) ,  'w'); 

i  =  0; 
x  =  0  ; 

Y  =  0; 
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16 

17 

if  baseline  ==  0 

18 

while  -ifeof  (f  idl ) 

19 

%  Spin  until  cluster  assignments 

20 

while  (strcmp (fgetl (fidl)  ,  '  (Data 

21 

if  feof (fidl ) 

22 

break; 

23 

end 

24 

end 

25 

26 

for  i  =  1 : node_count 

27 

if  feof ( fidl ) 

28 

break; 

29 

end 

30 

31 

%  Get  the  point  coordinates 

32 

temp_point  =  fscanf(fidl,  ' 

3)  '  ; 

33 

34 

x  =  (temp.point (2) +1) * (dim/2) ; 

35 

y  =  (temp.point (3) +1) * (dim/2) ; 

36 

37 

fprintf  ( f id2 ,  '%f  %f\n',  x,  y); 

38 

39 

points (temp .point  ( 1 ) +1 ,  1)  =  x; 

40 

points (temp_point ( 1 ) +1 ,  2)  =  y; 

41 

end 

42 

end 

43 

else 

44 

for 

i  =  l:node_count 

45 

if  feof (fidl ) 

46 

break; 

47 

end 

48 

49 

temp.point  =  fscanf(fidl,  ' %f  %f',  : 

50 

51 

x  =  temp.point ( 1 ) ; 

52 

y  =  temp  .point ( 2 ) ; 

53 

54 

fprintf ( fid2 ,  ' %f  %f\n',  x,  y) ; 

55 

56 

points (i,  1)  =  x; 

57 

points ( i ,  2 )  =  y ; 

58 

end 

59 

end 

60 

61 

f close ( f idl ) ; 

62 

f close ( f id2 ) ; 

63 

end 

i  function  run_kmeans (map.number ,  center.count) 
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heuristics  =  {'lloyd',  'swap',  'hybrid',  ' EZ-hybrid ' } ; 

output  =  sprintf ( ' . /output/output%d ' ,  map_number) ; 

for  h  =  1 : numel (heuristics ) 

input  =  sprintf ('. /commands/commands%d_%d_%s '  ,  ... 

map.number,  center.count,  char (heuristics (h) )) ; 

dos ( sprintf (' kmltest . exe  <  %s  »  %s',  input,  output)); 

end 

end 


function  map_assigned  =  set-assignments  (n_nodes ,  map,  ... 
assignments,  v_n_centers) 
n.heuristics  =  4; 

map_assigned  =  zeros  (numel  (v_n_centers)  ,  n_heuristics ,  3, 
nmodes )  ; 


for  j  =  1 : numel (v_n_centers ) 
for  i  =  1 : n_heuristics 
temp_map  =  map; 

temp_assignments  =  assignments ( j ,  i,  :); 
temprmap ( : , 3 )  =  temp_assignments  (  :  ) ; 


temp_map ( : ,  : ) 

temp_map_assigned  =  map_assigned ( j ,  i,  :,  :) 

map_assigned ( j ,  i,  :,  :)  =  temp_map ( : ,  :); 


for  k  =  l:n_nodes 

map.assigned ( j ,  i,  1,  k)  =  temp_map(k,  1); 

map.assigned  ( j ,  i,  2,  k)  =  temp_map(k,  2)  ; 

map_assigned ( j ,  i,  3,  k)  =  temprmap (k,  3) ; 

end 

end 

end 

end 
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Appendix  G:  VHDL  Code 


VHDL  files  represent  the  FPGA-based  circuitry  used  to  quickly  generate  an  adaptive 
hopset.  The  overall  hopset  selector  is  wired  together  in  adaptive  Jhopset  selector .  vhd, 
and  actual  hopset  selection  is  performed  in  hopset_selector_b .  vhd.  Files  beginning 
with  “tb  ”  are  test  bench  files  for  verifying  system  and  subsystem  functionality.  A  total  of 
1,430  lines  of  VHDL  code  were  implemented  for  this  research. 


1  library  ieee; 

2  use  ieee . stcLlogic_l 1 64 . all ; 

3  use  ieee . numeric_std . all ; 

4 

5  entity  adaptive_hopset_selector  is 

6  port ( 


7 

signal 

elk,  reset 

in 

std_logic ; 

8 

signal 

load_rem 

in 

std_logic ; 

9 

signal 

load-key 

in 

std_logic ; 

10 

signal 

new_rems 

in 

std_logic ; 

11 

signal 

new_key 

in 

std_logic ; 

12 

signal 

input 

in 

std_logic_vector  (31 

downto 

0)  ; 

13 

signal 

key 

in 

std_logic_vector  (31 

downto 

0)  ; 

14 

signal 

lower_bw_mask 

in 

std_logic_vector  (31 

downto 

0)  ; 

15 

signal 

upper_bw_mask 

in 

std_logic_vector  (31 

downto 

0)  ; 

16 

signal 

start 

in 

std_logic ; 

17 

signal 

next-hop 

in 

std_logic ; 

18 

signal 

hopset_ready 

out 

std_logic ; 

19 

signal 

finished 

out 

std_logic ; 

20 

signal 

hop_number 

out 

std_logic_vector  (10 

downto 

0)  ; 

21 

signal 

hop_channel 

out 

std_logic_vector  (10 

downto 

0) 

22 

) ; 

23  end  entity; 

24 

25  architecture  structural  of  adaptive_hopset_selector  is 

26  component  rem_merger  is 

27  port  ( 


28 

signal 

elk,  reset 

in 

std_logic ; 

29 

signal 

load_rem 

in 

std_logic ; 

30 

signal 

new_rems 

in 

std_logic ; 

31 

signal 

input 

in 

std_logic_vector  (31 

downto 

0)  ; 

32 

signal 

rem.rdat 

in 

std_logic_vector  (31 

downto 

0)  ; 

33 

signal 

rem_rsel 

out 

std_logic_vector  (5 

downto 

0)  ; 

34 

signal 

rem-wen 

out 

std_logic; 

35 

signal 

rem_wsel 

out 

std_logic_vector  (5 

downto 

0)  ; 

36 

signal 

rem.wdat 

out 

std_logic_vector  (31 

downto 

0) 

37  )  ; 
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38 

end  component; 

39 

40 

component 

cey_loader  is 

41 

port  ( 

42 

signal 

elk,  reset  : 

43 

signal 

load_key  : 

44 

signal 

new.key  : 

45 

signal 

input  : 

46 

signal 

key_wen  : 

47 

signal 

key_wsel  : 

48 

signal 

key_wdat  : 

49 

) ; 

50 

end  component; 

51 

52 

component 

hop set _s elector  _b 

53 

port  ( 

54 

signal 

elk,  reset  : 

55 

signal 

lower_bw_mask  : 

56 

signal 

upper_bw_mask  : 

57 

signal 

start  : 

58 

signal 

next-hop  : 

59 

signal 

new.rems  : 

60 

signal 

new_key  : 

61 

signal 

hopset.ready  : 

62 

signal 

finished  : 

63 

signal 

hop.number  : 

64 

signal 

hop_channel  : 

65 

66 

signal 

rem.rdat  : 

67 

signal 

rem_rsel  : 

68 

69 

signal 

key_rdat  : 

70 

signal 

key_rsel  : 

71 

72 

signal 

bram_write_en  : 

73 

signal 

bram_read_en  : 

74 

signal 

bram_fill  : 

75 

signal 

bram.addr  : 

0) 

:=  (others  =>  '0 

76 

signal 

bram_data_in  : 

:  = 

(others  =>  '  0  '  )  ; 

77 

signal 

bram_data_out  : 

(others  =>  ' 0 ' ) 

78 

) ; 

79 

end  component; 

80 

81 

component 

:>ram_s  ha  red_2  port 

82 

port  ( 

83 

signal  elk,  reset  : 

84 

signal  wdat  : 

85 

signal  wsel  : 

86 

signal  wen  : 

in  stcLlogic; 
in  stcLlogic; 
in  stcLlogic; 

in  std_logic_vector (31  downto  0)  ; 
out  std_logic; 

out  std_logic_vector (5  downto  0)  ; 
out  std_logic_vector (31  downto  0) 


is 

in  std-logic; 

in  std_logic_vector (31  downto  0); 

in  stddogic-vector  (31  downto  0); 

in  std_logic; 

in  std_logic; 

in  std.logic; 

in  std.logic; 

out  std.logic; 

out  std_logic; 

out  std_logic_vector (10  downto  0); 
out  std_logic_vector (10  downto  0); 

in  std_logic_vector (31  downto  0); 
out  stddogic-vector  (5  downto  0); 

in  std_logic_vector (31  downto  0); 
out  std_logic_vector (5  downto  0); 

buffer  std_logic  :=  'O'; 
out  std_logic  :=  'O'; 
out  std_logic; 

buffer  std_logic_vector ( 10  downto  .. 

r 

out  std_logic_vector (10  downto  0)  .. 

in  std-logic_vector ( 10  downto  0)  := 


is 


in 

std_logic; 

in 

std_logic_vector 

(31  downto  0) 

in 

std_logic -vector 

( 5  downto  0 ) ; 

in 

std_logic; 
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87 

signal 

rsell 

:  in 

std_logic_vector 

(5 

downto 

0)  ; 

88 

signal 

rsel2 

:  in 

std_logic_vector 

(5 

downto 

0)  ; 

89 

signal 

rdatl 

:  out 

std_logic_vector 

(31 

downto 

0)  ; 

90 

signal 

rdat2 

:  out 

std_logic_vector 

(31 

downto 

0) 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 


end  component; 

component  bram_channels_table  is 
port  ( 

signal  clka,  rsta  :  in  std_logic; 


signal  wea 
signal  ena  :  in 
signal  regcea  : 
signal  addra 
signal  dina 
signal  douta 


:  in  std_logic_vector  ( 0  downto  0 
std.logic; 
in  std.logic; 

in  std_logic_vector ( 10  downto 
in  std_logic_vector ( 10  downto 
out  std_logic_vector ( 10  downto 


end  component; 


106 

signal 

remrwdat 

stddogic-vector 

(31 

downto  0) ; 

107 

signal 

rem_wsel 

std_logic_vector 

(5 

downto  0) ; 

108 

signal 

rem.wen 

std_logic; 

109 

signal 

rem.rsell 

std_logic_vector 

(5 

downto  0) ; 

110 

signal 

rem_rsel2 

std_logic_vector 

(5 

downto  0) ; 

111 

signal 

rerrurdat  1 

stddogic-vector 

(31 

downto  0)  ; 

112 

signal 

rem_rdat2 

std_logic_vector 

(31 

downto  0)  ; 

1 13 

114 

signal 

key_wdat 

std_logic_vector 

(31 

downto  0)  ; 

115 

signal 

key_wsel 

std_logic_vector 

(5 

downto  0)  ; 

116 

signal 

key_wen 

std_logic; 

117 

signal 

key_rsel 1 

stddogic-vector 

(5 

downto  0) ; 

118 

signal 

key_rsel2 

std_logic_vector 

(5 

downto  0) ; 

119 

signal 

key.rdat 1 

std_logic_vector 

(31 

downto  0) ; 

120 

signal 

key_rdat2 

std_logic_vector 

(31 

downto  0) ; 

121 

122 

signal 

bram_write_en 

stddogic-vector  ( 0  downto  0)  :  = 

(others  =>  ' 0 ' ) ; 

123 

signal 

bram_read_en 

std_logic  : =  'O'; 

124 

signal 

bram_f ill 

std_logic  : =  'O'; 

125 

signal 

bram.addr 

std_logic_vector  (10 

downto  0)  : 

(others  =>  ' 0 ' ) ; 

126 

signal 

bram_data_in 

std_logic_vector  (10 

downto  0)  : 

(others  =>  ' 0 ' ) ; 

127 

signal 

bram_data_out 

std_logic_vector  (10 

downto  0)  : 

(others  =>  ' 0 ' ) ; 

128  begin 

129  key_rsel2  <  (others  =>  'Z'); 

130  key_rdat2  <  (others  =>  'Z'); 

131 


132  merger:  rem_merger 

133  port  map  ( 

134  elk  =>  elk. 
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135 


reset  =>  reset, 

136  loacLrem  =>  loacLrem, 

137  new_rems  =>  new_rems, 

138  input  =>  input, 

139  rem_rdat  =>  rem.rdatl, 

140  rem.rsel  =>  rem.rsell, 

141  rem_wen  =>  rem.wen, 

142  rem.Msel  =>  rem.wsel, 

143  rem_wdat  =>  rem_wdat 

144  )  ; 

145 

146  loader:  key_loader 

147  port  map  ( 

148  elk  =>  elk, 

149  reset  =>  reset, 

150  load-key  =>  load-key, 

151  new-key  =>  new_key, 

152  input  =>  key, 

153  key_wen  =>  key_wen, 

154  key_wsel  =>  key_wsel, 

155  key_wdat  =>  key.wdat 

156  )  ; 

157 

158  selector:  hopset_selector_b 

159  port  map  ( 

160  elk  =>  elk, 

161  reset  =>  reset, 

162  lower_bw_mask  =>  lower_bw_mask , 

163  upper_bw_mask  =>  upper_bw_mask , 

164  start  =>  start, 

165  next-hop  =>  next-hop, 

166  new_rems  =>  new.rems, 

167  new-key  =>  new_key, 

168  bram-fill  =>  bram_fill, 

169  hopset_ready  =>  hopset_ready , 

no  finished  =>  finished, 

m  hop.number  =>  hop.number, 

172  hop_channel  =>  hop_channel, 

173 

174  rem.rdat  =>  rem_rdat2, 

175  rem_rsel  =>  rem_rsel2, 

176 

177  key_rdat  =>  key_rdatl, 

ns  key_rsel  =>  key_rsell, 

179 

iso  bram-wr ite_en  =>  bram_write_en  ( 0 )  , 

181  bram_read_en  =>  bram_read_en , 

182  bram.addr  =>  bram.addr, 

183  bram_data_in  =>  bram_data_in , 

184  bram_data_out  =>  bram_data_out 

185  )  ; 

186 
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187  bram_rem:  bram_shared_2port 

188  port  map  ( 

189  elk  =>  elk, 

190  reset  =>  new_rems, 

191  wdat  =>  rem_wdat, 

192  wsel  =>  rem.wsel, 

193  wen  =>  rem_wen, 

194  rsell  =>  rem.rsell, 

195  rsel2  =>  rem_rsel2, 

196  rdatl  =>  rem_rdatl, 

197  rdat2  =>  rem_rdat2 

198  )  ; 

199 

200  bram_key:  bram_shared_2port 

201  port  map  ( 

202  elk  =>  elk, 

203  reset  =>  new_key, 

204  wdat  =>  key_wdat, 

205  wsel  =>  keyrwsel, 

206  wen  =>  keyrwen, 

207  rsell  =>  key_rsell, 

208  rsel2  =>  key_rsel2, 

209  rdatl  =>  key_rdatl, 

210  rdat2  =>  key_rdat2 

211  )  ; 

212 

213  bram:  bram_channels_table 

214  port  map  ( 

215  elka  =>  elk, 

216  rsta  =>  new_rems, 

217  wea  =>  bramrwr  ite_en , 

218  ena  =>  bram_read_en , 

219  regcea  =>  bram_fill, 

220  addra  =>  bram.addr, 

221  dina  =>  bram_data_in , 

222  douta  =>  bram_data_out 

223  )  ; 

224  end  structural; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


library  ieee; 

use  ieee . std_logic_l 1 64 . all ; 

use  ieee . numeric.std . all ; 

use  ieee . std_logic_unsigned . all; 


entity  bram_channels_table  is 
—  port ( 

—  signal  elk,  reset  :  in  std_logic; 

—  signal  write.en  :  in  std_logic; 

—  signal  read_en  :  in  std_logic; 

—  signal  fill  :  in  std_logic; 

—  signal  addr  :  in  integer; 
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13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 


—  signal  data_in  :  in  std_logic_vector  (10  downto  0); 

—  signal  data_out  :  out  std_logic_vector (10  downto  0) 

—  ) ; 

—  end  entity  bram_channels_table ; 

entity  bram_channels_table  is 
port  ( 

clka  :  in  std_logic; 
rsta  :  in  std_logic; 
ena  :  in  std_logic; 
regcea  :  in  std.logic; 

wea  :  in  std_logic_vector ( 0  downto  0); 
addra  :  in  std_logic_vector ( 10  downto  0); 
dina  :  in  std_logic_vector ( 10  downto  0); 
douta  :  out  std_logic_vector (10  downto  0) 

) ; 

end  bram_channels_table; 

architecture  structural  of  bram_channels_table  is 

type  channel-list  is  array  (0  to  2048)  of  std_logic_vector ( 10 
downto  0 ) ; 

signal  channels_table  :  channel_list ; 

signal  open.channels  :  integer  :=  0; 

signal  fill  :  std_logic  :=  'O'; 
signal  fill_table  :  std_logic  :=  'O'; 
signal  fill_start  :  std_logic  :=  'O'; 
signal  low_offset  :  integer  :=  0; 
signal  high_offset  :  integer  :=  0; 
signal  increment  :  integer  :=  0; 

signal  write_en_buf f er  :  std_logic  :=  'O'; 
signal  read_en_buf f er  :  std_logic  :=  'O'; 
begin 

channels_table (2048)  <  (others  =>  '0'); 

fill  <  regcea; 

bram:  process (clka,  rsta,  wea,  ena,  fill) 
begin 

if (rsta  =  '1')  then 

wr ite_en_buf f er  <  'O'; 

read_en_buf f er  <  'O'; 

open.channels  <  0; 
douta  <  (others  =>  '0'); 
elsif (rising_edge (clka) )  then 
write_en_buf fer  <  wea(0); 
read_en_buf f er  <  ena; 
open.channels  <  open.channels; 

if (wea (0)  =  not  write_en_buf f er )  then 


elk 

reset 

enable 

fill 

write  enable 
address 
data  in 
data  out 
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63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 


channels.table  (to.integer  (unsigned  (addra)  )  )  <  dina(10 
down to  0) ; 

open_channels  <  open_channels  +  1; 
end  if; 

elsif  (fill  =  '1')  then 
f  ill.table  <  '  1'; 
f ill_st art  <  ' 1 '  ; 
elsif ( fill.table  =  ' 1')  then 
if ( f ill_start  =  '  l1)  then 
low.offset  <  open.channels; 

high_offset  <  open_channels+open_channels  — 1; 
increment  <  open.channels; 
f  ill.start  <  '  0  f  ; 

elsif ( f ill.start  =  f0'  and  fill.table  =  '1’)  then 
if (high.offset  <  2047)  then 

channels.table  ( low.offset  to  high.offset)  <  ... 

channels.table  ( 0  to  open.channels  — 1)  ; 
low.offset  <  low.offset  +  increment; 
high.offset  <  high.offset  +  increment; 
else 

f ill.table  <  'O'; 

channels.table  ( low.offset  to  2047)  <  channels.table  ( 0 
to  (2047-low_of f set ) ) ; 
end  if; 
else 

low.offset  <  0; 
high.offset  <  0; 
increment  <  0; 
f ill.table  <  '0f; 
end  if; 

elsif ( rising.edge (ena)  and  f ill.table  =  '0f)  then 

douta  <  channels.table (to.integer (unsigned (addra) ) ) (10  .., 
downto  0) ; 
end  if; 

end  process  bram; 
end  structural; 


library  ieee; 

use  ieee . std.logic.l 1 64 . all; 

entity  bram_shared_2port  is 
port  ( 


signal 

elk. 

reset 

:  in  std.logic; 

signal 

wdat 

:  in 

std_logic_vector  (31 

downto  0 ) ; 

signal 

wsel 

:  in 

std_logic_vector  (5 

downto  0) ; 

signal 

wen 

:  in 

std_logic; 

signal 

rsell 

in  std_logic_vector 

(5  downto  0) ; 

signal 

rsel2 

in  std_logic_vector 

(5  downto  0) ; 

signal 

rdatl 

out  stddogic-vector 

(31  downto  0) 

signal 

rdat2 

out  stddogic-vector 

(31  downto  0) 
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15  end  entity; 

16 

17  architecture  structural  of  bram_shared_2port  is 
is  type  shmem64x32  is  array  (0  to  63)  of  std_logic_vector  (31 

downto  0 ) ; 

19  signal  mem  :  shmem64x32; 

20  begin 

21  shmem:  process  (elk,  reset) 

22  begin 


23 

if (reset 

= 

' 1 ' )  then 

24 

mem ( 0 ) 

< 

(others 

=> 

1') 

f 

25 

mem ( 1 ) 

< 

(others 

=> 

1') 

t 

26 

mem ( 2 ) 

< 

(others 

=> 

1') 

t 

27 

mem ( 3 ) 

< 

(others 

=> 

1') 

t 

28 

mem ( 4 ) 

< 

(others 

=> 

1') 

t 

29 

mem ( 5 ) 

< 

(others 

=> 

1') 

t 

30 

mem ( 6 ) 

< 

(others 

=> 

1' 

r 

31 

mem ( 7 ) 

< 

(others 

=> 

1' 

r 

32 

mem ( 8 ) 

< 

(others 

=> 

1' 

r 

33 

mem ( 9 ) 

< 

(others 

=> 

1') 

r 

34 

mem (10) 

< 

(others 

=> 

'  1 

) 

35 

mem (11) 

< 

(others 

=> 

'  1 

) 

36 

mem (12) 

< 

(others 

=> 

'  1 

) 

37 

mem (13) 

< 

(others 

=> 

'  1 

) 

38 

mem (14) 

< 

(others 

=> 

'  1 

) 

39 

mem (15) 

< 

(others 

=> 

'  1 

) 

40 

mem (16) 

< 

(others 

=> 

'  1 

) 

41 

mem (17) 

< 

(others 

=> 

'  1 

) 

42 

mem (18) 

< 

(others 

=> 

'  1 

) 

43 

mem (19) 

< 

(others 

=> 

'  1 

) 

44 

mem (20) 

< 

(others 

=> 

'  1 

) 

45 

mem (21) 

< 

(others 

=> 

'  1 

) 

46 

mem (22) 

< 

(others 

=> 

'  1 

) 

47 

mem (23) 

< 

(others 

=> 

'  1 

) 

48 

mem (24) 

< 

(others 

=> 

'  1 

) 

49 

mem (25) 

< 

(others 

=> 

'  1 

) 

50 

mem (26) 

< 

(others 

=> 

'  1 

) 

51 

mem (27) 

< 

(others 

=> 

'  1 

) 

52 

mem (28) 

< 

(others 

=> 

'  1 

) 

53 

mem (29) 

< 

(others 

=> 

'  1 

) 

54 

mem ( 30 ) 

< 

(others 

=> 

'  1 

) 

55 

mem ( 31 ) 

< 

(others 

=> 

'  1 

) 

56 

mem ( 32 ) 

< 

(others 

=> 

'  1 

) 

57 

mem (33) 

< 

(others 

=> 

'  1 

) 

58 

mem (34) 

< 

(others 

=> 

'  1 

) 

59 

mem (35) 

< 

(others 

=> 

'  1 

) 

60 

mem (36) 

< 

(others 

=> 

'  1 

) 

61 

mem ( 37 ) 

< 

(others 

=> 

'  1 

) 

62 

mem ( 38 ) 

< 

(others 

=> 

'  1 

) 

63 

mem (39) 

< 

(others 

=> 

'  1 

) 

64 

mem (40) 

< 

(others 

=> 

'  1 

) 

65 

mem (41) 

< 

(others 

=> 

'  1 

) ; 
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66 

mem (42)  < 

(others 

=> 

'i') ; 

67 

mem (43)  < 

(others 

=> 

'i') ; 

68 

mem (44)  < 

(others 

=> 

'i') ; 

69 

mem (45)  < 

(others 

=> 

'i') ; 

70 

mem (46)  < 

(others 

=> 

'i') ; 

71 

mem (47)  < 

(others 

=> 

’i'); 

72 

mem (48)  < 

(others 

=> 

’i'); 

73 

mem (49)  < 

(others 

=> 

'l'); 

74 

mem (50)  < 

(others 

=> 

'l'); 

75 

mem (51)  < 

(others 

=> 

’l'); 

76 

mem (52)  < 

(others 

=> 

'i') ; 

77 

mem (53)  < 

(others 

=> 

'i') ; 

78 

mem (54)  < 

(others 

=> 

'i') ; 

79 

mem (55)  < 

(others 

=> 

'i') ; 

80 

mem (56)  < 

(others 

=> 

'i') ; 

81 

mem (57)  < 

(others 

=> 

'i') ; 

82 

mem (58)  < 

(others 

=> 

’l'); 

83 

mem (59)  < 

(others 

=> 

’l'); 

84 

mem (60)  < 

(others 

=> 

’l'); 

85 

mem (61)  < 

(others 

=> 

’l'); 

86 

mem (62)  < 

(others 

=> 

’l'); 

87 

mem (63)  < 

(others 

=> 

'l') ; 

88 

elsif (rising_edge (elk) 

then 

89 

if (wen  = 

'  1  '  )  then 

90 

case  wsel  is 

91 

when 

"000000" 

=> 

mem ( 0 ) 

<  wdat ; 

92 

when 

"000001" 

=> 

mem ( 1 ) 

<  wdat ; 

93 

when 

"000010" 

=> 

mem ( 2 ) 

<  wdat ; 

94 

when 

"000011" 

=> 

mem ( 3 ) 

<  wdat ; 

95 

when 

"000100" 

=> 

mem ( 4 ) 

<  wdat ; 

96 

when 

"000101" 

=> 

mem ( 5 ) 

<  wdat ; 

97 

when 

"000110" 

=> 

mem ( 6 ) 

<  wdat ; 

98 

when 

"000111" 

=> 

mem ( 7 ) 

<  wdat ; 

99 

when 

"001000" 

=> 

mem ( 8 ) 

<  wdat ; 

100 

when 

"001001" 

=> 

mem ( 9 ) 

<  wdat ; 

101 

when 

"001010" 

=> 

mem (10) 

<  wdat ; 

102 

when 

"001011" 

=> 

mem (11) 

<  wdat ; 

103 

when 

"001100" 

=> 

mem ( 12 ) 

<  wdat ; 

104 

when 

"001101" 

=> 

mem (13) 

<  wdat ; 

105 

when 

"001110" 

=> 

mem (14) 

<  wdat ; 

106 

when 

"001111" 

=> 

mem (15) 

<  wdat ; 

107 

when 

"010000" 

=> 

mem (16) 

<  wdat ; 

108 

when 

"010001" 

=> 

mem (17) 

<  wdat ; 

109 

when 

"010010" 

=> 

mem (18) 

<  wdat ; 

110 

when 

"010011" 

=> 

mem (19) 

<  wdat ; 

111 

when 

"010100" 

=> 

mem (20) 

<  wdat ; 

112 

when 

"010101" 

=> 

mem (21) 

<  wdat ; 

113 

when 

"010110" 

=> 

mem (22 ) 

<  wdat ; 

114 

when 

"010111" 

=> 

mem (23) 

<  wdat ; 

115 

when 

"011000" 

=> 

mem (24) 

<  wdat ; 

116 

when 

"011001" 

=> 

mem (25) 

<  wdat ; 

117 

when 

"011010" 

=> 

mem (26) 

<  wdat ; 

Ill 


118 

when 

"011011" 

=> 

mem 

(27) 

< 

wdat ; 

119 

when 

"011100" 

=> 

mem 

(28) 

< 

wdat ; 

120 

when 

"011101" 

=> 

mem 

(29) 

< 

wdat ; 

121 

when 

"011110" 

=> 

mem 

(30) 

< 

wdat ; 

122 

when 

"011111" 

=> 

mem 

(31) 

< 

wdat ; 

123 

when 

"100000" 

=> 

mem 

(32) 

< 

wdat ; 

124 

when 

"100001" 

=> 

mem 

(33) 

< 

wdat ; 

125 

when 

"100010" 

=> 

mem 

(34) 

< 

wdat ; 

126 

when 

"100011" 

=> 

mem 

(35) 

< 

wdat ; 

127 

when 

"100100" 

=> 

mem 

(36) 

< 

wdat ; 

128 

when 

"100101" 

=> 

mem 

(37) 

< 

wdat ; 

129 

when 

"100110" 

=> 

mem 

(38) 

< 

wdat ; 

130 

when 

"100111" 

=> 

mem 

(39) 

< 

wdat ; 

131 

when 

"101000" 

=> 

mem 

(40) 

< 

wdat ; 

132 

when 

"101001" 

=> 

mem 

(41) 

< 

wdat ; 

133 

when 

"101010" 

=> 

mem 

(42) 

< 

wdat ; 

134 

when 

"101011" 

=> 

mem 

(43) 

< 

wdat ; 

135 

when 

"101100" 

=> 

mem 

(44) 

< 

wdat ; 

136 

when 

"101101" 

=> 

mem 

(45) 

< 

wdat ; 

137 

when 

"101110" 

=> 

mem 

(46) 

< 

wdat ; 

138 

when 

"101111" 

=> 

mem 

(47) 

< 

wdat ; 

139 

when 

"110000" 

=> 

mem 

(48) 

< 

wdat ; 

140 

when 

"110001" 

=> 

mem 

(49) 

< 

wdat ; 

141 

when 

"110010" 

=> 

mem 

(50) 

< 

wdat ; 

142 

when 

"110011" 

=> 

mem 

(51) 

< 

wdat ; 

143 

when 

"110100" 

=> 

mem 

(52) 

< 

wdat ; 

144 

when 

"110101" 

=> 

mem 

(53) 

< 

wdat ; 

145 

when 

"110110" 

=> 

mem 

(54) 

< 

wdat ; 

146 

when 

"110111" 

=> 

mem 

(55) 

< 

wdat ; 

147 

when 

"111000" 

=> 

mem 

(56) 

< 

wdat ; 

148 

when 

"111001" 

=> 

mem 

(57) 

< 

wdat ; 

149 

when 

"111010" 

=> 

mem 

(58) 

< 

wdat ; 

150 

when 

"111011" 

=> 

mem 

(59) 

< 

wdat ; 

151 

when 

"111100" 

=> 

mem 

(60) 

< 

wdat ; 

152 

when 

"111101" 

=> 

mem 

(61) 

< 

wdat ; 

153 

when 

"111110" 

=> 

mem 

(62) 

< 

wdat ; 

154 

when 

others  => 

mem (63)  < 

wdat ; 

155 

end  case; 

156 

end  if; 

157 

end  if; 

158 

end  process  i 

shmem; 

159 

160 

with  rsell  select 

161 

rdatl  <  mem(0)  when 

"  0  0  0  0  0 ' 

0", 

162 

mem ( 1 )  when 

1! 

0000 

01", 

163 

mem (2)  when 

II 

0000 

10", 

164 

mem (3)  when 

II 

0000 

11", 

165 

mem ( 4 )  when 

II 

0001 

00", 

166 

mem (5)  when 

II 

0001 

01", 

167 

mem (6)  when 

II 

0001 

10", 

168 

mem (7)  when 

II 

0001 

11", 

169 

mem (8)  when 

II 

0010 

00", 

178 


170 

mem  i 

(9) 

when  " 

001001" 

171 

mem  i 

(10) 

when 

"001010 

172 

mem  i 

(11) 

when 

"001011 

173 

mem  i 

(12) 

when 

"001100 

174 

mem  i 

(13) 

when 

"001101 

175 

mem  i 

(14) 

when 

"001110 

176 

mem  i 

(15) 

when 

"001111 

111 

mem  i 

(16) 

when 

"010000 

178 

mem  i 

(17) 

when 

"010001 

179 

mem  i 

(18) 

when 

"010010 

180 

mem  i 

(19) 

when 

"010011 

181 

mem  i 

(20) 

when 

"010100 

182 

mem  i 

(21) 

when 

"010101 

183 

mem  i 

(22) 

when 

"010110 

184 

mem  i 

(23) 

when 

"010111 

185 

mem  i 

(24) 

when 

"011000 

186 

mem  i 

(25) 

when 

"011001 

187 

mem  i 

(26) 

when 

"011010 

188 

mem  i 

(27) 

when 

"011011 

189 

mem  i 

(28) 

when 

"011100 

190 

mem  i 

(29) 

when 

"011101 

191 

mem  i 

(30) 

when 

"011110 

192 

mem  i 

(31) 

when 

"011111 

193 

mem  i 

(32) 

when 

"100000 

194 

mem  i 

(33) 

when 

"100001 

195 

mem  i 

(34) 

when 

"100010 

196 

mem  i 

(35) 

when 

"100011 

197 

mem  i 

(36) 

when 

"100100 

198 

mem  i 

(37) 

when 

"100101 

199 

mem  i 

(38) 

when 

"100110 

200 

mem  i 

(39) 

when 

"loom 

201 

mem  i 

(40) 

when 

"101000 

202 

mem  i 

(41) 

when 

"101001 

203 

mem  i 

(42) 

when 

"101010 

204 

mem  i 

(43) 

when 

"101011 

205 

mem  i 

(44) 

when 

"101100 

206 

mem  i 

(45) 

when 

"101101 

207 

mem  i 

(46) 

when 

"101110 

208 

mem  i 

(47) 

when 

"101111 

209 

mem  i 

(48) 

when 

"110000 

210 

mem  i 

(49) 

when 

"110001 

211 

mem  i 

(50) 

when 

"110010 

212 

mem  i 

(51) 

when 

"110011 

213 

mem  i 

(52) 

when 

"110100 

214 

mem  i 

(53) 

when 

"110101 

215 

mem  i 

(54) 

when 

"110110 

216 

mem  i 

(55) 

when 

"110111 

217 

mem  i 

(56) 

when 

"111000 

218 

mem  i 

(57) 

when 

"111001 

219 

mem  i 

(58) 

when 

"111010 

220 

mem  i 

(59) 

when 

"111011 

221 

mem  i 

(60) 

when 

"111100 

179 


222 

mem ( 61 

when 

"111101 

223 

mem ( 62 

when 

"111110 

224 

mem ( 63 

when 

others ; 

225 

226 

with  rsel2 

select 

227 

rdat2  < 

mem(O)  when  "000000", 

228 

mem 

1) 

when  " 

000001" 

229 

mem 

2) 

when  " 

000010" 

230 

mem 

3) 

when  " 

000011" 

231 

mem 

4) 

when  " 

000100" 

232 

mem 

5) 

when  " 

000101" 

233 

mem 

6) 

when  " 

000110" 

234 

mem 

7) 

when  " 

000111" 

235 

mem 

8) 

when  " 

001000" 

236 

mem 

9) 

when  " 

001001" 

237 

mem 

10 

when 

"001010 

238 

mem 

11 

when 

"001011 

239 

mem 

12 

when 

"001100 

240 

mem 

13 

when 

"001101 

241 

mem 

14 

when 

"001110 

242 

mem 

15 

when 

"001111 

243 

mem 

16 

when 

"010000 

244 

mem 

17 

when 

"010001 

245 

mem 

18 

when 

"010010 

246 

mem 

19 

when 

"010011 

247 

mem 

20 

when 

"010100 

248 

mem 

21 

when 

"010101 

249 

mem 

22 

when 

"010110 

250 

mem 

23 

when 

"010111 

251 

mem 

24 

when 

"011000 

252 

mem 

25 

when 

"011001 

253 

mem 

26 

when 

"011010 

254 

mem 

27 

when 

"011011 

255 

mem 

28 

when 

"011100 

256 

mem 

29 

when 

"011101 

257 

mem 

30 

when 

"011110 

258 

mem 

31 

when 

"011111 

259 

mem 

32 

when 

"100000 

260 

mem 

33 

when 

"100001 

261 

mem 

34 

when 

"100010 

262 

mem 

35 

when 

"100011 

263 

mem 

36 

when 

"100100 

264 

mem 

37 

when 

"100101 

265 

mem 

38 

when 

"100110 

266 

mem 

39 

when 

"100111 

267 

mem 

40 

when 

"101000 

268 

mem 

41 

when 

"101001 

269 

mem 

42 

when 

"101010 

270 

mem 

43 

when 

"101011 

271 

mem 

44 

when 

"101100 

272 

mem 

45 

when 

"101101 

273 

mem 

46 

when 

"101110 

180 


274 

mem 

47) 

when 

"101111 

275 

mem 

48) 

when 

"110000 

276 

mem 

49) 

when 

"110001 

277 

mem 

50) 

when 

"110010 

278 

mem 

51) 

when 

"110011 

279 

mem 

52) 

when 

"110100 

280 

mem 

53) 

when 

"110101 

281 

mem 

54) 

when 

"110110 

282 

mem 

55) 

when 

"110111 

283 

mem 

56) 

when 

"111000 

284 

mem 

57) 

when 

"111001 

285 

mem 

58) 

when 

"111010 

286 

mem 

59) 

when 

"111011 

287 

mem 

60) 

when 

"111100 

288 

mem 

61) 

when 

"111101 

289 

mem 

62) 

when 

"111110 

290 

mem 

63) 

when 

others ; 

291 

end  structural; 

library  ieee; 

use  ieee . std_logic_l  1 64 . all ; 

use  ieee . numeric.std . all ; 

use  ieee. std_logic_unsigned. all; 

entity  hopset_selector_b  is 
port  ( 

signal  elk,  reset  :  in  std_logic; 

signal  lower_bw_mask  :  in  std_logic_vector  (31  downto  0) 

signal  upper_bw_mask  :  in  std_logic_vector  (31  downto  0) 

signal  start  :  in  std_logic; 

signal  next_hop  :  in  std_logic; 

signal  new.reras  :  in  std_logic; 

signal  new_key  :  in  std_logic; 

signal  hopset_ready  :  out  std_logic; 

signal  finished  :  out  std_logic; 

signal  hop.number  :  out  std_logic_vector (10  downto  0); 
signal  hop_channel  :  out  std_logic_vector  ( 10  downto  0); 

signal  rem_rdat  :  in  std_logic_vector  (31  downto  0); 

signal  rem_rsel  :  out  std_logic_vector (5  downto  0); 

signal  key_rdat  :  in  std_logic_vector (31  downto  0); 

signal  key_rsel  :  out  std_logic_vector (5  downto  0); 

signal  bram_write_en  :  buffer  std_logic; 
signal  bram_read_en  :  out  std_logic; 
signal  bram_fill  :  out  std_logic; 


29 

signal 

bram.addr  :  buffer 

std_logic_vector 

(10 

downto 

0)  ; 

30 

signal 

bram_data_in  :  out 

stddogic-vector 

(10 

downto 

0)  ; 

31 

signal 

bram_data_out  :  in 

stddogic-vector 

(10 

downto 

0) 

) ; 
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33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 


end  entity; 

architecture  behavioral  of  hopset_selector Jo  is 
—  FSM  states 

type  ahg_state  is  (s_idle,  s.count,  s.ready,  s_finished) ; 
signal  current.state,  next.state  :  ahg_state; 


—  Loading  metrics 

signal  load-index  :  integer  :=  0; 
constant  load-increment  :  integer  :=  32; 
signal  read_count  :  integer  :=  0; 
constant  read-total  :  integer  :=  64; 
signal  key_index_count  :  integer  :=  0; 


—  Aggregate  rems 

signal  rem_buf f er.prev  :  std_logic_vector (31  downto  0) 
(others  =>  ' 0  ' ) ; 

signal  remJouf f er_curr  :  std_logic_vector (31  downto  0) 
(others  =>  ' 0 ' ) ; 

signal  remJouf f errnext  :  std_logic_vector  (31  downto  0) 
(others  =>  '  0  '  )  ; 

signal  remJouf fer_all  :  std_logic_vector ( 95  downto  0) 
(others  =>  ' 0 ' ) ; 

signal  remJouf f er.temp  :  std_logic_vector  ( 64  downto  0) 
(others  =>  ' 0 ' ) ; 

signal  remJouf f er_mask  :  std_logic_vector  ( 64  downto  0) 
(others  =>  ' 0 ' ) ; 


—  Key 

signal  keyJouffer  :  std_logic_vector (31  downto  0)  :=  (others  => 

'  o 1 ) ; 

signal  key_slice  :  std_logic_vector ( 10  downto  0)  :=  (others  =>  . 

’O'); 

—  Counting  metrics 

signal  open_channels  :  integer  :=  0; 
signal  channel-index  :  integer  :=  0; 
signal  count-total  :  integer  :=  0; 

—  Hop  index  tracker 

signal  single_complete  :  std_logic  :=  'O'; 
signal  hop_output_index  :  integer  :=  0; 
signal  hop_output_index_buf fer  :  integer  :=  0; 

signal  run_of_ones_mask  :  std_logic_vector (10  downto  0)  :=  ... 

"11101111110"; 

begin 

key_slice  <  keyJouf fer (10  downto  0) ; 
bram_read_en  <  next-hop; 
hop_channel  <  bram_data_out ; 

remJouf fer_all  <  remJouf f errnext  &  remJouf f erxurr  &  ... 
remJouf  f  errprev; 


182 


75 


rem_buf f er_temp  <  rem_buf fer_all (channel_index+64  downto 
channeLindex)  ; 

76  remlbuf fer_mask  <  upper_bw_mask  &  '  1  '  &  lower _bw_mask; 

77 

78  selector:  process  (elk,  reset,  next-state) 

79  begin 

so  current-state  <  current-state; 

81 

82  if  (reset  =  '1')  then 

83  current-state  <  s_idle; 

84  elsif  (rising_edge  (elk)  )  then 

85  current-state  <  next.state; 

86  end  if; 

87  end  process  selector; 

88 

89  fsm:  process  (elk,  reset,  current-state,  start,  next-hop, 

new_rems,  new_key) 

90  — variable  channel-lookup  :  integer  :=  0; 

91  begin 

92  if  (reset  =  '1')  then 

93  —  State  logic 

94  open.channels  <  0; 

95  channel-index  <  0; 

96  key_buffer  <  (others  =>  '0'); 

97  hop_output_index  <  0; 

98  single-complete  <  'O'; 

99  run_of  _ones_mask  <  "11111111110"; 

100 

101  —  BRAM  signals 

102  bramrwr ite_en  <  'O'; 

103  bram_fill  <  'O'; 

104  bram.addr  <  (others  =>  '0'); 

105  bram_data_in  <  (others  =>  '0'); 

106  rem.rsel  <  (others  =>  '0'); 

107  key_rsel  <  (others  =>  '0'); 

108 

109  —  Output  signals 

no  hopset.ready  <  'O'; 

in  finished  <  'O'; 

112  hop.number  <  (others  =>  '0'); 

113  — hop_channel  <  (others  =>  '  0  ' )  ; 

114 

ns  —  Next  state  logic 

ii6  next_state  <  s_idle; 

m  elsif (rising_edge (elk) )  then 

ns  —  Default  logic 

119  open-channels  <  open-channels  ; 

120  channel-index  <  channel-index; 

121  key_buffer  <  key_buffer; 

122  hop_output_index  <  hop_output_index; 

123  single_complete  <  single_complete ; 

124 
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125 


126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 


case  current_state  is 
when  s_idle  => 

—  Current  state  logic 
open_channels  <  0; 
channel-index  <  0; 
key_buffer  <  (others  =>  f0f); 
single_complete  <  'O'; 
hop_output_index  <  0; 
run_of _ones_mask  <  "11111111110"; 

—  BRAM  signals 
bram.wr it e_en  <  'O'; 

bram_f ill  <  'O'; 
bram_addr  <  (others  =>  ’O'); 
bram_data_in  <  (others  =>  f0'); 

—  Output  signals 
hopset_ready  <  'O'; 
finished  <  '  1 '  ; 
hop_number  <  (others  =>  '0'); 

— hop_channel  <  (others  =>  '0'); 

—  Next  state  logic 
if  (start  =  '1')  then 
next_state  <  s_count; 

rem_rsel  <  std_logic_vector  (to.unsigned  ( read_count  +  l ,  ... 

6)  )  ; 

rem.buf f er.prev  <  rem_buf f er.curr ; 
rem.buf  f  er.curr  <  rem_buf f er_next ; 
rem_buf f er_next  <  rem_rdat; 
else 

next-state  <  s.idle; 
end  if; 

when  s.count  => 

—  Current  state  logic 

rem_rsel  <  std_logic_vector  ( to_unsigned  ( read_count  +  l ,  6)  )  ; 

count-total  <  count-total  +  1; 

if (channel-index  =  31)  then 
channel-index  <  0; 

read-count  <  read-count  +  1; 
rem_buf f er.prev  <  rem_buf f er.curr ; 
rem.buf f er.curr  <  rem_buf f er_next ; 
rem.buf f er.next  <  rem.rdat; 

if  (  ( remlbuf  f  er_temp  and  remlbuf  f  er_mask )  =  ... 
rem.buf  f  er_mask )  then 
open.channels  <  open.channels  +  1; 

bram_wr it e_en  <  not  bram_write_en; 
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174 

175 

176 

177 

178 

179 

180 
181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 
219 


bram_addr  <  ... 

std_logic_vector (to_unsigned (open.channels,  11) ) ; 
— bram_fill  <  'O'; 
bram_data_in  <  ... 

std_logic_vector  (to.unsigned  ( count-total ,  11)  )  ; 
end  if; 

elsif ( read_count  <  64)  then 

channel-index  <  channel-index  +  1; 

if  (  (  rem_buf  f  er.temp  and  rem_buf  f  er_mask )  =  ... 
rem.buf  f  erjnask )  then 
open.channels  <  open.channels  +  1; 

bram.wr  ite.en  <  not  bram.wr ite.en ; 
bram.addr  <  ... 

std.logic.vector (to.unsigned (open.channels,  11) ) ; 
bram.f ill  <  'O'; 
bram.data.in  <  ... 

std.logic.vector  (to.unsigned  ( count  .total ,  11)  )  ; 
end  if; 
end  if; 

—  BRAM  signals 
— bram.write.en  <  'O'; 

— bram.addr  <  0; 

— bram_data_in  <  (others  =>  '0'); 


—  Output  signals 
— hopset.ready  <  'O'; 
finished  <  'O'; 

— hop.number  <  (others  =>  '0'); 
— hop.channel  <  (others  =>  '0'); 


1'; 

0; 


—  Next  state  logic 
if ( count-total  <  2047)  then 
next_state  <  s.count; 
else 

bram.fill  < 
read.count  < 
key.buffer  <  key.rdat; 
next_state  <  s.ready; 
end  if; 

when  s.ready  => 

—  Current  state  logic 
bram.fill  <  'O'; 


if (next.hop  =  '1'  and  single.complete  =  'O'  and  ... 
read.count  <  64)  then 

key.rsel  <  std.logic.vector (to.unsigned (read.count ,  6)); 
hop.output.index  <  hop.output.index  +  1; 
single.complete  <  '1'; 
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if  (key_index_count  =  31)  then 
key_index_count  <  0; 
key_buffer  <  key_rdat; 
read_count  <  read_count  +  1; 
else 

key_index_count  <  key_index_count  +  1; 
key_buffer  <  key_buf f er  ( 30  downto  0)  &  key_buf f er  ( 31 ) 
hop_output_index  <  hop_output_index  +  1; 
end  if; 

if (key_slice  =  "11111111111")  then 

run_of  _ones_mask  <  run.of  _ones_mask  ( 9  downto  0)  &  ... 
run_of _ones_mask  (10)  ; 

bram_addr  <  (key_slice  and  run_of _ones_mask )  xor  ... 
std_logic_vect or  (to_un signed  ( open_channels ,  11)); 

else 

bram.addr  <  key_slice  xor  .  .  . 

std_logic_vect or  (to.unsigned  (open.channels,  11)  )  ; 
end  if; 
end  if; 

if  (next_hop  =  '0')  then 
single_complete  <  'O'; 
end  if; 

—  Output  signals 
hopset_ready  <  f 1 '  ; 
finished  <  '0f; 
hop.number  <  ... 

std_logic_vector  (to_un signed  (hop_out put -index ,  11)); 

—  BRAM  signals 

— bram.wr it e_en  <  'O'; 

— bram.addr  <  0; 

— bram_fill  <  'O'; 

— bram_data_in  <  (others  =>  '0'); 

—  If  all  ones,  rotate  mask  left  one  bit  and  use  mask 

—  Next  state  logic 

if  (hop_output_index  <  2048)  then 
next-state  <  s.ready; 
else 

next.state  <  s_finished; 
end  if; 

when  s_finished  =>  —  Equivalent  to  s_finished 

—  Current  state  logic 

—  (NONE) 

—  Output  signals 
hopset_ready  <  'O'; 
finished  <  ' 1 ' ; 
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— hop.number  <  (others  =>  '0'); 

— hop_channel  <  (others  =>  '0'); 

—  Next  state  logic 

if (new_rems  =  ’ 1'  or  new_key  =  '1')  then 
next_state  <  s_idle; 
else 

next_state  <  s.finished; 
end  if; 
end  case; 
end  if; 

end  process  fsm; 
end  behavioral; 


library  ieee; 

use  ieee  .  std_logic_l  1  64  .  all; 

use  ieee . numeric_std . all; 

use  ieee . std_logic_un signed . all ; 

entity  key_loader  is 
port  ( 

signal  elk,  reset  :  in  std_logic; 
signal  load_key  :  in  std_logic; 
signal  new_key  :  in  std_logic; 

signal  input  :  in  std_logic_vector (31  downto  0); 
signal  key_wen  :  out  std_logic; 

signal  key_wsel  :  out  st d_logic_vect or  ( 5  downto  0)  ; 
signal  key_wdat  :  out  st d_logic_vect or  ( 31  downto  0) 

)  ; 

end  entity; 

architecture  behavioral  of  key_loader  is 
signal  index  :  integer  :=  0; 
begin 

load_maps :  process (elk,  reset) 
begin 

if  (reset  =  ' 1’)  then 
key_wen  <  'O'; 
key_wsel  <  "000000"; 
key_wdat  <  (others  =>  '0'); 

index  <  0; 

elsif ( rising_edge (elk) )  then 
index  <  index; 

if (index  =  64  or  new.key  =  '1')  then 
key_wen  <  'O'; 
keyrwsel  <  "000000"; 
key_wdat  <  (others  =>  '0'); 
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1  '  )  then 


38  index  <  0; 

39  elsif  ( load_key  =  ' 

40  key.wen  <  '  1 '  ; 

41  key.wsel  <  std.logic.vector (to.unsigned (index,  6)  )  ; 

42  key_wdat  <  input; 

43 

44  index  <  index  +  1; 

45  end  if; 

46  end  if; 

47  end  process  load-maps; 

48  end  behavioral; 


1  library  ieee; 

2  use  ieee . std.logic.l  1 64 . all ; 

3  use  ieee . numeric.std . all ; 

4  use  ieee . std_logic_unsigned . all; 

5 

6  entity  rem_merger  is 


7 

port  ( 

8 

signal 

elk,  reset 

:  in  std.logic; 

9 

signal 

load.rem  : 

in  std.logic; 

10 

signal 

new.rems  : 

in  std.logic; 

11 

signal 

input  :  in 

std-logic-vector (31  downto  0) ; 

12 

13 

signal 

rem.rdat  : 

in  std_logic_vector (31 

downto 

0)  ; 

14 

signal 

rem_rsel  : 

out  std_logic_vector  (5 

downto 

0)  ; 

15 

signal 

rem.wen  : 

out  std.logic; 

16 

signal 

rem_wsel  : 

out  std.logic.vector  (5 

downto 

0)  ; 

17 

signal 

rem_wdat  : 

out  std.logic.vector  (31 

downto 

0) 

18 

) ; 

19  end  entity; 

20 

21  architecture  behavioral  of  rem_merger  is 

22  signal  index  :  integer  :=  0; 

23  begin 

24  load_maps  :  process  (elk,  reset) 

25  begin 

26  if  (reset  =  '1')  then 

27  rem.rsel  <  "000000"; 

28  rem.wen  <  'O'; 

29  rem.wsel  <  "000000"; 

30  rem.wdat  <  (others  =>  '1'); 

31 

32  index  <  0; 

33  elsif  (rising.edge  (elk)  )  then 

34  index  <  index; 

35 

36  if  (index  =  64  or  new.rems  =  '1')  then 

37  rem_rsel  <  "000000"; 

38  rem-wen  <  'O'; 

39  rem_wsel  <  "000000"; 
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rem.wdat 


(others  =>  '  0  '  )  ; 


index  <  0; 

elsif ( load_rem  =  '1')  then 

rem.rsel  <  std_logic_vector  (to_unsigned  (index+1,  6)); 
rem.wen  <  '  1 '  ; 

rem.wsel  <  std_logic_vector  (to_unsigned  ( index,  6)); 
rem.wdat  <  rem.rdat  and  input; 

index  <  index  +  1; 
end  if; 
end  if; 

end  process  load_maps; 
end  behavioral; 


library  ieee; 

use  ieee  .  std_logic_l  1  64  .  all; 
use  ieee . numeric_std . all; 
use  std. textio . all; 

use  ieee . math_real . all ;  —  for  UNIFORM,  TRUNC  functions 

entity  tb_adapt ive_hopset_selector_b  is 
generic ( 

period:  time  :=  10  ns 

) ; 

end  entity; 


architecture  behavioral  of  tb_adapt ive_hopset_selector_b  is 
—  Component  signals 
signal  elk,  reset  :  std.logic; 
signal  load.rem  :  std_logic  :=  '0f; 
signal  load_key  :  std_logic  :=  'O'; 
signal  new_rems  :  std_logic  :=  'O'; 
signal  new_key  :  std_logic  :=  'O'; 

signal  input  :  std_logic_vect or (31  downto  0)  :=  (others  =>  '0'); 

signal  key  :  st d_logic_vect or ( 31  downto  0)  :=  (others  =>  ' 0 ' ) ; 

signal  lower_bw_mask  :  st d_logic_vect or  ( 31  downto  0); 

signal  upper_bw_mask  :  st d_logic_vect or  ( 31  downto  0); 

signal  start  :  std_logic  :=  'O'; 

signal  next_hop  :  std_logic  :=  'O'; 

signal  hopset_ready  :  std_logic; 

signal  finished  :  std_logic; 

signal  hop.number  :  std_logic_vect or  ( 10  downto  0); 
signal  hop_channel  :  st d_logic_vect or  ( 1 0  downto  0); 


component  adapt ive_hopset_selector_synthesis  is 
port  ( 

signal  elk,  reset  :  in  std_logic; 
signal  load_rem  :  in  std_logic; 
signal  load_key  :  in  std_logic; 
signal  new.rems  :  in  std.logic; 
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signal  new.key  :  in  stcLlogic; 

signal  input  :  in  std_logic_vect or ( 31  downto  0); 

signal  key  :  in  std_logic_vect or (31  downto  0); 

signal  lower_bw_mask  :  in  std_logic_vect or  ( 31  downto  0); 

signal  upper_bw_mask  :  in  std_logic_vect or  ( 31  downto  0); 

signal  start  :  in  std.logic; 

signal  next_hop  :  in  std.logic; 

signal  hopset_ready  :  out  std_logic; 

signal  finished  :  out  std_logic; 

signal  hop_number  :  out  std_logic_vector  ( 10  downto  0); 
signal  hop_channel  :  out  st d_logic_vector  ( 1 0  downto  0) 

) ; 

end  component; 

—  function  st r_to_stdvec  ( input :  string)  return  std_logic_vect or  is 

—  variable  temp:  std_logic_vect or ( input ' range)  :=  (others  =>  ... 

’O'); 

—  begin 

—  for  i  in  input ' range  loop 

—  if (input (i)  =  '1')  then 

—  temp (i)  : =  ' 1 ' ; 

—  elsif (input (i)  =  '0')  then 

—  temp (i)  : =  'O'; 

—  end  if; 

—  end  loop; 

—  return  temp; 

—  end  function  str_to_stdvec; 

begin 

selector :  adaptive.hopset.selector.synthesis 
port  map ( 
elk  =>  elk, 
reset  =>  reset, 
load.rem  =>  load_rem, 
load.key  =>  load.key, 
new.rems  =>  new.rems, 
new.key  =>  new.key, 
input  =>  input, 
key  =>  key, 

lower_bw_mask  =>  lower.bw.mask , 
upper _bw_mask  =>  upper.bw.mask , 
start  =>  start, 
next.hop  =>  next.hop, 
hopset.ready  =>  hopset.ready , 
finished  =>  finished, 
hop.number  =>  hop.number, 
hop.channel  =>  hop.channel 


test:  process 

file  infile  :  text; 
file  outfile  :  text; 
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88 

—  variable  inline  :  line; 

89 

variable  outline  :  line; 

90 

—  variable  indata  :  string (64  downto  1); 

91 

variable  outdata  :  string (64  downto  1); 

92 

variable  maps  :  integer  :=  3; 

93 

94 

variable  seedl,  seed2  :  positive; 

—  Seed 

values  for  random  generator 

95 

variable  rand  :  real; 

—  ... 

Random  real— number  value  in  range  0  to 

1 . 0 

96 

variable  int_rand  :  integer; 

—  ... 

Random  integer  value  in  range  0..4095 

97 

variable  randVector  :  std_logic_vector (31 

downto  0 ) ;  —  ... 

Random  input  vector 

98 

variable  index  :  integer; 

99 

begin 

100 

—  Reset  the  system 

101 

reset  <  ' 1 ' ; 

102 

wait  for  period; 

103 

reset  <  'O'; 

104 

105 

—  Toggle  new  rems 

106 

new.rems  <  ' 1 '  ; 

107 

wait  for  period; 

108 

new_rems  <  'O'; 

109 

110 

—  Toggle  new  keys 

111 

new_key  <  '  1 '  ; 

112 

wait  for  period; 

113 

new_key  <  'O'; 

114 

115 

—  Load  the  REMs 

116 

for  m  in  1  to  maps  loop 

117 

for  i  in  0  to  63  loop 

118 

UNIFORM ( seedl ,  seed2,  rand);  . . . 

generate  random  number 

119 

int_rand  :=  INTEGER (TRUNC (rand*42 94 9672 95 . 0 )) ;  ... 

—  rescale  to  0..2~32-l,  .. 

find  integer  part 

120 

randVector  :=  std_logic_vector (to.unsigned ( int_rand,  ... 

randVector ' LENGTH) ) ;  —  convert  to 

std_logic_vector 

121 

122 

input  <  randVector; 

123 

— input  <  x"FFFFFFFF" ; 

124 

125 

load_rem  <  '  1 '  ; 

126 

wait  for  period; 

127 

load.rem  <  'O'; 

128 

wait  for  period; 

129 

end  loop; 

130 

end  loop; 
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—  Load  the  key 

for  i  in  0  to  63  loop 
index  :=  i  *  32; 

UNIFORM ( seedl ,  seed2,  rand) ;  ... 

—  generate 

random  number 

int.rand  :=  INTEGER (TRUNC (rand* 4 2  94  9  67295.0) )  ;  ... 

—  rescale  to  0..2~32— 1,  find 

integer  part 

randVector  :=  std_logic_vector (to_unsigned ( int.rand,  ... 
randVector ' LENGTH) ) ;  —  convert  to  std_logic_vector 

key  <  randVector; 

load-key  <  '  1 '  ; 
wait  for  period; 
load.key  <  'O'; 
wait  for  period; 
end  loop; 

—  Set  half— bandwidth 
lower_bw_mask  <  x" 80000000"; 
upper_bw_mask  <  x " 00000001" ; 

—  Toggle  start  flag 
start  <  '  1 ' ; 

wait  for  10*period; 
start  <  '0f; 

—  Wait  until  the  hopset  has  been  calculated 
wait  until  hopset.ready  =  '1'; 

—  Wait  a  clock  cycle  for  good  measure 
wait  for  15*period; 

—  Open  the  file  for  reading 

—  file.open (outfile,  (  " . \hopsets\hopset "  &  integer ' image (i)  & 

" . t xt " )  ,  WRITE  _MODE )  ; 

file_open (outfile,  (".\hopsets\hopset.txt"),  WRITE.MODE) ; 

—  Retrieve  all  hopsets  in  order 
for  i  in  0  to  2047  loop 

next.hop  <  ' 1 ' ; 
wait  for  5*period; 
next-hop  <  'O'; 
wait  for  5*period; 

write (outline,  .  .  . 

integer  '  image  ( t  o.integer  (unsigned  (hop.channel )  )  )  )  ; 
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176  writeline  (outfile,  outline); 

177  end  loop; 

178 

179  — report  "Wrote  hopset  "  &  integer '  image  (i)  &  "  to  file"; 

iso  report  "Wrote  hopset  to  file"; 

181 

182  file_close (outfile)  ; 

183 

184  wait  for  10*period; 

185 

186  new.rems  <  '  1 '  ; 

187  wait  for  5*period; 

188  new_rems  <  'O'; 

189 

190  new_key  <  '  1 '  ; 

191  wait  for  5*period; 

192  new.key  <  'O'; 

193 

194  —  End  test 

195  wait; 

196  end  process  test; 

197 

198  clock:  process 

199  begin 

200  elk  <  '  1 '  ; 

201  wait  for  period  /  2; 

202  elk  <  'O'; 

203  wait  for  period  /  2; 

204  end  process  clock; 

205  end  behavioral; 


1  library  ieee; 

2  use  ieee . std_logic_l 1 64 . all ; 

3  use  ieee .math_real . all ; 

4  use  ieee . numeric_std . all ; 

5 

6  entity  tb_hopset_selector  is 

7  generic ( 

8  period:  time  :=  10  ns 

9  )  ; 

10  end  entity; 

11 

12  architecture  behavioral  of  tb_hopset_selector  is 

13  signal  elk,  reset  :  std_logic; 

14  signal  start  :  std_logic  :=  'O'; 

15  signal  aggregate_rem  :  std_logic_vector  (2047  downto  0)  :=  ... 

(others  =>  ' 0 ' ) ; 

16  signal  key  :  std_logic_vector  (2047  downto  0)  :=  (others  =>  '0'); 

n  signal  next_hop  :  std_logic  :=  'O'; 

is  signal  hopset_ready  :  std_logic  :=  'O'; 

19  signal  finished  :  std_logic  :=  'O'; 


—  for  UNIFORM,  TRUNC  functions 

—  for  TO-UNSIGNED  function 
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47 
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signal  hop.number  :  std_logic_vector ( 10  downto  0)  := 

'  0  ' )  ; 

signal  hop_channel  :  std_logic_vector  ( 10  downto  0)  :=  (others  .. 

=>  ' 0 ' ) ; 

component  hopset_selector_b  is 
port  ( 

signal  elk,  reset  :  in  std_logic; 
signal  start  :  in  std_logic; 

signal  aggregate.rem  :  in  std_logic_vector (2047  downto  0); 

signal  key  :  in  std_logic_vector  (2047  downto  0); 

signal  next_hop  :  in  std_logic; 

signal  hopset.ready  :  out  std_logic; 

signal  finished  :  out  std_logic; 

signal  hop_number  :  out  std_logic_vector  ( 1 0  downto  0) ; 
signal  hop_channel  :  out  std_logic_vector ( 1 0  downto  0) 

)  ; 

end  component; 
begin 

selector:  hopset_selector_b 
port  map ( 
elk  =>  elk, 
reset  =>  reset, 
start  =>  start, 

aggregate_rem  =>  aggregate_rem, 
key  =>  key, 
next_hop  =>  next_hop, 
hopset.ready  =>  hopset.ready , 
finished  =>  finished, 
hop_number  =>  hopmumber, 
hop_channel  =>  hop_channel 

) ; 

test:  process 

variable  seedl,  seed2  :  positive;  —  Seed 

values  for  random  generator 

variable  rand  :  real;  —  .  .  . 

Random  real— number  value  in  range  0  to  1.0 
variable  int_rand  :  integer;  —  . . . 

Random  integer  value  in  range  0..4095 
variable  randVector  :  std_logic_vector  (31  downto  0);  —  ... 

Random  input  vector 
variable  index  :  integer; 
begin 

—  Reset  the  system 
reset  <  '  1 ' ; 
wait  for  period; 
reset  <  'O'; 

—  Load  the  aggregate  REM 
for  i  in  0  to  63  loop 
index  :=  i  *  32; 
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UNIFORM ( seedl ,  seed2,  rand) ;  ... 

—  generate 

random  number 

int.rand  :=  INTEGER (TRUNC (rand*4294967295 . 0)  )  ;  ... 

—  rescale  to  0..2"32— 1,  find 

integer  part 

randVector  :=  std_logic_vector (to_unsigned ( int.rand,  ... 
randVector ' LENGTH) ) ;  —  convert  to  std.logic.vector 

aggregate_rem (index+31  downto  index)  <  not  randVector; 

wait  for  period; 
end  loop; 


—  Load  the  key 
for  i  in  0  to  63  loop 
index  :=  i  *  32; 

UNIFORM ( seedl ,  seed2,  rand);  ... 

—  generate 

random  number 

int.rand  :=  INTEGER (TRUNC (rand*4294967295 . 0) )  ;  ... 

—  rescale  to  0..2~32— 1,  find 

integer  part 

randVector  :=  std.logic.vector  (to.unsigned  ( int.rand,  ... 
randVector ' LENGTH) ) ;  —  convert  to  std.logic.vector 

key (index+31  downto  index)  <  randVector; 


87  wait  for  period; 

88  end  loop; 

89 

90  start  <  '  1 '  ; 

91  wait  for  2*period; 

92  start  <  'O'; 

93 

94  wait  until  hopset.ready  =  '  1 ' ; 

95 

96  wait  for  5  *  period; 

97 

98  for  i  in  0  to  2047  loop 

99  next.hop  <  '  1 '  ; 

100  wait  for  period; 

101  next.hop  <  'O'; 

102  wait  for  period; 

103  end  loop; 

104 

105  —  End  test 

106  wait; 

107  end  process  test; 
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108 

109 

clock : 

process 

110 

begin 

111 

elk 

<  '  1 '  ; 

112 

wait 

for  period 

/  2; 

113 

elk 

<  'O'; 

114 

wait 

for  period 

/  2; 

115 

end  process  clock; 

1 16  end  behavioral; 


1  library  ieee; 

2  use  ieee . std_logic_l 1 64 . all ; 

3 

4  entity  tb_rem_merger  is 

5  generic ( 

6  period:  time  :=  50  ns 

7  )  ; 

8  end  entity; 

9 

10  architecture  behavioral  of  tb_rem_merger  is 

11  —  Component  signals 

12  signal  elk,  reset  :  std_logic; 

13  signal  load  :  std_logic  :=  'O'; 

14  signal  newunaps  :  std_logic  :=  'O'; 

15  signal  input_upper,  input_lower  :  std_logic_vector  (31  downto  0)  ... 

: =  (others  =>  ' 0 ' ) ; 

16  signal  output_upper ,  output.lower  :  std_logic_vector  (31  downto  .  .  . 

0)  :=  (others  =>  '0'); 

17 

is  —  Internal  signals 

19  signal  output  :  std_logic_vector  ( 63  downto  0); 

20 

21  component  remjnerger  is 

22  port  ( 

23  signal  elk,  reset  :  in  std_logic; 

24  signal  load  :  in  std_logic; 

25  signal  new_maps  :  in  std_logic; 

26  signal  input_upper,  input_lower  :  in  std_logic_vector  (31  ... 

downto  0 )  ; 

27  signal  output-upper ,  output_lower  :  out  std_logic_vector  (31  .  .  . 

downto  0) 

28  )  ; 

29  end  component; 

30  begin 

31  merger:  remjnerger 

32  port  map  ( 

33  elk  =>  elk, 

34  reset  =>  reset, 

35  load  =>  load, 

36  newunaps  =>  newunaps, 

37  input_upper  =>  input-upper. 
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38  input_lower  =>  input_lower, 

39  output_upper  =>  output_upper , 

40  output_lower  =>  output_lower 

41  )  ; 

42 

43  output  <  output_upper  &  output.lower ; 

44 

45  test:  process 

46  begin 

47  —  Reset  the  system 

48  reset  <  '  1 ' ; 

49  wait  for  period; 

so  reset  <  'O'; 

51 

52  —  Test  #1 

53  input-upper  <  x"FFFF0000"; 

54  input_lower  <  x"0000FF00"; 

55  load  <  '  1 '  ; 

56  wait  for  period; 

57  load  <  'O'; 

58 

59  wait  for  period; 

60 

61  —  Check 

62  assert  (output  =  x"FFFF00000000FF00" ) 

63  report  "Test  #1  failed!" 

64  severity  error; 

65 

66  —  Test  #2 

67  input-upper  <  x"0FF00000"; 

68  input_lower  <  x"FFFFFFFF"; 

69  load  <  '  1  '  ; 

70  wait  for  period; 

71  load  <  'O'; 

72 

73  wait  for  period; 

74 

75  —  Check 

76  assert  (output  =  x"  0FF000000000FF00 "  ) 

77  report  "Test  #2  failed!" 

78  severity  error; 

79 

so  —  Test  #3 

si  input_upper  <  x"0FFF0000"; 

82  input_lower  <  x"0000FFFF"; 

83  load  <  '  1 '  ; 

84  wait  for  period; 

85  load  <  'O'; 

86 

87  wait  for  period; 

88 

89  —  Check 
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90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 


assert  (output  =  x"0FF000000000FF00" ) 
report  "Test  #3  failed!" 
severity  error; 

—  Test  #4  (Soft-reset  for  existing  output) 
newjnaps  <  '  1 '  ; 
wait  for  period; 
newjnaps  <  'O'; 

wait  for  period; 

—  Check 

assert  (output  =  x"FFFFFFFFFFFFFFFF" ) 
report  "Test  #4  failed!" 
severity  error; 

—  End  test 
wait  ; 

end  process  test; 

clock:  process 
begin 

elk  <  ' 1 ' ; 

wait  for  period  /  2; 
elk  <  'O'; 

wait  for  period  /  2; 
end  process  clock; 
end  behavioral; 
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